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Abstract 

We give an overview of the Integer Quantum Hall Effect. We propose a mathematical 
framework using Non-Commutative Geometry as defined by A. Connes. Within this 
framework, it is proved that the Hall conductivity is quantized and that plateaux occur 
when the Fermi energy varies in a region of localized states. 

PACS numbers: 73.40.Hm, 02.30.Wd, 72.15.Rn 

This article will appear in the October 94 special issue of the Journal of Mathematical 
Physics. 



Contents 



1 Introduction 



2 IQHE: experiments and theories 



2.1 The classical Hall effect 



4 

4 



2.2 The Quantum Hall Effect) 6 



2.3 The Hall effect for the free Fermi gas 



2~4 'i'he role of localization! 10 



2.5 The Laughlin argumenti 12 

2.6 'I'he (Jhern-Kubo relationship] 13 



3 The Non-Commutative Geometry of the IQHE 



16 

O The non-commutative Brillouin zonel 17 



3.2 



Hall conductance and non-commutative Chern character] 19 

20 



3.3 Localization and plateaux of the Hall conductance 



3.4 Summary of the main result^ 22 

3.5 Homogeneous Schrodinger's operators] 23 



3.6 Observables and calculus 25 



*e-niail: jeanbel@siberia.ups-tlse.fr 
t e-mail: andreas@tolosa.ups-tlse.fr 
■I- e-mail: hermann@tolosa.ups-tlse.fr 

§URA 505, CNRS and Laboratoire de Physique Quantique, 118, Route de Narbonne, 31062-Toulouse Cedex, 
France 



4 The Kubo-Chern formula 



4.1 



30 

30 

Kubo's tormulal 34 



The relaxation time approximation 



4.3 Estimating the deviations from the IQHE limit| 36 



4.4 Dixmier trace and Sobolev space| 39 



15 Non-commutative Chern character! 44 



4.6 Connes formula 46 



O Chern character and I'redholm indexl 50 



4.8 Quantization, Fredholm and relative index] 51 



5 Localization and non-commutative Sobolev space 



52 

5.1 'I'he Anderson-Pastur localization criterion] 53 

5.2 Non-commutative localization criterion and localization length] 55 



5.3 Localization in physical models] 62 



6 



64 



Applications and complements] 

]6.1 Low-Iyiug states do not contribute to the IQHE 64 



]6.2 Where and how does the localization length diverge?] 65 



]6.3 Chern numbers and localization in Harper's equation 



7 Introduction to the FQHE 



fn~a 



verview 



7.2 Laughlin's ansatz for the u 



O 'I'he elementary theory of the v 



66 

69 

69 

groundstatel 70 

= i:-FQHE] 75 



7.4 The role of gauge invariance and incompressibility] 77 



1 Introduction 

In 1880, E.H. Hall [|T| undertook the classical experiment which led to the so-called Hall effect. 
A century later, von Klitzing and his co-workers [Q] showed that the Hall conductivity was 
quantized at very low temperatures as an integer multiple of the universal constant e'^/h. Here 
e is the electron charge whereas h is Planck's constant. This is the Integer Quantum Hall Effect 
(IQHE). For this discovery, which led to a new accurate measurement of the fine structure 
constant and a new definition of the standard of resistance , von Klitzing was awarded the 
Nobel price in 1985. 

On the other hand, during the seventies, A. Connes 0, |^ extended most of the tools of 
differential geometry to non-commutative C*-algebras, thus creating a new branch of math- 
ematics called Non-Commutative Geom,etry. The main new result obtained in this field was 
the definition of cyclic cohomology and the proof of an index theorem for elliptic operators 
on a foliated manifold. For this work and also his contribution to the study of von Neumann 
algebras, Connes was awarded the Fields Medal in 1982. He recently extended this theory to 
what is now called Quantum Calculus 0. 

After the works by Laughlin and especially by Kohmoto, den Nijs, Nightingale and 
Thouless [H (called TKN2 below), it became clear that the quantization of the Hall conduc- 
tance at low temperature had a geometric origin. The universality of this effect had then an 



explanation. Moreover, as proposed by Prange [^ |T^ and Thouless [|TI|, the plateaux of the 
Hall conductance which appear while changing the magnetic field or the charge-carrier density, 
are due to localization. Neither the original Laughlin paper nor the TKN2 one however could 
give a description of both properties in the same model. Developing a mathematical framework 
able to reconcile topological and localization properties at once was a challenging problem. At- 



tempts were made by Avron et al. [0 who exhibited quantization but were not able to prove 



that these quantum numbers were insensitive to disorder. In 1986, H. Kunz |]T3| went further on 
and managed to prove this for disorder small enough to avoid filling the gaps between Landau 
levels. 



But in [|14|, |T^, 0, one of us proposed to use Non-Commutative Geometry to extend the 
TKN2 argument to the case of arbitrary magnetic field and disordered crystal. It turned out 
that the condition under which plateaux occur was precisely the finiteness of the localization 



length near the Fermi level. This work was rephrased later on by Avron et al. |T^ in terms 
of charge transport and relative index, filling the remaining gap between experimental obser- 
vations, theoretical intuition and mathematical frame. 

Our aim in this work is to review these various contributions in a synthetic and detailed 
way. We will use this opportunity to give proofs that are missing or scattered in the literature. 
In addition, we will discuss the effect of disorder from two complementary aspects. 

On the one hand, we will develop our point of view on localization produced by quenched 
disorder. This is crucial for understanding the IQHE. We review various localization criteria 
and formulate them in terms of Non-Commutative Geometry. With the Dixmier trace, A. 
Connes introduced a remarkable technique into Quantum Calculus. In our context, it allows 
us to give the precise condition under which the Hall conductance is quantized; this condition 
is shown to be a localization condition. 

On the other hand, we also propose a model for electronic transport giving rise to the 
so-called "relaxation time approximation" and allowing to derive a Kubo formula for the con- 
ductivity. This approach allows to describe the effect of time-dependent disorder in a phe- 
nomenological way. This latter has quite different consequences from those of the quenched 
disorder such as a non-zero finite direct conductivity. Even though this approach is not orig- 
inal in its principle, the non-commutative framework allows us to treat the case of aperiodic 
crystals and magnetic fields when Bloch theory fails. Therefore, strictly speaking, our Kubo 
formula is new. We also show, without proofs, how to justify the linear response theory within 
this framework, leaving the formal proofs for a future work. The advantage of this approach 
is to give control of the various approximations that have to be made to fit the ideal result 
with experiments. For this reason, we discuss the effects of temperature, of non-linear terms 
in the electric field, of the finite size of samples and finally those of collisions and disorder. In 
particular, we argue that the discrepancy dan between the measured Hall conductivity and the 
ideal one, given by a Chern number, is dominated by the collision terms. In the center of a 
plateau, we get the rough estimate 

< const. V- , (11 

an huc 

where z/ is the filling factor, A is the localization length and fic is the charge- carrier mobility. 
e/h is a universal constant, u about unity and the localization length typically of the order 
of the magnetic length. Inserting measured values for the mobility, one obtains 10~^ for the 
right-hand side expression. This estimate does not take into account the Mott conductivity. 
However, it shows why both a large quenched disorder (in order to have small localization 



lengths) and a large mobility (namely a low collision rate) are necessary in order to get ac- 
curate measurements. Such a compromise is realized in hetero junctions and to less extend 
in MOSFETs. The estimate (0) also permits to understand intuitively why the plateaux in 
the Fractional Quantum Hall Effect (FQHE) are less precise, since the localization length of 
Laughlin quasiparticles is probably larger than that of electrons at integer plateaux, and their 
mobility probably lower. 

No attempt will be made however to extend our non-commutative approach to the FQHE. 
We will only give some insight and a short review of works that we feel relevant in view of a 
mathematically complete description of the FQHE. 

This rest of the article is organized as follows. In Chapter ^ we give the conventional 
explanations of the IQHE. In particular, we discuss the Laughlin argument, the topological 
aspect introduced by TKN2 and the effects of localization in a qualitative way. Chapter ^ is 
devoted to the mathematical framework needed for Non-Commutative Geometry. In particular 
we describe how to overcome the difficulty of not having Bloch's theorem for aperiodic media. 
We then show that the Brillouin zone still exists as a non-commutative manifold. We also give 
the main steps of our strategy leading to a complete mathematical description of the IQHE. In 
Chapter | we discuss transport theory leading to Kubo's formula. We show that in the IQHE 
idealization, the Hall conductance is a non-commutative Chern number. We also relate this 
Chern number to a Fredholm index which leads to the quantization of the Hall conductance. 
Through the notion of relative index we show in which sense this approach is a rigorous version 
of the Laughlin argument. Chapter ^is devoted to localization theory. We give various criteria 
and define various localization lengths which are commonly used in the literature. We also show 
how to express these notions in the non-commutative language. This part allows us to explain 
on a rigorous basis the occurrence of plateaux of the Hall conductance. Finally, we show that 
such criteria are in fact satisfied in models such as the Anderson model. In Chapter ^ we give 
some consequences of this theory for practical models. In particular we show that low-lying 
states do not contribute to the IQHE. We also discuss the open question where the jumps of 
the Hall conductance occur. The last Chapter |^ is a short review of available results on the 
FQHE. 
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2 IQHE: experiments and theories 

2.1 The classical Hall effect 

Let us consider a very fiat conductor placed in a constant uniform magnetic field in the z 
direction perpendicular to the plane Oxy of the plate (see Fig.l). If we force a constant current 
in the x direction, the electron fiuid will be submitted to the Lorentz force perpendicular to the 




Figure 1: The classical Hall effect: the sample is a thin metallic plate of width 6. The magnetic 
field B is uniform, and perpendicular to the plate. The current density j parallel to the x-axis 
is stationary. The magnetic field pushes the charges as indicated creating the electric field £ 
along the y direction. The Hall voltage is measured between opposite sides along the y-axis 

current and the magnetic field. Hall realized that the electron fluid is incompressible so that 
the Lorentz force must produce a pressure, namely a potential difference perpendicular to the 
current. 

Let j be the current density, B the magnetic fleld and S the Hall electric field. In a stationary 
state, the electric forces acting on the charges are opposite to the Lorentz forces. This leads to 
the equation 



nq£ + j X B = 



(2) 



where n is the charge- carrier density and q is the charge of the carriers. Since the magnetic 
field B is perpendicular to both j and S, solving @ for j gives 



J 



^BxS 



a-S 



where B is the modulus of the magnetic field and cr is the conductivity tensor. The anti-diagonal 
components of the tensor are the only non-vanishing ones and can be written as ±0"/^ 5, where 
6 is the plate width and an is called the Hall conductance. Thus 



C^H 



qn6 



We remark that the sign of an depends upon the sign of the carrier charge. In particular, the 
orientation of the Hall field will change when passing from electrons to holes. Both possibil- 
ities were already observed by Hall using various metals. This observation is commonly used 
nowadays to determine which kind of particles carries the current. 

Let i be the plate width in the y direction (see Fig.l). The current intensity inside the plate 
is then given by / = j6£ where j is the modulus of j. The potential difference created by the 
Hall field is Vh = —iS.u if u is the unit vector along the y axis. Using (|^) we find: 



V, 



H 



BI 

nq6 



I 

cth 



(3) 



In particular, for a given current intensity /, the thinner the plate the higher the potential 
difference. For example, for a good conductor like gold at room temperature, the charge carrier 




Figure 2: Schematic representation of the experimental observations in the IQHE. The Hall 
conductivity an is drawn in units of e^/h versus filling factor v. The dashed line shows the 
Hall conductivity of the Landau Hamiltonian without disorder. The direct conductivity a// is 
shown in arbitrary units. 



see |T^ chap. 1). Thus, for a magnetic field of IT, a current 



density is of order of 6 x lO'^^m' 
intensity of lA and a potential difference of ImV the plate width is about Ifim. These numbers 
explain why the effect was so difficult to observe. It forced Hall to use very thin gold leaves 
in the beginning. In modern devices, much thinner "plates" with thickness of about lOOA are 
produced in inversion layers between two semi-conductors. 

In view of @, the Hall conductance has the dimension of the inverse of a resistance. Since 
the product n6 is the number of charge carriers per unit area, the dimensionless ratio 

nSh 



Be 



called the filling factor, represents the fraction of a Landau level filled by conduction electrons 
of the thin plate. In terms of this parameter, we obtain for a free electron gas: 



(^H 



(4) 



Rh e^ 

where Rh is called the Hall resistance. It is a universal constant with value Rh = 25812. 80O. 
Rh can be measured directly with an accuracy better than 10~^ in QHE experiments. Since 
January 1990, this is the new standard of resistance at the national bureau of standards 0. 



2.2 The Quantum Hall Effect 

Let us concentrate upon the dependence of the Hall conductance (in units of e^ /h) on the 
filling factor v. In the classical Hall effect, these two quantities are just equal (eq.^). Lowering 



the temperature below IK leads to the observation of plateaux for integer values of the Hall 
conductance (see Fig.2). In von Klitzing's experiment 0, the variation of u was obtained by 
changing the charge carrier density, whereas in later experiments one preferred varying the 
magnetic field. The accuracy of the Hall conductance on the plateaux is better than lO"*^. For 
values of the filling factor corresponding to the plateaux, the direct conductivity a//, namely 
the conductivity along the current density axis, vanishes. These two observations are actually 
the most important ones. The main problems to be explained are the following: 

(i) Why do the plateaux appear exactly at integer values ? 

(ii) How do the plateaux appear ? 

(ill) Why are these plateaux related to the vanishing of the direct conductivity ? 

To observe the QHE, physicists have used conduction electrons trapped in the vicinity of an 
interface between two semiconductors. The local potential difference between the two sides 
produces a bending of the local Fermi level. Near the interface, this Fermi level meets with the 
valence band creating states liable to participate in the conductivity. This bending occurs on 
a distance of the order of lOOA from the interface, so that the charge carriers are effectively 
concentrated within such a thin strip. In addition, by changing the potential difference between 
the two sides, the so-called gate voltage, one can control the charge- carrier density. 

The samples used in QHE experiments belong to two different categories. The first one 



is called MOSFET |H], for metal-oxide silicon field effect transistors. The interface separates 
doped silicon from silicon oxide. This device was common in the beginning of the eighties and 
was the one used in von Klitzing's experiment. However, the electron mobility is relatively low 
because the control of the flatness of the interface is difficult. 

The samples of the other category are heteroj unctions. The interface separates GaAs from 
an alloy oi Al^Gai^x^s. This kind of device nowadays makes available interfaces almost without 
any defects. Moreover, electrons therein have a high mobility. These devices are most commonly 
used in modern quantum Hall experiments. 

In both kinds of samples, there are many sources of defects producing microscopic disorder. 
The first comes from the doping ions. Even though they are usually far from the interface 
(about lOOOA), the long range Coulomb potential they produce is strong enough to infiuence 
the charges on the interface. It is not possible to control the position of these ions in the crystal. 
The second source of defects is the roughness of the interface. This is an important effect in 
MOSFET's, much less in heteroj unctions. In the latter the accuracy is better than one atomic 
layer in every lOOOA along the interface ^^. Finally, long range density modulation of the 



compounds may produce visible effects. This is the case especially for heterojunctions where 
the aluminium concentration may vary by a few percent on a scale of Ifim [ PH| . 

It is important to notice that the observation of plateaux supposes several conditions. 

(i) The effect is more easily seen if the electron fluid is concentrated in such a thin region that it 
can be considered as two-dimensional. In fact, owing to the trapping effect of the potential 
interface, the motion perpendicular to the interface is quantized. For good samples in 
high magnetic flelds, the energy difference between two corresponding eigenvalues is big 
compared to ksT (where T is the temperature and ks the Boltzmann constant), so that 
only the lowest such level has to be considered. Hence the problem becomes effectively 
two-dimensional. 



(ii) The plateaux disappear beyond a temperature of a few Kelvin. More exact, the inelastic 
relaxation time has to be large enough; otherwise corrections will be needed in formulae 
calculating the current; this will lead to the destruction of the plateaux. This is the reason 
why the IQHE is seen more easily in heterojunctions than in MOSFETs: the electron 
mobility in the former is higher than in the latter. 

(ill) We will see that some quenched disorder producing only elastic scattering is necessary for 
the appearance of the plateaux (see Section |2.4j). In practice, the disorder that occurs is 



strong enough to produce a filling of the gaps between Landau levels [21, 22 



(iv) Clearly the sample size must be big enough as to allow the use of the infinite volume limit. 
Mesoscopic systems exhibit conductance fluctuations from sample to sample which may 
partially distroy the effect. Finite volume effects however have been shown to decrease 
exponentially fast with the sample size. 

(v) Finally the electric field needed to produce the current has to be small. If it is too high, 
non-linear phenomena may distroy the plateaux ||23||. 



Provided the previous conditions are satisfied, the quantum Hall effect is a universal phe- 
nomenon. It is quite independent of the specific shape of the sample. Hall plateaux have 
also been observed in microwave experiments where the topology of the sample is trivial [^ 



Note, however, that the centers of the plateaux need not be located near integer values of the 
filling factor. It depends upon which kind of doping ion is used |]2S[. The values of the Hall 



conductance, however, on the plateaux are independent of the nature of the used sample. 

2.3 The Hall effect for the free Fermi gas 

Let us first consider a very simple model in which the charge carriers are spinless, free, two- 
dimensional fermions with charge q. Our aim is to show that no quantization of the Hall 
conductance is observed in this case. 

Since the particles are independent, the quantum motion is described by the one-particle 
Hamiltonian. Let A = (^41,^2) be the vector potential given by 

d,A2 - ^2^1 = B , 

where B is the modulus of the magnetic field. The energy operator is then given by the Landau 



Hamiltonian [EH] 



_ (P - qAf 
^' - 2m. ' ^^^ 

where P is the 2D momentum operator and m^ is the effective mass of the particle. This 
operator is not translation invariant owing to the symmetry breaking produced by the vector 
potential. However, if one replaces the usual representation of the translation group by the 
so-called magnetic translations ||2^, the Landau Hamiltonian becomes translation invariant. 
Let us introduce the quasimomentum operators K = [P — qA)/h; they fulfill the canonical 
commutation relations 



Here qB/h plays the role of an effective Planck constant. We see that the Landau Hamiltonian 
describes a harmonic oscillator so that its spectrum is given by the Landau levels, namely 

En = hujc{n + -) , 

where ujc = qB/m^, is the cyclotron frequency and n G N. Each Landau level is infinitely 
degenerate owing to translation invariance. The degeneracy per unit area is finite however and 
given by qB/h. Perturbing this operator will give rise to a band spectrum. 

Let us now compute the current. We assume the system to be in a thermodynamical equi- 
librium. In order to allow non-zero current, we describe the system in the grand-canonical 
ensemble; in practice such a fluid is open. Since the particles are independent fermions, the 
thermal averaged density per unit volume of a one-particle translation invariant extensive ob- 
servable O at temperature T and chemical potential /i is given by 

<0 >T,,= ^m^^Trj,{fTAHL)0) , 

where fT,^{E) = (1 + e^^^~^')^^ is the Fermi-distribution function and (3 = l/ksT. In the 
limit above, A denotes a square box centered at the origin, and |A| is its area. The current is 
represented by the operator 

where X is the position operator. It is easy to check that this current operator commutes 
with the magnetic translations. Obviously the thermal average of the current vanishes since 
no current can flow without an external source of energy. To make it non-zero, we switch 
on an electric fleld S at time t = 0. For simplicity, we assume this fleld to be uniform and 
time-independent. After the fleld has been switched on, the time evolution is given by 

dJ I ^ 

-^ = J- [Hl/^ J] ^ with H^,^ =Hl- q£X . 

The solution of this equation is elementary. Using complexifled variables, that is M = M1 + 1M2 
whenever M = (Mi, M2) is a 2D vector, we flnd 

Jit) = -iq^ + e-*-=* Jo , 

where Jo is some initial datum. This solution consists of a time-independent part and of an 
oscillating part with the cyclotron frequency. The time average is just the constant part; it is 
the system's response to the applied electric fleld. Taking the thermal average of this constant 
part, we have the observed current, namely, if j = ji + zJ2; 

This is nothing but the classical formula (0), written in complex notation. 

Therefore we see that the classical Hall formula still holds in quantum mechanics for the 
free fermion gas at all temperatures. There is no way to see any trace of quantization of the 
Hall conductance, neither is there any kind of plateaux of the Hall conductance! 



2.4 The role of localization 

The main experimental property we have not yet considered is the vanishing of the direct 
conductivity when the filhng factor corresponds to the plateaux. In the previous argument 
we ignored all effects leading to a finite non-zero direct conductivity. Among these effects 
there are several sources of dissipation such as phonon scattering or photoemission. At very 
low temperature, these sources usually have very limited influence and impurity scattering 
dominates. This is why we are led to consider non-dissipative effects like Anderson localization 
in disordered systems. As explained in Section p.2| , several kinds of defects can influence the 
electron motion. These defects are usually distributed in a random way so that the forces they 
create on the charge carriers are actually represented by a random potential. In many cases 
one considers these defects as isolated and of small influence, so that a first-order perturbation 
theory based upon one-electron scattering on one impurity already gives a good account of the 
observed effects. For 2D systems however, it turns out that the low-density limit for impurities 
does not give the relevant contribution. We will see in Section ^]5| how to define properly a 
potential representing the effects of a high density of random scatterers. For the moment let 
us stay at an intuitive level. 

As explained by P.W. Anderson |2^, the occurrence of a random potential in a one-particle 
Hamiltonian may lead to the quantum localization of particles. More precisely, the quantal 
wave representing these particles reflects on the potential bumps producing interferences. The 
Bragg condition is necessary for building a constructive interference pattern throughout the 
crystal. This requires some regularity of the crystal, such as periodicity or quasiperiodicity, 
in order to allow the wave to propagate. If the potential creating these bumps exhibits some 
randomness, the probability for the Bragg condition to be satisfled everywhere in the crystal 
eventually vanishes. Therefore, the wave function will vanish at inflnity leading to the trapping 
of particles in the local minima of the potential. This is Anderson localization. 



It has been argued ||29[ that one-dimensional systems of non-interacting particles in a 
disordered potential exhibit localization at any strength of disorder. Mathematicians have 
proved such a claim under relatively mild conditions on the randomness of the potential 
PP] , pi| , P^ P5| , 0. They have also proved localization in any dimension for strong disor- 



der |35, 36, 37, 38, B9l HO]. A flnite-size scaling argument, proposed by Thouless [41], has been 



used in [^ to show that two-dimensional systems are also localized at any disorder unless 



there are spin-orbit couplings [0. The same argument also shows that in higher dimensions, 
localization disappears at low disorder. We remark that the notion of localization we are using 
does not exclude divergence of the localization length at an isolated energy. These results were 
supplemented by many numerical calculations [^, |^, ^, |46[] . 

On the other hand, the occurrence of a random potential will create new states with energies 
in the gaps between Landau levels. This can be measured by the density of states (DOS). To 
deflne it, let M{E) be the integrated density of states (IDS), namely the number of eigenstates 
of the Hamiltonian per unit volume below the energy E: 

M{E) = lim ^#{eigenvalues of H\a < E}, (6) 

A--+00 ]A] 

where degenerate eigenvalues are counted with their multiplicity. This is a non-decreasing 
function of E. Therefore its derivative p{E) = dJ\f{E)/dE is well deflned as a Stieljes-Lebesgue 
positive measure and is called the DOS. Under mild conditions on the distribution of the random 



potential, it is possible to show ^, ^ |33[ that the DOS is a smooth function (see Figure ^^) . 



If the potential strength exceeds the energy difference huc between two consecutive Landau 
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Figure 3: The density of states for the Landau Haniiltonian with a random potential. 

levels, the gaps will be entirely filled. This is actually the physical situation in most samples 
used up to now. In heterojunctions, for instance, the minimum of the density of states may 
represent as much as 30% of its maximum l21f| , although in modern samples it is usually very 



small HI. 

Recall that the chemical potential at zero temperature is called the Fermi level Ep. Since 
the charge carriers are spinless fermions, their density at T = is given by 



n 



f ^ dE p{E) = M{Ef) 

J —oo 



(7) 



The absence of spectral gaps means that the IDS is monotone increasing, so that changing 
continuously the particle density is equivalent to changing continuously the Fermi level. More 
generally, if a magnetic field is switched on, changing continuously the filling factor is equivalent 
to changing continuously the Fermi level. On the other hand, while the Fermi level crosses a 
region of localized states, the direct conductivity must vanish whereas the Hall conductivity 
cannot change. This last fact is not immediately obvious, but more justification will be given 
later on. Changing the filling factor therefore will force the Fermi level to change continuously 
within this region of localized states while the Hall conductance will stay constant. This is the 
main mechanism leading to the existence of plateaux. 

In contrast, if the spectrum had no localized states, the Hall conductance would change 
while changing the filling factor, as long as the Fermi level would move within the spectrum. 
Moreover, if a spectral gap occurred between two bands, let us say between energies E^ and 
E^, the IDS would equal a constant uq on that gap and changing the filling factor from n^ — e 
to no + e would cause the Fermi level to jump discontinuously from E^ to E^ with the value 
{E^ + E^)/2 whenever n = uq. In this case, once again, no plateaux could be observed. This 
is why there is no quantized Hall effect in the free fermion theory. 

One of the consequences of this argument is that between two different plateaux, there must 
be an energy for which the localization length diverges |3^, |^. Even though there should be 
no extended states in 2D disordered systems, the localization length need not be constant in 
energy. As it happens, if the impurities are electrically neutral in the average, the localization 
length diverges exactly at the Landau levels. Therefore, the influence of impurities decreases 
near the band centers, and other sources of interactions, like the Coulomb potential between 
charge carriers, may become dominant. This is actually the basic observation leading to the 
understanding of the FQHE. 

We have argued that the direct conductivity should vanish whenever the Fermi level lies 
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Figure 4: The chosen geometry for the Laughhn Gedanken experiment. 

in a region of localized states. One may wonder how it is then possible to have a non-zero 
Hall conductance. The answer is that the Hall current is actually carried by edge states. This 
has been seen in numerical calculations [E3, M and there are the first theoretical hints [pOll. 



The experimental situation is not as clear as it looks from theory. This is probably due to the 
existence of zones in which states are localized surrounded by filamentary regions in which the 
Hall current can actually flow. 



2.5 The Laughlin argument 

The first attempt to explain the integer quantization of the Hall conductance in units of e'^/h 
was proposed by Laughlin 0. Laughlin originally chose a cylinder geometry for his argument. 
His justification for this was that the effect seemed to be universal and should therefore be 
independent of the choice of the geometry. Here, we present a Laughlin argument in the plane 



in form of a singular gauge transformation. In Section [4.8| , this presentation and an observation 
of Avron, Seller and Simon [|1^ will allow us to discuss the links between the Laughlin argument 
and the Chern-character approach presented in this article. 

Let us consider a free gas of non-interacting electrons in the two-dimensional plane subjected 
to an exterior magnetic field perpendicular to the plane. We choose an origin and then pass 
an infinitely thin flux tube through it. A radial electromotive force is created by means of a 
slowly varying flux (f){t) (see Fig |2.5|) . In polar coordinates, the vector potential with symmetric 
gauge for the constant field is A{r,9) = {—^rsinO — 2^sin6' 4>(t), ^rcos9 + ^.cos9 (pit)); the 
Hamiltonian then reads: 



Hr. 



1 



2m* 



-fi'd. 



2q2 



r \ I r I ZTir 



(8) 



Recall that the electron charge is — e. We now assume that the time evolution of the fiux 
through the cylinder is so slow that the adiabatic approximation describes the evolution of the 
states. The eigenstates may then be computed explicitly: 



Wn,m\^) "i t-J ^n,m,4> 



^-^me 



27re0(t) 



B Ln 



2we<p(t) 7^2 



(9) 
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where i^ = j2h/qB is the magnetic length; ra G N is the principal quantum number corre- 
sponding to the energy En = fiujciji + 1/2) and m G Z is the orbital quantum number so that 
the angular momentum is ml%\ Cn,m,<f> are easily computable normalization constants and the 
polynomials that appear are the Laguerre polynomials given by: 

L°(x) = ^e^a;-°a^(e-^x'^+°) n G N, a G R. 

As the flux (j){t) = — is varied from t = to t = r, one flux quantum is forced through the 
tube and a state ipn,m evolves to a state ipn,m+i up to a phase factor e~*^. Note that this reflects 
the fact that the singular gauge transformation introduces an extra angular momentum —h into 
the system. Now we assume that the filling factor is an integer N . The net effect of the above 
process is that the state with lowest angular momentum of each Landau level is transported to 
infinity. Let us fix a large circle C around the origin of radius R. The current density during 
the process is then (approximatively) given by j = 2^~T^'^ ^^^ strength of the electric field on 
the circle is S = —dt{(t>it) / {2ttR)) = —h/^reR); therefore, the Hall conductivity is calculated 
as an = ]/£■ = Ne'^/h. 

What did we learn from this argument ? Until now, not so much in fact. For had we taken 
another non-integer filling factor u and supposed that the electron density were uniform, we 
would obtain an = i^e'^/h; that is the classical result already calculated in Section ^]4[ 

Using ideas of Prange, Halperin and Joynt [|], |T0], ^ however, the picture may be completed 
to furnish a qualitative understanding of the IQHE in the following way: if the sample has 
impurities, these transform some of the extended states (^ into localized states; moreover, the 
explicit form of the remaining extended states changes. Now the localized states will not change 
as a flux quantum is forced through the flux tube, and it can be argued that the remaining 
extended states of one Landau level still carry the same current, that is one charge e through 
the circle C. As some of the localized states caused by the impurities are situated below the 
energy of the Landau levels and others above, we come to the desired qualitative explanation 
with help of the argument explained in Section p.4| . 

There are two ingredients in this Gedanken-experiment which ought to be emphasized. The 
first is the gauge invariance which produces the periodicity of the Hamiltonian with respect 
to the varying flux, the period being the flux quantum. The second is that the Hall effect 
corresponds to a charge transport of one unit for each filled Landau level. The conductance 
quantization is therefore likely to be connected to the charge quantization. 

2.6 The Chern-Kubo relationship 

The geometrical origin of the Hall conductance quantization was revealed by TKN2 P] and 



Avron et al. [|I^. TKN2 considered an electron gas submitted to a uniform magnetic field on 
a 2D square lattice in the tight-binding approximation. In the Landau gauge the Hamiltonian 
reads: 

{Hip){mi,m2) = ^{rrii -1,1712) + ^Ij {mi + I,m2) 

(10) 
+ e-2™"^i7/>(mi, ma - 1) + e2™™i?/'(mi, m2 + I) , 

where a = 0/0o is the ratio of the magnetic flux (p in the unit cell to the flux quantum 0o = h/e. 
Using the translation invariance along the y-axis, the solution of the eigenvalue equation can 
be written in the form 
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with 

ip{mi — 1) + f{rni + 1) + 2 cos (A;2 — 2TTami)(f{mi) = E{k2)ip{mi) . (11) 



This is Harper's equation |S^. It can also be obtained by adding a weak periodic potential to 
the Landau Hamiltonian (|^) and then projecting onto one Landau level. Now TKN2 assumed 
that the ratio a is a rational number of the form p/q where p and q are relatively prime integers. 
By translation invariance of period q, the solution of eq. (0) can also be written in Bloch's 
form 

so that Harper's equation now becomes 

{H{ki,k2)0{n) ■■= e*'=i^(n-l) + e-*'=i^(n + l) + 2cos(A;2-27r2n)^(n) 

= Eik,,k2)an). 

This is the secular equation for the q x q hermitian matrix H{ki,k2)- The energy spectrum 
is given by energy bands corresponding to the eigenvalues ei{ki, ^2), I < I < q. Let Pi{ki, /C2) 
be the corresponding eigenprojection. We will assume for simplicity that all bands are well 
separated: for the Harper equation, all gaps but the central one are open [^ Q ^ . 

It is quite clear that H{ki,k2) is a trigonometric polynomial in {ki,k2), implying that 
the eigenprojections and the eigenvalues are smooth periodic functions of (^1,^2). Moreover, 
introducing the two unitary q x q matrices defined by 

iui^)in) = ^{n - 1) , {U2i){n) = e"^*''i"^(n) , 

we derive the covariance relation 

u,H{ki, k2)u^^ = H(ki, k2 + 27r^) , U2H(ki, ^2)^^' = H{k^ - 27r^, ^2) . (12) 

q q 

These relations show that the eigenvalues and the trace of any function of H and its derivatives 
are actually 27r/g-periodic in (/ci, /C2). 



The Hall conductance can now be calculated by Kubo's formula |£6[. It gives the current as 
the velocity- velocity correlation function. The formula may be deduced from classical arguments 
using the Boltzmann equation |^, but a deduction closer to quantum mechanics using linear 
response theory will be presented in Chapter ^. The result derived by TKN2 is then the 
following: at zero temperature, if the Fermi level belongs to a gap of the Harper Hamiltonian 
of eq. ([XOp, the transverse conductivity is given by the Kubo-Chern relation (see |T 



an = ^Ch(Pp) 
h 



where 



Ch{PF) = j^^'^A2m-^l{PF{k^,k2)[d^PF{k^,k2),d2PF{Kk2)])\ . (13) 



/t2 47r2 [ q 
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Here da = d/dka,a = 1,2 and Pp is the eigenprojection on energies smaller than the Fermi 
level. The integrand is a complex 2-form on the torus; its de Rham-cohomology class is called 
the first Chern class or the first Chern character of the projection Pp. We shall also refer to 
the integral over this form as Chern character, because the integral will retain meaning in the 
non- commutative context; the integral will also be called Chern number. 

For further informations on Chern classes and characters we refer to ||57| , ^ . The Chern 



character of a projection is a topological invariant. To see this, we remark that the data of 
the family k = (fci,/c2) G T^ i-^ -Pf(^) of projections defines a complex fiber bundle Vp over 
the 2-torus in a natural way. More precisely, Vp is the set of pairs {k, ^) G T^ x C*^ such that 
-Pf(^)^ = ^- Since the Fermi level lies in a gap, Pp is smooth and its dimension is constant. 
Therefore Vp is a smooth vector bundle over the 2-torus. It can be shown that Ch(Pp) depends 
only upon the equivalence class of Vp provided the equivalence is isomorphism of vector bundles 
modulo adding trivial bundles. In particular, it is a homotopy invariant quantity. This implies 
that adding a small perturbation to the Harper Hamiltonian will not change the value of 
Ch(Pp), at least as long as the Fermi level does not cross the spectrum while turning on the 
perturbation! This is the argument that was needed to explain the robustness of the Hall 
conductance. 

In addition to being robust, the Chern number Ch(Pp) is actually an integer. To see this 
explicitly, we will show (see Section |4.5| ) that the Chern character Ch is additive with respect 
to the direct sum of two orthogonal projections P and Q, that is Ch(P©(5) = Ch(P) + Ch.{Q). 
Even though this is not obvious from eq. ([T3|) , all cross-terms in the left hand side vanish. It 
is thus sufficient to show that if P is a one dimensional smooth projection over T^, then its 
Chern character is an integer. This can be shown as follows. 

Let C,{k) be a unit vector in the vector space C^ such that P{k)C,{k) = C,{k), V/c. It is easy 
to show that 

Ch(P) = / ^TiiP{k)[diP{k),d2P{k)]) = - r dk^ r dk2'^m < dii{k)\d2m > , 

JT2 ZlTT TT Jo Jo 

(14) 
where Qm is the imaginary part. Using Stokes formula, this double integral can be written as 

ch(p) = l'^ ^ < m\dim > isr -r^< ^md,m > i^-r . (15) 

Even though P is doubly periodic in k, ^ need not be periodic. The obstruction to the period- 
icity of C, is precisely the non-vanishing of the Chern character. However, we can find two real 
functions 9i and 62 such that 

aku 27r) = e^'^^'-^^ah, 0) , e(2vr, ^2) = e^'^^'^^ao, h) , (16) 

provided < ki < 27t and < k2 < 27t respectively. The reason is that after one period the 
new vector ^ defines the same subspace as the old one, so they are proportional. Since ^{k) is 
a unit vector for any k, the proportionality factor must be a phase. Writing ^(27r, 27r) in two 
ways we get ^i(27r) — ^i(O) = 02(2vr) — 62(0) (mod 27r). Replacing in eq. (|l^) leads to 

Cj^.p. ^ (gi(27r)-gi(0))-(g2(27r)-g2(0)) ^ 
^ ^ 27r 

Thanks to this argument, we see that the Chern character of a smooth projection valued 
function over the 2-torus is an integer. 
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We remark, however, that the expression used in eq. (0) is not exactly the same as the one 
used in eq. (|TB|): in the latter, we have divided by q. However, the covariance relation eq. (0) 
shows that the integrand in eq. (0) is 27r/g-periodic in k so that it is sufficient to integrate 
over the square [0,27r/g]^^ and multiply the result by g^. On the other hand, the periodicity 
condition (^6[) has to be modified to take the covariance into account. Now the argument goes 
along the same line and gives rise to an integer as well. 

This series of arguments shows that indeed the Hall conductance (in units of e^//i)is a 
very robust integer, as long as the Fermi level remains in a spectral gap. However, the previous 
argument does not work for any kind of perturbation. It is valid provided the perturbation of the 
Hamiltonian is given by a matrix valued function 6H{k) satisfying the same covariance condition 



(|T^). Otherwise, the previous construction does not tell us anything about the homotopy 
invariance of the Hall conductance. This is exactly the limitation of the TKN2 result. More 
precisely, no disorder is allowed here and, in addition, only a rational magnetic field can be 
treated in this way, a quite unphysical constraint. For indeed, in most experiments the variation 
of the filling factor is obtained through changing continuously the magnetic field. Even though 
rational numbers are dense in the real line, the previous expression (|13|) depends in an explicit 
way on the denominator q of the fraction so that we have no way to check whether the integer 
obtained is robust when the magnetic field is changed. 

The drastic condition of having no disorder is also disastrous: the periodicity of the Hamil- 



tonian excludes localized states, and the argument in Section ^^ shows that no plateaux can 
be observed in this case. Therefore we have not yet reached the goal. 

Nevertheless, the main result of TKN2 is the recognition that the Hall conductance can 
be interpreted as a Chern character. This is a key fact in understanding why it is quantized 
and robust as well. Moreover, the topology underlying this Chern character is that of the 
quasimomentum space, namely the periodic Brillouin zone which is a 2-torus. From this point of 
view, the Laughlin argument in its original form may be misleading since it gives the impression 
of emphasizing the topology of the sample in real space, which has nothing to do with the IQHE. 

The solution to the TKN2 limitation is given by the Non-Commutative Geometry of the 
Brillouin zone where both of the above restrictions can be dropped. We will see that the Kubo- 
Chern relation still holds, that the Chern character is still a topological invariant and that it 
is an integer. But we will also discover something new, characteristic of Non-Commutative 
Geometry, namely that localization can be expressed in a very clear way in this context and 
that the conclusion will be valid as long as the Fermi level lies in a region of localized states. 



3 The Non- Commutative Geometry of the IQHE 

The four-traces way 

In this chapter the strategy and the main steps of the non- commutative approach are given. In 
particular we will introduce the four different traces that are technically needed to express the 
complete results of this theory. The first one is the usual trace on matrices or on trace-class 



operators. The second one, introduced in the Section |3.1| below, is the trace per unit volume 
which permits to compute the Hall conductance by the Kubo formula. The third one is the 
graded trace or supertrace introduced in Section |3]^. This is the first technical tool proposed 
by A. Connes IQ to define the cyclic cohomology and constitutes the first important step in 
proving quantization of the Hall conductance [1^. The last one is the Dixmier trace defined 
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by Dixmier in 1964 |^ and of which the importance for Quantum Differential Calculus was 



emphasized by A. Connes [^, |^, |§]. It will be defined in Section O but we already explain in 



Section |3.3| below how to use it in order to relate localization properties to the validity of the 
previous results. 

3.1 The non-commutative Brillouin zone 

In a D-dimensional perfect crystal, the description of observables is provided by Bloch's theory. 
Owing to the existence of a translation group symmetry, each observable of interest commutes 
with this group. It is a standard result in solid state physics that such an observable is a matrix 
valued periodic function of the quasi-momenta k. The matrix indices usually label both the 
energy bands and the ions in the unit cell of the crystal whenever it is not a Bravais lattice. 
The period group in the quasi-momentum space is the reciprocal lattice. 

If the crystal suffers disorder or if a magnetic field is turned on, this description fails. In 
most situations however, physicists have found ways to overcome this difficulty. For instance, 
impurities are treated as isolated objects interacting with the electron Fermi sea using per- 
turbation theory. Magnetic fields in 3-dimensional real crystals are usually so small that a 
quasi-classical approximation gives a very good account of the physical properties. 

Still the conceptual problem of dealing with the breaking of translation symmetry remains. 
It has not been considered seriously until new physical results forced physicists to face it. 
One important example in the past was Anderson localization due to a random potential 



28| , P9| , p3i |32| , |3J]. Even though important progress has been made, the solution is still in a 
rather rough shape and numerical results are still the only justification of many intuitive ideas 
in this field. When arriving at the QHE, the role of localization became so important that 
there was no way to avoid it. 

The main difficulty in such cases is that translates of the one-electron Hamiltonian no 
longer commute with the Hamiltonian itself. Nevertheless, the crystal under study is still 
macroscopically homogeneous in space so that its electronic properties are translation invariant. 
Thus choosing one among the translates of the Hamiltonian is completely arbitrary. In other 
words, the one-electron Hamiltonian is only known up to translation. Our first proposal to 
deal with this arbitrariness is to consider all the translates at once. More precisely, our basic 
object is the observable algebra A generated by the family of all translates of the one-electron 
Hamiltonian in the space G, where G is either R^ or Z^. 

It is remarkable that the algebra A exists in the periodic case as well and coincides then 



with the algebra of matrix- valued, periodic, continuous functions of the quasi-momenta |£T|. In 
this case, due to the periodicity with respect to the reciprocal lattice, it is sufficient to consider 
quasi-momenta in the first Brillouin zone B. Then B is topologically a D-dimensional torus. 
Therefore the matrix elements of a typical observable are continuous functions over B. 

Proceeding by analogy, we will consider the algebra A in the non-periodic case as the 
non- commutative analogue of the set of continuous functions over a virtual object called the 
"non-commutative Brillouin zone" p^ , |61|] . This is, however, ineffective as long as we have not 
defined rules for calculus. In the periodic case physical formulae require two kinds of calculus 
operations, namely integrating over the Brillouin zone and differentiating with respect to quasi- 
momenta. 

If A is an observable in a periodic crystal, its average over quasi-momenta is given by 
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<A>= f -T^Tr{A{k)) 



Here |B| represents the volume of the Brillouin zone in momentum space. It turns out that 
such an average coincides exactly with the "trace per unit volume" , namely we have 

<A>= lim ^Tta{Aa) = T{A) , (17) 

where A is a square centered at the origin, |A| is its volume, TrA and A^ are the restrictions of 
the trace and the observable A respectively to A. 

In a non-periodic, but homogeneous crystal it is still possible to define the trace per unit volume 
T{A) of an observable A. Formula (|l^ however becomes ambiguous in this case, because the 



limit need not exist. We will show in Section p.6| how to define T in general. It will satisfy the 
following properties: 

(i) T is linear, like the integral; 

(ii) T is positive, namely if A is a positive observable, T{A) > 0; 

(ill) T is a trace, namely, even though observables may not commute with each other, we still 
have r{AB) = r{BA). 

Differentiating with respect to quasi-momenta can be defined along the same line. In the 
periodic case, we remark that the derivative djA = dA/dkj is also given by 

djA = -t[Xj,A] , (18) 



where X = (Xi, . . . , Xn) is the position operator. We will use the formula ([T8| ) as a definition 
of the derivative in the non-periodic case. Owing to the properties of the the commutator it 
satisfies: 

(i) dj is linear; 

(ii) dj obeys the Leibniz rule dj{AB) = {djA)B + A{djB); 

(iii) dj commutes with the adjoint operation, namely dj{A*) = (djA)*. 

Such a linear map on the observable algebra is called a "*-derivation" . In our case it is moreover 
possible to exponentiate it, for, if 6* = (^i, . . . , 6d), we find 

e'-^{A)=e-''-'^Ae''-'^ , 

where 9- X = (6'iXi H h 9dXd) and 9 ■ V = {9idi H h 9DdD)- Equipped with the 

trace per unit volume and with these rules for differentiations, A becomes a "non- commutative 
manifold", namely the non- commutative Brillouin zone. 
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3.2 Hall conductance and non- commutative Chern character 

Using the dictionary created in the previous section between the periodic crystals and the non- 



periodic ones and in view of eq. (|T3[), one is led to the following formula for the Hall conductance 
at zero temperature 

an = ^Ch(Pp) , Ch(Pp) = ^T{PF[d,PF, d^Pp]) , (19) 

where Pp is the eigenprojection of the Hamiltonian on energies smaller than or equal to the 
Fermi level Ep. We will justify this formula in Section ( |4.2| ) below within the framework of the 
relaxation time approximation. However this formula is only valid if: 

(i) the electron gas is two-dimensional (so D = 2); 

(ii) the temperature is zero; 

(ill) the thermodynamic limit is reached; 

(iv) the electric field is vanishingly small; 

(v) the collision time is infinite; 

(vi) electron-electron interactions are ignored. 



We will discuss in Section (4^) the various corrections to that formula whenever these conditions 
are not strictly satisfied. 

The main result of this paper is that the non- commutative Chern character Ch(Pp) of the 
Fermi projection is an integer provided the Fermi level belongs to a gap of extended states (see 
Theorem Q). Moreover, we will show that it is a continuous function of the Fermi level as 
long as this latter lies in a region of localized states (see Theorem |l|). This last result implies 
that changing the filling factor creates plateaux corresponding to having the Fermi level in an 
interval of localized states (see Proposition [l^). In addition, whenever the Hall conductance 
jumps from one integer to another one, the localization length must diverge somewhere in 
between (see Corollary (^). 

In order to get this result we will follow a strategy introduced by A. Connes in this context 
^, but originally due to M. Atiyah |6^. Before describing it, we need some notations and 



definitions. We build a new Hilbert space which is made of two copies 7i+ and Ti.- of the 
physical Hilbert space Ti of quantum states describing the electron. This is like working with 
Pauli spinors. In this doubled Hilbert space 7i = TC+ © 7i_ , the grading operator G and the 
"Hilbert transform" F are defined as follows: 



6= r , . f= X i-^' . (20) 





where X = X1 + 1X2 (here the dimension is D = 2). It is clear that F is selfadjoint and satisfies 
F^ = 1. An operator T onTi will said to be of degree if it commutes with G and of degree 1 
if it anticommutes with G. Then, every operator on 7i can be written in a unique way as the 
sum T = Tq + Ti where Tj has degree j. The graded commutator (or supercommutator) of two 
operators and the graded differential dT are defined by 
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[T, T']s = TT' - ^^)de9{T)de9iT')r^,j, ^ ^j, ^ ^p^ y|^ _ 

Then, (PT = 0. The graded trace Tr5 (or supertrace) is defined by 



Tr5(T) = -Tr^(GF[F, T]s) = Tr^(r++ - uT__u) , (21) 



1, 
2 

where u = X/\X\ and T++ and T are the diagonal components of T with respect to the 

decomposition of Ti. It is a hnear map on the algebra of operators such that Trs'(TT') = 
Trs'(T'r). Moreover, operators of degree 1 have zero trace. However, this trace is not positive. 
Observables in A will become operators of degree 0, namely A & A will be represented by 
A = y4+ © A_ where A± acts as A on each of the components T-C± = Ti. 

With this formalism, A. Connes H] gave a formula which was extended to the present 



situation in [|lg] to the following one (see Section 



Ch(Pp) = / dP{uj)TTs{PFdpFdPF) , (22) 

Jn 

where the space Q represents the disorder configurations while the integral over dP is the 

average over the disorder. We will show that the integrand does not depend upon the disorder 

configuration. Moreover, we will show that the right hand side of (|22|) is P-almost surely equal 

to the index of a Fredholm operator, namely 

Ch(P^) = Ind(P^F+_ Tph^)) ■ (23) 

The Chern character is therefore an integer, at least whenever it is well-defined. 

One may wonder what is the physical meaning of this integer. Actually, an answer was 
provided quite recently by Avron et al. [0. Considering the definition of the graded trace in 
formula (|2T|) , they interpreted the index found for the Chern character as the relative index 
of the projections Pp and uPpu. This is to say even though these projections are infinite- 
dimensional, their difference has finite rank giving rise to a new index called the "relative 
index". Then they argue that u represents the action of the singular Laughlin gauge (see 
Section [27^ ) and that this difference can be interpreted as a charge transported to infinity as 
in the original Laughlin argument [0]. 

3.3 Localization and plateaux of the Hall conductance 



In the last Section ( ^.2] ), we have explained the topological aspect of the Hall conductance 
quantization. However, these formulae only hold if both sides are well defined. For instance, it 
is not clear whether the Fermi projection is differentiable. If it is not, what is the meaning of 



formula (19) ? This is the aspect that we want to discuss now. 

First of all, using the Schwarz inequality, a sufficient condition for formula (|T9|) to hold is 
that the Fermi projection satisfy T(|VPpp) < oo if we set V = (9i,92). It is important to 
remark that this condition is much weaker than demanding Pp to be differentiable. This is 
rather a non-commutative analogue of a Sobolev norm (or the square of it). By definition, this 
expression can be written as 

T{\VPp\^) = I dP{u) I d^x\ < 0\Pp\x > Hfp . (24) 

Jn Jx&G 



Note that equation (24) can be defined in any dimension D. From this expression, we see 



that the finiteness of the Sobolev norm is related to the finiteness of some localization length. 
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Indeed, we will show that only the part of the energy spectrum near the Fermi level gives a 
contribution to formula (|2^) so that Pp can be replaced by the projection Pa onto energies 
in some interval A close to Ep. On the other hand, the matrix element < 0|Pa|3; > gives a 
measure of how far states with energies in A are localized. Thus, the Sobolev norm is some 
measurement of the localization length. We will develop this fact in Section (|5.2| ) below. One 
way of defining a localization length consists in measuring how far a wave packet goes in time. 
This gives 

f{A) = hm U^ dt f dP{u) < 0||(X(t) -X)no > , 
T^oo T Jt=o Jn 

where the time evolution is governed by if Pa and where |Ap = J2i ^i^i if ^ = (^i, • • • , Ad). 
In this expression, X is the position operator and X{t) is the evolution of X under the one- 
electron Hamiltonian after time t. Let now A/" be the density of states (DOS) already introduced 



in Section |2.4| . It can be defined equivalently as the positive measure Af on the real line such 



that for any interval A 

T(Pa) = / dAfiE) . 
Ja 



That these two definitions coincide is guaranteed by a theorem of Shubin |61, p3]. We will show 



that if /^(A) is finite, there is a positive A/'-integrable function 1{E) over A such that 

l\A') = [ dU{E)l{Ef , 

for any Borel subset A' of A (see Theorem pTSl) . The number 1{E) will be called the "localization 
length" at energy E. Moreover, we also find 

r(|VPAf)< / dM{E)l{E)\ 

JA' 

such that the mapping Ep E A \—>- Pp is continuous with respect to the Sobolev norm (Theorem 



The previous argument shows that, whenever the localization length at the Fermi level is 
finite, the Chern character of the Fermi projection is well defined. Now let us consider the right 
hand side of formula (^21) . In a recent work A. Connes [^0|, |^ proposed to use another trace 



which was introduced in the sixties by J. Dixmier |B9. This Dixmier trace Trnj^ will be defined 



in Section ( [4.4| ) below. Let us simply say that it is a positive trace on the set of compact 
operators, such that any trace-class operator is annihilated by it whereas if T is a positive 
operator with finite Dixmier trace, then T^^*^ is trace class for any e > 0. In [Q, A. Connes 



proved a formula relating the so-called Wodzicky residue of a pseudodiferential operator to its 
Dixmier trace. Adapting this formula to our situation leads to the following remarkable result 
(see Theorem ^: 

T{\VPp\') = -TiU\dpF\') . 

Note that this formula only holds for electrons in two dimensions. In higher dimensions, things 
are more involved. Therefore, we see that as soon as the localization length is finite in a 
neighborhood of the Fermi level, the operator dPp is square summable with respect to the 
Dixmier trace implying that its third power is trace-class. In view of the formute (|2T|j2^ ), 
we see that this is exactly the condition for the existence of the right-hand side of eq. (|2^) . 
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Moreover, the continuity of Pp in Ep with respect to the Sobolev norm iniphes that the index 
in formula (|23| ) is constant as long as the Fermi level stays in a region in which the localization 
length is finite. Therefore the Chern character is an integer and the Hall conductance is 
quantized and exhibits plateaux. 

In this way, the mathematical frame developed here gives rise to a complete mathematical 
description of the IQHE, within the approximations that have been described previously. 



3.4 Summary of the main results 

Let us summarize our mathematical results in this section. 

Theorem 1 Let H = H* be a Hamiltonian affiliated to the C*-algehra A = C*{^ x G,B) 
defined in Section ^^. Let P be a G-invariant and ergodic probability measure on fl. Then we 
have the following results: 

i) (Kubo-Chern formula) In the limit where 

a) the volume of the sample is infinite, 

b) the relaxation time is infinite and 
a) the temperature is zero, 

the Hall conductance of an electron gas without interaction and described by the one- 
particle Hamiltonian H is given by the formula 

2 2 

aH = jChiPp) = j2m riPF[diPF,d2PF]) , 

where Pp is the eigenprojection of H on energies smaller than or equal to the Fermi level 
Ep and T is the trace per unit volume associated to P. 

ii) (Quantization of the Hall conductance) If in addition Pp belongs to the Sobolev space 
S associated to A, then Ch(Pp) is an integer which represents the charge transported at 
infinity by a Laughlin adiabatic switching of a flux quantum. 

ill) (Localization regime) Under the same conditions as the ones in ii), the direct conduc- 
tivity vanishes. 

iv) (Existence of Plateaux) Moreover, if A is an interval on which the localization length 
/^(A) defined in Section \373{ is finite, then as long as the Fermi energy stays in A and is 
a continuity point of the density of states of H, Ch(Pp) is constant and Pp belongs to 
the Sobolev space S. 



The following corollary is an immediate consequence. It was proved in pT] , |13 



Corollary 1 // the Hall conductance cth jumps from one integer to another in between the 
values Vi and v^ of the filling factor, there is an energy between the corresponding values of the 
Fermi levels at which the localization length diverges. 

Strictly speaking, this theorem has been completely proved only for the case of a discrete 
lattice (tight-binding representation). Most of it is valid for the continuum, but parts of the 
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proofs require extra technical tools so that the proof of this theorem is not complete in this 
paper. We postpone the complete proof of it for the continuum case to a future work |7^ . 



As a side result, we emphasis that we have developed a non-commutative framework valid 
to justify the transport theory for aperiodic media (see Section ^]T] below). It allows us to give 
rigorous estimates on the error terms whenever the conditions of the previous theorem are not 
strictly satisfied (see Section |4.3| below). We will not give the mathematical proofs that these 
errors are rigorous bounds here even though they actually are. This will also be the main topic 
of a future work. But we have estimated them and we show that they are compatible with the 
accuracy observed in the experiments. 

Another result which is actually completely developed here due to its importance in the 
Quantum Hall Effect, concerns the definition and the properties of the localization length. We 
give a non- commutative expression for it and show how it is related to the existence of plateaux. 



The main results in this direction are the Theorems |T3| and in Section p.2| . We also show 
that the localization length is indeed finite in models like the Anderson model for disordered 
systems for which proofs of exponential localizations are available. 

The proof of Theorem |l] will be divided up into a number of partial results; there will be 
no explicit paragraph 'Proof of Theorem ^', let us therefore outline the main steps. In 
Sections ^^ and pl6| , we give a precise mathematical description of a homogeneous Schrodinger 
operator and its hull and we construct the observable algebra. In Chapter ^ the Kubo formula 
for the Hall conductance is derived and we present the (non-commutative) geometrical argument 
for the integer quantization of the Hall conductance (compare Theorem |Tl|). That the proven 
index theorem is precisely valid under the weak localization condition Pp & S results from 
Theorem p. Point iv) is proved in Chapter p. 

3.5 Homogeneous Schrodinger's operators 



Most of the results of the next two sections have already been proved in ||6T|, Section 2]. We 
will only give the main steps that are necessary in this paper for the purpose of the IQHE. 
In earlier works [111, BTI, one of us has introduced the notion of a homogeneous Hamilto- 



nian. The purpose of this notion is to describe materials which are translation invariant at a 
macroscopic scale but not necessarily at a microscopic one. In particular, it is well suited for 
the description of aperiodic materials. In such a medium, there is no natural choice of an origin 
in space. If if is a Hamiltonian describing one particle in this medium, we can replace it by 
any of its translates Ha = U{a)HU{a)~^, a G R^; the physics will be the same. This choice 
is entirely arbitrary, so that the smallest possible set of observables must contain at least the 
full family {Ha, a G R^}; this family will be completed with respect to a suitable topology. 
We remark that H need not be a bounded operator, so that calculations are made easier if we 
consider its resolvent instead. Let us define precisely what we mean by "homogeneity" of the 
medium described by H. 

Definition 1 Let TC be a separable Hilbert space. Let G be a locally compact group (for instance 
'RP or TP ). Let U be a unitary projective representation of G, namely for each a E G, there 
is a unitary operator U{a) acting on H such that the family U = {f/(a);a G G} satisfies the 
following properties: 

(i) U{a)U{b) = U{a + b)e^'^^"''^^ Va, b E G, where (f){a, b) is some phase factor. 

(a) For each ip ETi, the map a E G ^-^ U{a)il) E 7i is continuous. 



23 



Then a self adjoint operator H on Ti is homogeneous with respect to G if the family {Ra{z) = 
U{a){zl — H)~^U{a)~^; a G G} admits a compact strong closure. 

Remark. A sequence (A„)„>o of bounded operators on H converges strongly to the bounded 
operator A if for every tp & 7i, the sequence {Anip} of vectors in Ti converges in norm to Aijj. 
The set considered in the definition has therefore a strong compact closure if, for given e > 
and for a finite set ipi, . . . , ipN of vectors in Ti, there is a finite set ai, . . . , a^ in G such that for 
every a & G and every I < j < N, there is 1 < i < m such that || {Ra{z) — RaX^))'^] 11^ ^• 
In other words, the full family of translates of R{z) is well approximated on vectors by a finite 
number of these translates; this finite number then repeats itself infinitely many times up to 
infinity. 

The virtue of this definition comes from the construction of the "hull" . Let z belong to the 
resolvent set p{H) of H and let H be homogeneous. We denote by Q{z) the strong closure of 
the family {Ra{z) = U{a){zl — H)^^U{a)~^\ a G G}; it is therefore a compact space. Moreover, 
it is endowed with an action of the group G by means of the representation U . This action 
defines a group of homeomorphisms of VL{z). Thanks to the resolvent equation, it is quite easy 
to prove that if z' is another point in the resolvent set p{H), the spaces Q{z) and Q{z') are 
homeomorphic [^. Identifying them gives rise to an abstract compact space Q endowed with 



an action of G by a group of homeomorphisms. li lu E il and a G G, we will denote by T^cj the 
result of the action of a on uj, and by Ruj{z) the representative of uj in VL{z). Then we have 

U{a)R^{z)U{ay = Rt^^z) , 

RM') - RM) = (^ - z')R^{z')RM) ■ (25) 

In addition, z ^-^ Ru){z) is norm-holomorphic on p{H) for every cj G fi, and the application 
UJ \-^ Ru){z) is strongly continuous. 

Definition 2 Let H he an operator, homogeneous with respect to the representation U of the 
locally compact group G on the Hilbert space Ti. Then the hull of H is the dynamical system 
{Q,G,T) where fl is the compact space given by the strong closure of the family {Ra{z) = 
U{a){zl — H)^^U{a)^^] a G G}, and G acts on Q through T. 



In general, the equation ( pS]) is not sufficient to ensure that R^{z) is the resolvent of some 
self-adjoint operator if^^, for indeed, one may have Ruj{z) = if no additional assumption is 
demanded. A sufficient condition is that H be given by Hq + V where: (i) Hq is self-adjoint 
and G-invariant, (ii) V is relatively bounded with respect to Hq, i.e. || {z — Hq)^^V \\< oo, 
(iii) lim|,|^„o II (z - Ho)-^V 11= 0. Then, R^{z) = {1 - (zl - Ho)-^V^}-\zl - Hq)-^ where 
{zl — Hq)^^V^ is defined by the strong limit of {zl — HQ)^^Va^, which obviously exists. So 
Rui{z) is the resolvent of Hq + K,. 

If if is a Schrodinger operator, the situation becomes simpler. Let us consider the case of a 
particle in R^ with mass m and charge g, submitted to a bounded potential V and a uniform 
magnetic field B with vector potential A. We will describe the vector potential in the symmetric 
gauge, namely A = {— 8x2/2, Bxi/2). The Schrodinger operator is given by 

H=7^ Y.iPj-QA,r + V = HL + V . (26) 
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The unperturbed part Hi is translation invariant, provided one uses magnetic translations p7 
defined by (if a e R^, ^ G L'^{I{?) ) 

U{a)%p{x) = exp < X A a> iIj{x — a) , (27) 

where a; A a = Xia2 — X2ai. It is easy to check that the operators U{a) form a projective unitary 
representation of the translation group R^. The main result in this case is given by 



Theorem 2 Let H he given by eq. ( ^2d^ ) with V a measurable essentially bounded function on 
R2. Then 



(i) H is homogeneous with respect to the magnetic translations (]F/[j; 



(a) the hull uj of H is homeomorphic to the hull ofV, namely the weak closure of the family 
{1/(. -a);aeR2} m L^(R2); 

(Hi) there is a Borelian function v on Q such that, if we denote by V^ the bounded function 
representing the point u & Q, then Vi^{x) = v{T~^uj) for almost every x G R^ and all u G VL. 
If in addition V is uniformly continuous and bounded, then v is continuous. 



The proof of this theorem can be found in |61|, Section 2.4]. 



In many cases, it is actually better to work in the tight-binding representation. The reason 
is that only electrons with energies near the Fermi level contribute to the current. One usually 
defines an effective Hamiltonian by reducing the Schrodinger operator to the interval of energies 
of interest |]n[. The Hamiltonian can then be described as a matrix H{x,x') indexed by sites 



in the lattice Z^ acting on elements of the Hilbert space £^(Z^) as follows 

HiIj{x) = y H{x, x')exp l tn—x A x \ ipix') , V e f{Z'^) , (28) 

where is the flux in the unit cell, whereas 0o = h/q is the flux quantum. In most examples, 
one can find a sequence / such that I]aez2 f{3) < oo and \H{x,x')\ < f{x — x'). In particular, 
H is bounded and there is no longer a need to consider its resolvent. Let now U be the unitary 
projective representation of the translation group Z^ given by 

U{a)ilj{x)=exp\nT^xAa\tlj{x-a), ^p e f{Z'^) . (29) 



Theorem 3 Let H be given by eq. (^23(). Then H is homogeneous with respect to the projective 
representation U (eq. ^2^)) of the translation group. Moreover, iffl is the hull of H, there is a 
continuous function h on Q x Z"^ , vanishing at infinity, such that H^{x, x') = h{T~^uj, x' — x), 
for every pair (x, x') G Z^ and uj eVL. 



Again the proof can be found in |^T], Section 2.4] 



3.6 Observables and calculus 

In the previous section we have constructed the hull of a Hamiltonian. It is a compact space that 
represents the degree of aperiodicity of the crystalline forces acting on the charge carriers. For 
disordered systems, the hull is nothing but the space of disorder configurations. In principle, the 
algebra of observables should be constructed from the Hamiltonian H by taking all functions of 
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H and its translates. However, we proceed in a different way giving a more explicit construction 
which may give a bigger algebra in general, but will be easier to use. 

Let f2 be a compact topological space equipped with a R^-action by a group {T"; a G R^} 
of homeomorphisms. Given a uniform magnetic field B, we can associate to this dynamical 
system a C*-algebra C*(n x R^,S) as follows. We first consider the topological vector space 
Ck{^ X R^) of continuous functions with compact support on fi x R^. It is endowed with the 
following structure of a *-algebra by 



AB{uj,x) 



(iqB 
< f. 



2 d yA{uj, y)B{T ^w, x - y)exp < -^x A y 



(30) 



A*{uj,x) = A{T-^uj,-x) 



where A, B & C^{Q x R^), u & Q and x G R^. For u & Q, this *-algebra is represented on 
L2(R2) by 



7r^{A)ij{x) 



r iqJS 
^ £yA{T-''uj,y- x)exp | ^^ A x \ %l){y) 



2/'-d2\ 



V^ G L^(R- 



(31) 



namely, vr^ is linear, 7r^(yli?) = -n^^{A)-nt^{B) and HuiiA)* = ni^{A*). In addition, 7r^(A) is a 
bounded operator for || 7r^(A) ||<|| / ||oo,i where 



A lloo 1= max < sup / (fy\A{uj,y)\, sup (fy\A*{uj 
I ujen "'R.2 wen JR^ 



This defines a norm which satisfies || AB ||oo,i<|| A ||oo,i|| -B lloo.i, || A ||oo,i=|| A* ||oo,i- Since 
A has compact support, its (cxd, l)-norm is finite. We will denote by L°°'^[Q x R^; i3) the 
completion of Ci^{Q x R^) under this norm. We then remark that these representations are 
related by the covariance condition 



U{a)7r^{A)U{a)-' = n^^a^iA) . 



Now we set 



II A 11= sup II TT^{A) II , 

which defines a C*-norm on L°°'^{Q. x B?;B). This permits us to define C*{Q x R^, ;B) as the 
completion of Cfi{Q x R^) or of L°°'^{Q x R^; B) under this norm. Clearly, the representations 
TT^ can be continuously extended to this C*-algebra. We remark that this algebra has no unit. 



In particular, for the Hamiltonian (|26|) we get ||6l|. Section 2.5]: 



Theorem 4 Let H be given by eq. ^d^ ) where A is the vector potential of a uniform magnetic 
field in the symmetric gauge and let V be in L°°(R^). We denote by Vt its hull. Then for each 
z in the resolvent set of H and for every x G R^, there is an element R{z; x) E C*{fi x R^, B), 
such that for each co E Q, n^^R^z; x)) = {zl — Hrp-s^)~^. 

In the discrete case (see eq. (p8| j3|), the same construction holds provided we replace R^ by Z^ 
and the integral over R^ by a sum over Z^, namely 
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Theorem 5 Let H be given by eq. (^^ where is the magnetic flux through the unit cell and 
00 = h/q is the flux quantum. We denote by Q its hull as in Theorem |^. Then the function h 
appearing in Theorem ^ belongs to C*{Q x Z^,B) and for each u & Q, TTi^{h) = H^. 

We remark that in the discrete case, the function l{u, x) = (5-g is a unit of C*{Q x Z^, B); this 
is the main difference between the continuous and the discrete case. In the non-commutative 
terminology, the discrete case corresponds to a compact non-commutative manifold, whereas 
the continuous case corresponds to a locally compact, but not compact, non- commutative 
manifold. In the sequel, many results will hold for both of the constructed C*-algebras. We 
therefore introduce the notation A for both of them. The algebras of functions with compact 
support which are at the base of the construction are denoted by ^o- Moreover, let G be the 
physical space; it is either R^ or Z^. 

If the Hamiltonian H is unbounded, it does not belong to the C*-algebra A. However, we 
have seen in Theorem ^ that its resolvent belongs to A. As a consequence, f{H) belongs to A 
for every continuous function / on the real line vanishing at infinity. This leads to the following 
definition. 

Definition 3 Let A be the C*-algebra C*{Q x G,B) defined above and let Ti be the Hilbert 
space L'^{G) endowed with the projective representation U defined in eq.(^^, ^ ). We will say 



that a selfadjoint homogeneous operator H is affiliated to A whenever its hull can be embedded 
in the dynamical system {Q, G, T) and if for each z in its resolvent set, there is an element 
R{z) e A such that 

{zl-H^)-' = 7i^{R{z)), 
for all uj E Q. 

By 'embedded' we mean that the hull is a T-invariant subset of Q. 

The two rules of calculus, namely integration and derivations, are now constructed as follows. 
First of all, let P be a G-invariant probability measure on Q. For most of the results of this 
paper, the choice of such a measure is irrelevant. We will discuss its physical relevance in 
Section ^.7| and Section |S.2| . A trace on A is constructed as follows. If A G ^o we set: 

Tp{A) = f dP{uj)A{uj,6) . 
Jn 

This defines a positive trace (see Section |3.1|). It is faithful (namely Tp{AA*) = implies 



A = 0) whenever the support of P is Q. The trace is normalized if G = Z^ (namely 7p(l) = 1) 
and is unbounded if G = R^. It actually coincides with the trace per unit area. More precisely 



we have |T3, |61| : 



Proposition 1 Let A belong to Aq. Then for P -almost all uj 's 

rp(A) = hm^TrA(vr^(A)), (32) 

where A denotes a sequence of squares in G centered at the origin and covering G and Tt\ is 
the restriction to A of the usual trace. 



27 



In the sequel, we will drop the subscript P. 

Given a selfadjoint element H & A, we define its DOS as the positive measure Af on the 
real line such that for any continuous function / with compact support on R 



TifiH)) = f dM{E)f\E) 
Jr. 



In view of ([3^), this definition agrees with the definition given in eq. (|^). 
If p > 1, we denote by Lp{A, T) the completion of ^o under the norm 

\\A\\L.= (T{{AA*Y'^f'' 



In particular, one can check that the space L'^{A,T) is a Hilbert space (GNS construction) 
identical to L'^iVt xG.dP® dg). The map G L'^{A, T) ^ Acj) e L^{A, T) for A e ^ defines a 
representation tigns oi A. The weak closure L°°(^, T) = ncNsi.-^)" is a von Neumann algebra. 
By construction, the trace T extends to a trace on this algebra. We remark that if if is a self- 
adjoint element of A, its eigenprojections are in general elements of the von Neumann algebra 

L^iA,r). 

Let us give another characterization of the von Neumann algebra L°° {A, T) which can be 



found in ^Sf. Let W be the set of weakly measurable families tu G fi i-^ A^ of bounded 
operators on L?'{G) which are covariant 

U{a)A^U{ay^ = At^^ , a eG, cu eVl . 

and P-essentially bounded. This latter means that the norm of A^ is bounded in u except 
possibly on a subset of zero probability with respect to P. We endow W with the norm 

II A ||p= P-esssup II A^ ||£(L2(G)) , AeW. 

Sum, product and adjoint of elements of W are defined pointwise in the obvious way. Then 
Connes |^5[| proved that W is a von Neumann algebra, namely a C*-algebra with predual ||6^] , 
and that 

Theorem 6 L°°{A,T) is canonically isomorphic to the von Neumann algebra W of P-essen- 
tially bounded measurable covariant families of operators on L'^{G). 

Actually, this isomorphism is obvious if we realize that the Hilbert space L^{A,T) of the 
GA^S'-representation of A with respect to the trace T can be written as the direct P-integral 
of L'^{G). We will not give details. A consequence of this result is that the family {tt^; uj ^ Q} 
of representations of A extends to a (faithful) family of (weakly measurable) representations of 
W. In particular, any spectral projection P of a Hamiltonian H & A can be seen as a covariant 
P-essentially bounded family P^ of projections, where P^j is an eigenprojection of H^. 
The differential structure is obtained through the derivations defined on ^o by 

djA{uj,x) = iXjA{u,x) . (33) 



It is a *-derivation in the sense given in Section |3.1| . By exponentiation it defines a two- 
parameter group of *-automorphisms given by 



p^{A){uj,x) = e'^"'A{uj,x) , 
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where A; G R^ and k ■ x = kiXi + k2X2. We notice that in the discrete case, k is defined modulo 
27rZ^. Introducing on L'^{G) the position operator X = (Xi,X2), namely the multiplication by 
X, we get 

nUPki^)) = e-'^-^r^MY''-'^ , ^^i.d,A) = -i[X^, ti^)] . (34) 

We remark that in the case of periodic media, this derivation is just differentiation in quasi- 
momentum space. We will denote by C^{A) the set of elements A G ^ for which the map 
k G R^ HH> p^{A) G ^ is iV-times continuously differentiable. If N is an integer, this is 
equivalent to say that || (9^92^4 ||< cxd for any pair a, b of integers such that a + b < N. In much 
the same way non-commutative Sobolev spaces can be defined. For the purpose of this work 
we will use the Sobolev space S = li?{A, T), namely the Hilbert space obtained by completion 
of ^0 under the Hilbert norm given by the inner product 

< A\B >^2= T{A*B) + TiWA* ■ VB) , A,B eAo , 

where V = ((9i,92) is the non-commutative gradient operator. 

We will finish this section by giving a technical result which will be used later on. 

Proposition 2 Let H be a self adjoint element in C^{A) where A = C* {Q x Z"^; B) and let M 
be its density of states. Then for any function f G L^{Il,dJ\f) we have: 

r{f{H)d,H) = 0. 

Sketch of the proof. By density, it is enough to prove this result for a smooth function / on 
R with compact support. Then one can write 



f{H) = f dt e-*^/(t) , 



where / is the Fourier transform of /. Since H is bounded, this integral converges in norm. 
Classical results on Fourier transform show that / is a smooth rapidly-decreasing function over 
R. Thus, it is sufficient to show the result for f{E) = e~**^. In this case, Duhamel's formula 
gives 



Jo 
Taking the trace, the left hand side vanishes, since the trace is p^-invariant. The right hand 
side is given by T(e~**^Vif). The proof is then complete. □ 

In the continuum case, a similar result holds, but it is technically more complicated. We will 
restrict ourselves to the case of Schrodinger operators. 



Proposition 3 Let H be given by eq. ( ^Bj ) with V G L°°(R^). We denote by Q the hull of 
V . Then V H is a well-defined self adjoint operator bounded from above by {H + bY'^ for some 
positive b. Moreover, if M denotes the DOS of H , for every M -measurable function f such that 
the map E i-^ Ef{E) is in -^^(R, dj\f), we have 

T{f(H)VH) = 0. 
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Sketch of the proof. We use the results of ||6T1 , Section 2.5] to conclude that e *^ admits an 



integral kernel F{t; uj, x) for 3?e(t) > because H is bounded from below. This kernel is jointly 
continuous in cj, x and holomorphic in t in this domain. Moreover, it decays rapidly at infinity 
in X (uniformly on compact subsets for the other variables). Then we use the same argument 
as in Proposition provided we replace e~**^ by e-i^+^^)^ for any e > 0. □ 



4 The Kubo-Chern formula 

This chapter is devoted to the Kubo formula and its relation to the Chern character. The first 
three sections are devoted to transport theory and are not treated on a completely rigorous 
footing. This is because technical proofs would require too much details and not shed more 
light on the physics. 



4.1 The relaxation time approximation 

The theory of transport is an essential tool in dealing with electronic properties of solids. There 
are numerous theoretical approaches with complexity varying from the Drude-Sommerfeld the- 
ory to the N-body framework. Whatever the starting point, the Greenwood-Kubo formulae 
for the transport coefficients, such as the electric or thermal conductivity, are the main conse- 
quences. They are widely used and accepted because of their correspondence with experimental 
results. Still their derivation from first principles is questionable. One does not really know the 
precise domain of validity of the linear response approximation. The occurrence of many new 
devices in solid state physics liable to test these ideas, such as mesoscopic systems, has raised 
this question again. 

It is not the purpose of this work to investigate that problem. However, we have seen in 
Section 12. 11 that in the classical Hall effect, the relation between the Hall current and the Hall 



voltage is linear although there is no dissipation mechanism. For this reason, one might expect 
that the derivation of Kubo's formula for the Hall conductance in a quantum system should 
not require a many particle theory. This is actually not true. We will see in this section that 
such a point of view is very naive and gives rise to paradoxa. Without dissipation, a quantum 
theory leads either to a vanishing conductivity or to an inffnite one in most physically sound 
situations. 

Moreover the great accuracy of the IQHE has been used in metrology for the deffnition of 
a new standard of resistance [^. It is thus necessary to derive a formula allowing the control 
of deviations from the ideal QHE. We have already indicated in Section |3.2| what the physical 



conditions are under which the IQHE is exact, namely under which it can be stated as a theorem 
within a well-deffned mathematical framework. 

This is why we propose a one-particle model including collision effects, such as interaction 
with phonons or other electrons, realizing thereby the so-called relaxation time approximation 
(RTA). We derive a Kubo formula for the conductivity tensor which allows the linear response 
approach to be justified and gives control on the order of magnitude of the deviation from the 
ideal theory. The construction of this model is based upon the phenomenological approach 
which can be found in standard books in solid state physics such as [1T8| , Chap. 1,2, 12, 13]. 

We consider the electron fluid of our system as a gas of independent fermions. Neglecting 
interactions between electrons is actually very rough and in many cases completely wrong. 
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Nevertheless the Landau theory of Fermi hquids shows that such an approximation is quite 
acceptable if the particles are actually quasi-particles "dressed" by the interactions |6^. In 



particular, their mass need not necessarily be the electron mass. Moreover their lifetime is 
finite, but it goes to infinity as their energy gets closer to the Fermi level. Quasi-particles 
carrying current are therefore stable at zero temperature. In the sequel, we continue to talk of 
electrons having in mind that we are actually treating quasi-particles. 

The advantage of this independent-electron approximation is that one can avoid using second 
quantization and can restrict oneself to a one-particle description. The constraint given by 
Pauli's principle is then represented by the use of the Fermi-Dirac distribution function. 

For this reason, our starting point will be the one-particle Hamiltonian H of the form already 



described in Section p.5| . It includes kinetic energy of the electrons as well as whatever time- 
independent forces there are acting on them in the crystal. These latter include the periodic 
potential created by the ions and the aperiodic corrections due to impurities, defects and other 
kinds of deformations. This Hamiltonian goes beyond the band theory since it may include 
Anderson localization for instance. We will assume that H is affiliated to the observable algebra 
A = C*{Q, G, B) where G = R^ or Z^ and Q is some compact space (see Section WM - In the 
grand canonical ensemble the thermal equilibrium at inverse temperature jS = l/ksT (where 
ks is the Boltzmann constant and T is the temperature) and chemical potential /i is described 
by means of the Fermi-Dirac distribution. In the algebraic language it means that, if A G ^ is 
a one-particle extensive observable, its thermal average per volume is given by 

< A >p,,= T{Afp,,{H)) , with fp,,{H) = (1 + e^(^-^))-i . (35) 

One important example for an observable is the current operator 

Using the differential structure (|55|), this can be written as 



n 
It is physically obvious that the average current vanishes since the velocity distribution is usually 
symmetric under the change of its sign. This is actually what happens within our framework 
because of the Propositions ^ and ^ and the formula (|35| ) above. Producing a non-vanishing 
average current requires imposing an external force such as an electromagnetic field. Let us 
consider the simplest case in which we superimpose a constant uniform electric field S on our 
system for time t > 0. Then the evolution of an observable is no longer given by H^^ but rather 
by H^^ = Hi^ — qS ■ X. We notice that while H is affiliated to A, this is certainly not true 
for this Hamiltonian since the position operator is not homogeneous. However, the evolution 
under this operator leaves A invariant. For indeed, whenever A is smooth in A, Heisenberg's 
equation reads 

f^4|ff„.^J + |.(VA)„. (36) 

If in particular H belongs to A, the right hand side of ( ^6|) stays in A so that we expect this 
equation to be integrable within A. This can be proved by use of the Trotter product formula 
)9|. We denote by r/j the evolution given by H, namely 
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,tH . ,tH 



rn{A) = e'-Ae"- . (37) 

Since H is affiliated to A, this evolution leaves A invariant and defines a one-parameter group 



of *-autoniorphisnis of A. Because of formula (|3^ , we also have 



/ \ qt£'X qte-X 

T^u> (Pf f (^)j = e "^T^u>{A)e'^r- . 
This last evolution also leaves A invariant. By means of the Trotter product formula we find 

e^^^^.%^(A)e-*^^-^^ = .- hm tt J(m p,,,^)^(A)) . 

Here, s-lim is the limit in the strong operator topology. So this defines a one-parameter group of 
automorphisms rj^ of the von Neumann algebra W. It represents the evolution of the observables 
after the electric field has been turned on. We will not investigate here whether this group of 
automorphisms leaves A itself invariant. Let us only notice that, whenever H is bounded, a 
Dyson expansion shows that r|^ is an automorphism of A. Moreover, we only have to work in 
the Hilbert space L'^{A, T) so that the previous result will be sufficient if H is unbounded. The 
new current at time t > is then formally given by 

At) = vfiJ) ■ 

Since the Hamiltonian H^ ^no longer commutes with if^^, the thermal average j(t) =< J{t) >^^^ 
of the current will no longer vanish in general. The macroscopic response we expect from a 
constant uniform electric field is a constant and uniform current, but the microscopic quantal 
forces lead to fluctuations of j(t) in time. To extract the response, it is therefore necessary to 
consider the time average of this current which is actually what one measures in experiment. 
For indeed, the typical relaxation time of an apparatus measuring the current is of the order of 
Ims to lyus (unless very short times are needed), and this has to be compared with the typical 
collision time of the order of 10~^'^s for the best conductors. Thus, we set 

j' , .-= lim - / dtj(t) = \im5 dte~^^j(t) . (38) 

However, our model does not take collisions into account. Indeed we have: 

Proposition 4 // the Hamiltonian H is bounded in A, the projection of the time and thermal 
averaged current ja g along the direction of the electric field £ vanishes. 

Proof. Let us compute £ ■ J it). Using the Heisenberg equation it is easy to see that 

£./(t) = ^, H{t)=r^^{H). 

Taking the time average gives us 

1 r ..c Tu^ H{T)-H 



-j^ dt£.J{t) 



T 

Since H is bounded in norm and since || H{t) \\ = \\ H ||, the right hand side vanishes as T ^ oo. 
D 

This result is easy to understand in the one-band approximation. In this case, the electric field 
produces a time shift of the quasi-momentum, namely A; G B i-^ k{t) = k — q£t/h. Taking 
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the time average is therefore equivalent to averaging over quasi-momenta, and this is exactly 
zero. We point out that this result applies also to models with disorder of the form given in 
(P5|). The Anderson model is a special example of such Hamiltonians. More generally, if we 
accept that transitions to bands far from the Fermi level are essentially forbidden, the effective 
Hamiltonian is always bounded and the Proposition ^ then leads to the vanishing of the current 
component parallel to £. 

For these reasons, the presented approach is definitely too naive. As we already pointed 
out, collisions occur after time periods very short compared to the time over which we average 
the current. These collisions produce a loss of memory in the current evolution and are actually 
responsible for the non- vanishing of the time average. Everything looks like as if time evolution 
has to be considered over short intervals only. We propose to take collision effects into account 
by means of the following model: the time evolution is described by the new time-dependent 
Hamiltonian Hcoiiit) = H — qS ■ X + Wcoiiit) where 

neZ 

In this expression, the t„'s represent collision times. They are labeled such that . . . , t_i < to = 
< ti < . . . < t„ < tn+i < . . .. Because these times occur randomly, we will assume that the 
time delays r„ = t„ — t„_i are independent identically distributed random variables distributed 
according to Poisson's law with mean value E(r„) = r 

The WnS are the collision operators. Their main effect is to produce a loss of memory 
during the time evolution of the current. This process should enforce thermal equilibrium. 
In particular, it should not modify the Fermi-Dirac distribution. The only way to respect 
this constraint is to force the W^s to commute with the Hamiltonian H, but to be random 
otherwise. More precisely, we will assume that the WnS are independent identically distributed 
random operators, commuting with the Hamiltonian H. Their distribution is supposed to be 
symmetric under the change of sign Wn ^^ —Wn- Let then k be the operator acting on A as 

k{A) = E(e r^ Ae « ) . 

It can be extended to a bounded operator on L^{A,T). Moreover, because of the change-of 
sign symmetry, it is selfadjoint. Let then Ah be the closed subspace of L'^{A, T) generated by 
bounded functions of H . We will then assume that there is < /« < 1 such that 

II k{A) \\<n\\A\\ , \/ Ae {Ah)^ 

Since Wn commutes with H it follows that k leaves Ah invariant so that it forces any operator 
A & A along the direction of H. We will not give explicit examples of such random operators 
because we will only use the collision efficiency operator k later on, so that it solely characterizes 
our model. The parameter k is purely phenomenological and represents an efficiency coefficient 
of the collision process. The smaller k, the more efficient are the collisions. We see that 
the relaxation time, namely the time after which there is a complete loss of memory of the 
initial evolution, has to be renormalized by the efficiency coefficient. We will see below that 
Trei = r/(l — k) is a good measure for this relaxation time. We will discuss later on how to 
choose this parameter to fit with real systems. 
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4.2 Kubo's formula 

We now follow the strategy defined in the previous section and compute the current average 
with the collisions taken into account. This requires to calculate the evolution operator S^{t) 
between time and t > 0, where ^ represents the random variables ^ = (r^, Wn)n>Q- It is well 
known that a kick term like WnS{t — tn) produces a contribution e*^"/^ in the evolution between 



times t„ — and t„ + 0, namely at the kick time ||70|, 0. Therefore, if n > 1 and t„„i < t < t^ 
we find: 



n-l 



i(t-tr,-i)Hg/h TT tWj/h i{tj-tj.i)H^/h 



S^{t) = eH*-*n-i)%/^ Y[ 



e ^' e 



We will now set Ch{A) = {i/h)[H,A] for A in a suitable dense subalgebra of A. This is a *- 
derivation of ^ because of (|37|). Therefore, it defines an anti-selfadjoint operator on L'^{A,T). 
Moreover, the evolution of observables is given by the automorphism 

n-l 
J=l 

where Cwji^) = (^/^)[^j) ^]- The operator 7]^^ may also be seen as a unitary on L'^{A,T). In 
view of eq. (|38|), the time average of this evolution requires the calculation of 



where E denotes the average over ^ and 6 > 0. After a straightforward calculation we find: 

This expression is meaningful because the operator in parentheses has a real part bounded from 
below by 6. The average current is then given by: 

J,,,,e<^) = J-'^{faH)vk^H)) . (39) 

We then remark that, for S = 0, the right hand side of (^) vanishes because k{f{H)) = 
f{H) and CnifiH)) = for any function /, whereas Ch is anti-selfadjoint. Subtracting this 
expression with S = and using the inner product < A\B >= T{A*B) in L'^{A, T), gives 

i./(^) = (|)' E ^^ < 9J,AH)\ _, \ ^H > . (40) 

If we assume that the non-linear term in the electric field is negligible, we can let 6 converge to 
zero. For indeed, VH and V fp^H) are orthogonal to Ah- Thus it is enough to consider the 
restriction of the operator 5+ ^—^ — Ch to the subspace orthogonal to Ah- This restriction has 
a real part bounded from below by (1 — K)/r > 0. Thus it is an invertible operator. If we do 
not neglect the non-linear term, we have to investigate more thoroughly what happens as 5 | 0. 
We will not discuss that matter here and postpone it to a future work. Letting 5 converge to 
zero, we will get the Kubo formula for the conductivity tensor: 
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Kubo's formula 

W-=^^^ with <y^, = {j^)<d,fpAH)\ ^^_^_^^ d.H> . (41) 

Let us reniark that the quantity in the bracket is dimensionless in 2D: the trace T is a trace 
per unit volume (so that it has the dimension of the inverse of an area) while the derivative dj 
has the dimension of a length. In addition hCn has the dimension of an energy. Therefore, we 
get a dimensionless in units of q^/h. 

If, in addition, the system is symmetric under rotation by 7r/2 in space, namely if the 
Hamiltonian H and the efficiency operator k are invariant under such a rotation, then the 
conductivity tensor can be written in matrix form as 

where the Hall conductance is the off-diagonal term an, while a// is the direct conductivity. 

Let us now consider the limit for which the IQHE is valid. Namely the electric field is 
vanishingly small, the temperature is zero and the relaxation time is infinite. In this limit, the 
Fermi distribution becomes 

lim^,^(/7) = Pp, (42) 

/3Too 

where Pp, the Fermi projection, is the spectral projection onto energy levels lower than the 
Fermi energy. Here the limit is taken with respect to the norm in L'^{A,T). Actually, this 
result is correct only if the Fermi level Ep is a continuity point of the DOS of H. Otherwise the 
eigenprojection P{Ef} corresponding to the eigenvalue Ep satisfies T{P^Ef}) > and therefore 
defines a non-zero element of L'^{A,T). This is what happens for instance in the case of the 
Landau Hamiltonian whenever Ep = {n + l/2)hujc. Moreover, we need the Fermi projection 
to be Sobolev differentiable, namely it has to to satisfy VPp G L'^{A,T), otherwise the formal 
limit is meaningless. We will see in Section |5.2| below that such a condition is related to the 
finiteness of localization length at the Fermi level. We also need to show that the limit of the 
derivative of the Fermi distribution is the derivative of its limit. Even though we know this to 
be true for physically reasonable conditions on H, we will not give the proof here but postpone 
it to a future work. 

On the other hand, in the limit of infinite relaxation time, we are left with the formal 
expression Cj^diH which is meaningless in general. If we formally diagonalize the Hamiltonian 
H, matrix elements of that expression are given by 

In particular it diverges whenever E ^ E' unless the numerator vanishes for some reason. This 
divergency does, however, not occur in the expression of the conductivity tensor. The reason 
is that the derivative VPp has non-vanishing matrix elements only between energies E and E' 
such that E < Ep < E' or E' < Ep < E. For indeed, if P is a projection and d is a derivation, 
we see that 

dP = d{P^), =^ dP={l-P)dPP + PdP{l-P). 
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So we need to consider the operators PpC^diH {\ — Pp) and (1 — PF)Cf/diH Pp only. We 
then obtain 

Lemma 1 // the Fermi level is not a discontinuity point of the DOS of H , and if the Fermi 
projection is Sobolev difjerentiable, the following formulcB hold 

PpCj^'m (1 - Pf) = -thPFd^Ppil - Pf), 

(1 - PF)C-H^d,HPF = lh{l - PF)d,PFPF, 

Proof. Let us consider the first formula only, because the other can be treated in the same 
way (notice however the change of sign). Let -B+ be the right hand side. Then 

h£H{B+) = PF[H,d,PF]{l - Pf). 

Since H commutes with Pf, we find [iJ, SjPp] = — [(9j//, Pp] (after using basic properties of 
derivations). This gives immediately: 

hCH{B+)=PFd,H{l~PF). 

Since 5+ connects only energies below the Fermi level with energies above it, and since the 
Fermi level is not an eigenvalue of H, Ch is invertible on the subspace of such operators, giving 
the result of the lemma. □ 

Corollary 2 (IQHE-Kubo Formula) // the Fermi level is not a discontinuity point of the 
DOS of H , in the zero temperature and infinite relaxation time limit and provided the Fermi 
projection is Sobolev difjerentiable, namely VPf € L'^{A,T), the conductivity tensor is given 
by 

a,, = ^2mr{PF[d,PF,d,PF]) • 
In particular the direct conductivity vanishes. 



4.3 Estimating the deviations from the IQHE Umit 

Before returning to the IQHE, we want to give some idea of the accuracy of the IQHE-Kubo 
formula given in the Corollary |^. We recall that it is obtained under the following conditions: 

(i) the area of the sample is considered as infinite; 

(ii) we work within the relaxation time approximation; 

(ill) the electric field is vanishingly small; 

(iv) the temperature is zero; 

(v) the relaxation time is infinite. 
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In this section we will evaluate roughly the size of the correction terms whenever one of these 
conditions is relaxed. We know that the relative error on the IQHE measurement of the 
universal constant e^ /h is of order of 2 x 10~® at best PD|. Therefore the neglected terms 



should be smaller than this number in order that the experiment is reliable. 

We will not estimate seriously the finite-size effects even though they are in principle ac- 
cessible to a mathematical estimate within the non-commutative framework. However, it is 
generally accepted that these effects decrease to zero as e~^/^ where L is the sample size and 
r some typical length. We will choose r to be of the order of magnitude of the localization 
length. We shall see that it diverges precisely whenever the Hall conductance jumps from one 
plateau to another. Then as this localization length increases, there is a critical value beyond 
which the Coulomb interaction between electrons can no longer be neglected. This is the situ- 
ation in which the FQHE occurs. To estimate this value, one can define r to be such that the 
electrostatic energy of a pair of electrons separated by a distance r is of the order of magnitude 
of the cyclotron energy. This gives r k. Ifim namely an overestimated large distance compared 
to the magnetic length. Thus as soon as L > 20/im, the finite size effects are negligible on the 
integer plateaux. Fluctuations from sample to sample due to finite size and disorder will then 
be negligible. Let us indicate that breakdown of the IQHE due to finite size effects has been 



observed |23|1 . 



The relaxation time approximation should actually be enough to estimate other effects. 
Corrections to such an approximation should not be effective at zero temperature, since the 
relaxation time depends almost only upon the energy level. Only if transitions between different 
bands are taken into account, it is necessary to go beyond this approximation. This problem 
is too difficult to be investigated here. 

Let us consider the effect of a non-zero electric field. Non-linear effects such as bistability 
or hysteresis have indeed been observed in such devices [^. If formula ( ^OD is correct, the 
correction term due to a non-zero electric field is given by 

Sa, , = (h^ < dJsJHM---^ "^^i^ n -/ diH > . 

The relative error is measured by the size of the ratio qS-V /hCn between the electric energy 
and the level separation. If we choose a level separation of the order of the distance between 
Landau levels, namely huc, and let the electric energy be of the order of eSa (where a is the 
atomic distance, namely lA), we find Safer ^ 10~^ for an electric field of IV/m. Thus non- 
linear effects require higher electric fields and it is very easy to choose the electric field in such 
a way as to make this correction negligible. 

Non-zero temperature effects at infinte relaxation time can be estimated as follows. Coming 
back to eq. (^), using eq. ( p]) and the Lemma ^ we can write the Kubo formula at finite 
temperature as 

dE 4^(E)^2m T {P^E[diP<E, d2P<E]) , (43) 

provided we assume the relaxation time to be infinite. Here P<e is the eigenprojection of the 
Hamiltonian on energies smaller than or equal to E. The spectral theorem allowed us to write 
(with an integration by parts): 

/OO /-OO 

dP<E f{E) = - dE f\E) P^E. 

-OO J —OO 
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whenever / is a bounded different iable function vanishing rapidly at +00. We will show that 
the integrand is quantized and equal to ne^ /h whenever hudn — 1/2) < E < hujciji + 1/2) 
with n = 0,1,2,.... This is actually true for the most common devices. If they have charged 
impurities the jump does not necessarily occur at the Landau level [^. But this is of no 
importance for the present discussion since we only need orders of magnitude for the correction 



terms. Incidentally eq. (^3]) shows that the direct conductivity vanishes at finite temperature 
as well (if the relaxation time is infinite). The error term can easily be computed and is given 

by 



^H n'<0,n'^n ^ J^^{n>-l/2) 4 COsh {(3{E - /i)/2) 

provided hujc{n — 1/2) < Ep < %ujc{n + 1/2). We know that the chemical potential equals Ep 
up to small terms of order T^ (where T is the temperature). If Ep is located at a distance 
rfiujc from the nearest Landau level, where < r < 1/2, this correction is thus of the order of 
g-r/3aa;c_ p^j, ^ magnetic field of lOT and a charge carrier mass of .Irrig, Tiojc/kB ~ 140°i^. This 
term is smaller than 10^^ for temperature lower than 4°K. Accurate experiments are often 
performed at 50mK and magnetic fields can be as high as 18T. Moreover, in heterojunctions, 
the effective mass of the charge carriers is one order of magnitude smaller. Therefore, the pure 
temperature effect (at infinite relaxation time) is indeed negligible. 

The infinite relaxation time approximation is in fact the most important effect. The relevant 
parameter is the relaxation time Trei = t/(1 — n)', it can be estimated by means of Drude's 
formula by measuring conductivities. Assuming the relaxation time to be large, the lowest 
order contribution to the Hall conductivity is given by (see eq. (|4l|) ) 

an ~ ^—r (Pp[CH\d,Pp), d2PF]) . (45) 

Estimating ll,~^ by l/cUc, we see that the error term due to this contribution is controlled by 
e = {tuJc)~^. For semiconductors used in the IQHE, the mobility fic (c denotes the type of charge 
carriers) at zero magnetic field and 4.2°i^ varies from about lO'^cm'^/Vs for the MOSFET to 
lO^cm'^/Vs for the AlGaAS or InAs-GaSb heterojunctions |2^[pp.40& 41]. These high-mobility 
devices are used for the FQHE. On the other hand, the effective mass of charge carriers varies 
from 0.2r7ie for the MOSFET to O.OSrrie for InAs-GaSb heterojunctions. Since tUc = p-c^, we 
find £ Ri 0.1 ~ 0.005 for B = lOT. We remark, however, that our relaxation time Trei only 
incorporates time-dependent disorder and may therefore be significantly smaller than the one 
calculated with Drude's formula. 

It is usually estimated that the deviation of the Hall conductance from its ideal value is 
linked to the direct conductivity cr// by 



a 



II 



ban 

Measurements give minio j j^ < 10~'^max{a//) ^ 10~'^e'^/h [Q so that the relative error on the 
Hall conductance is indeed of the order of 10~^. 

Why is a// so small ? Looking at the equivalent of equation (^5]) for the direct conductivity, 
we see that one way to estimate it consists of replacing Cj^ by l/cOc whereas the remaining 
terms are related to the localization length A. But since T{Pp) = n is the charge carrier density, 
we expect this term to be of the order of n\^ /T^^iuOf.. The mobility /Xc of the charge carriers 
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c is TreiOJc/B. Moieover we introduce the filling factor v = nh/qB leading to the very rough 
estimate 

6aH e A^ 

< const. z/- . 

(Th h He 

For the heterojunction, v is typically between 1 and 10 at most. The ratio q/h is given by the 
electron charge e/h and is thus universal. The localization length is always bigger than 80A. 
The mobility is at most equal to 2-l{]^cm? /Vs for the best heterojunctions. We then find a ratio 
of the order of 10~^. Therefore, this estimate is too crude to explain the high accuracy observed 
in experiment. Nevertheless, it shows that collision effects are dominant whereas localization 
is a necessity. 



In any way, physical arguments indicate [71, 20 1 that because of the small value of the 



density of states between Landau levels, the leading contribution to the direct conductivity is 
given by phonon-assisted hopping, at least while the Fermi energy is not too close to a Landau 
level. Estimation (|4^) only includes conductivity by thermal acctivation, whereas (^5|) only 
incorporates effects due to a finite collision rate. Here, we shall not treat the interplay of the 
two phenomena, but only indicate that Mott theory leads to a// ~ (Tq exp(— (Tq/T)^/^) in two 
dimensions, while including Coulomb interaction would give a// ~ ctq exp(— (Tq/T)^/^) [^]; 



this latter is in better agreement with experiment [71, 73 . 



4.4 Dixmier trace and Sobolev space 

In this section, we introduce the Dixmier trace and we will prove a formula that can be found 



in a similar form in a paper by A. Connes [^. We will see that this tool is a key point both 



in proving the integrality of the Chern number of a projection and in the study of localization. 
Most of the material presented her can be found in [|6^, ^ , so that we will give no proof unless 
absolutely necessary. 

First, let us recall that, given a separable Hilbert space 7i, /C(7i) is the C*-algebra of compact 
operators on 7i, namely the norm-closure of the set of finite rank linear operators on Ti. The 
Schatten ideal D'iTi) is defined as the set of compact operators T such that Z^^i /^n < oo, 
where the /i^'s are the eigenvalues of (TT*)^/"^ labeled in the decreasing order. The following 



proposition summarizes the main properties of Schatten ideals; it can be found in [|7^, ^, for 
example. 

Proposition 5 Let CiTi) he the algebra of hounded operators on Ti and Tr the usual trace on 
C{7i). Then we have the following: 

(i) CP{n) = {Te an) \ Tr([T[P) < oo}. 

(ii) CiTi) is a two-sided ideal in C{7{). 

(iii) CiTi) is a Banach space for the norm \\ T \\^= (Tr([T[^))p = {J2n{.l^n{TY))'p . 

(iv) CPiU) C aiH) forp< q. 

(v) Let p,q,s G [1, oo) with ^ = - + -, S* G D'iTi) and T G C^iTi). Then, Holder's inequality 
holds: \\ ST |[r<|[ S \\p\\ T \\q. 

Now we introduce the Macaev ideals O'^iTi) and CP~{l-i) and the Dixmier trace. 
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Definition 4 Let Ti he a separable Hilbert space and K, he the ideal of compact operators on 
H. For p G [1, oo), the Macaev ideal O'^iTi) <Z JC is the set of compact operators T of which 
the characteristic values satisfy 



N 
n=l 



Ar-.oo In A^ 



where the characteristic values are the eigenvalues (/i„) of {TT*)^ labeled in decreasing order. 
C^~{7i) is defined in the much same way but with the limsup equal to zero. We will also set 



1 ^ 

n=l 



Af->oo In A^ 



Theorem 7 i) O'^ and CP are two-sided ideals in C{7i). 

ii) For p e [1, cx)) one has: C? C C?' C C?^ C £p+" Ve > 0. 

Hi) The expression || T ||p-(. defines a norm on D'^ , making it into a Banach space. 



Next, the Dixmier trace is constructed as follows |^9[. Let Lim be a positive linear functional 
on the space of bounded sequences /^(N) of positive real numbers which is translation and 
scale invariant. If a G /^(N) converges, then the functional Lim satisfies: 

Limfa) = lim a„ . (46) 

Scale invariance means that Lim(a) = Lim(ai, ai, 0:2, a2, • • •)• To construct Lim, Dixmier uses 
an invariant mean on the Euclidean group of R (the existence of such means results from a 
theorem of von Neumann) . 

Definition 5 For positive T G C^^ and a fixed scale-invariant and positive linear functional 
Lim on /^(N) satisfying (^^), the Dixmier trace is defined by: 

1 ^ 
TtuUT) = Lim(7— - J2 /^n) ■ 

"^^^ n=l 

Remark that T G £^+ if and only if TrDi^(|T|) < 00. Moreover, if the sequence (j^ Z^iLiA^") 
converges, then all functionals Lim of the sequence are equal to the limit and the Dixmier trace 
is given by this limit. From this definition, one can show that Tr^ix is a trace in the following 



sense 



Proposition 6 The functional Tr^ix defined in Definition \^ can be extended as a linear form 

on C^^ such that: 

(i) positivity: if T E C^^ is a positive operator, then Tij^i^{T) > 0, 

(a) trace property: if S,T & C^~^ then Tij^i^^ST) = Tij^i^iTS), 

(Hi) unitary invariance: ifT& C^'^ and U is unitary then Ttj^i^IUTU^-^) = TiDi^iT), 

(iv) continuity: it is continuous with respect to the semi-norm \\ T ||i+. Moreover, Tr^x vanishes 

on L}~ . 
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Let us introduce the operator 6 acting on a linear operator A on P{G) as 

SA=lu.A]. u = §±^ 

where Xi,X2 are the components of the position operator. The main result of this section is 
given by the following proposition. 

Proposition 7 Let Q be a compact metrizable space on which G acts by homeomorphisms. Let 
P be a G-invariant ergodic probability on Q. One then denotes by A the C* -algebra of this 
dynamical system and by T the trace on A corresponding to P. Let Aq be the dense subalgebra 
of continuous functions with compact support on Q x G. 
Then for every A G ^O; the following formula holds: 

T{\VA\^) = -Troi.(|M<^n , for P-almost all uo . (47) 

If S denotes the Sobolev space associated to T , this formula can be continued to elements A & S. 
In particular, if A E S, then SA,^ G £^^ for P-almost all u. 

In the remaining part of this section we present the proof of Proposition |^ for the case of 
discrete physical space G = Z^; the continuous case will be treated in future work [^. The 
first step in this proof is the following lemma: 

Lemma 2 Let T be a bounded operator on £^(Z^) such that 

i) 3 r such that < n\T\m >= \f \n — m\ > r. 

ii) There exists a positive constant C such that \ < n\T\m > \ < C/{1 + \n\^) V m G Z^ . 

T is then in the Macaev ideal C'^^ and, for any linear functional Lim, its Dixmier trace can be 
calculated by 

Tro,.(r) = TrDi.(Diag(T)) 
where DiagfTj is the diagonal matrix such that < 'n,|Diag(T)|m >= 5n,m < n\T\m >. 

Proof. Since T has only a finite number of diagonals, it can be written as a finite sum of 
operators having only one non-zero diagonal. Using the additivity of the Dixmier trace there 
is no loss of generality in assuming that T has only one non-zero diagonal, namely that it acts 
on £2(z^) as 

TiP{n) = t{n)ilj{n - a) , ip e f{Z^), a e Z^ , 

where t is a sequence on Z^ such that \t{n)\ < C{1 + |?7,|)~'^. It is thus enough to prove that 
T G C^~^ whatever the value of a and that, if a 7^ 0, its Dixmier trace vanishes. 

It is clear that the modulus \T\ of T is a diagonal operator dominated by CR where R is the 
multiplication operator by (1 + |?2|)^'^. Let us show that R G C^^ which implies T G £^+. Its 
eigenvalues are 1/j^ with a multiplicity 0{j^~^), therefore labeling them in decreasing order 
/ii > . . . > /is > . . . with their multiplicity, we get: 
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I N ^ 0(iVi/0) ^ 

s^P l:7^^^s<sup— V const.- < oo . 

Thus ReC^+. 

Let us now assume a 7^ 0. We will show that T is then unitarily equivalent to —T. Since 
the Dixmier trace is invariant by unitary transformations (Proposition ^) it will follow that 
Troix(T) = 0. We remark that for any n G Z'^ the subspace £"„ = /^(n + aZ) is invariant 
under the action of T. Clearly £^„ is isomorphic to /^(Z) and through this isomorphism T 
acts as Tnf{j) = t{n + ja)ip{j — 1) on ip E l^CZi). Let us define the unitary operator U on 
/^(Z) as the multiplication by (—)■'• Then one easily finds UTnU^^ = —Tn- Lifting U to £n 
gives a unitary denoted by [/„. Now ["^{Z^) is the direct sum of the Sn whenever n runs in a 
fundamental domain of the subgroup aZ acting on Z^ by translation. Taking the direct sum 
of the corresponding U^s gives a unitary U on /^(Z^) such that UTU^^ = —T. □ 

The next step in the proof of Proposition |^ is the following 

Lemma 3 Let T, be a subset of TP not including the origin and with finite density DensfSj, 
namely 

Dens(S) = lim -—7- T^ 1 , 

neS,\n\<N 

where \n\ is the euclidean norm of the vector n. Then if Ry, is the restriction to S of the 
operator of multiplication by l/\n\^ in P{Z^), we find 

TtoURs) = — Dens(S) , 

where ujd is the area of the D — 1 unit sphere o/R^. In particular, for D = 2, the geometrical 
constant in the right hand side is it. 

Proof. The eigenvalues of i? are 1/j^. The multiplicity 5'j(S) of such an eigenvalue is therefore 
given by the number of n's in S such that \n\ = j. Let E^v is the subset of S of elements n 
with \n\ < N. Since the eigenvalues of R are already labeled in decreasing order, we obtain 

TrUR^) = hm -±- J: ^ , (48) 

where \T^n\ denotes the number of points in Ejv. Using the definition of the density of a subset, 
we see that, as j — >• 00, the multiplicity QjiX,) is asymptotically given by the product of the 
density of E by the volume between the balls of radii j — 1 and j . Namely 

^,(E)^'^°°Dens(E)a;Bj^-^ 

In much the same way, we get lE^rl ~°° Dens(E) VLj:, N^ if VLj:, is the volume of the unit ball 
of R^. Taking the logarithm we are left with DlnA^ + 0(1) in the expression of the Dixmier 
trace of R^,. Plugging all these relations in eq. (^8]), we get the result. □ 

The last technical step in the proof of Proposition ^ is provided by the following lemma 
where the dimension is D = 2. 
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Lemma 4 Let f be a continuous non-negative function on Q and a G Z^, a 7^ 0. Let F^ the 
operator on P^ZZ"^) defined by: 

\n\ \n — a\ 
Then F^ G C^^ and its Dixmier trace is given by 

TrDi„(FJJ) = — \a\^ I dP{uj)f{uj) , for P-almost every u E il . 

Proof. As \n\ — > 00 the function 



n n — a 

0W = 1^ 



2 



\n\ \n — a\ 

admits the asymptotics 0(n) ~ |apsin^a„/|np modulo terms of order l/|n|^, where «„ is 
the angle between the directions of a and n. 

Let us now slice the space Q according to the finite partition {Qj s}, where 5 > is small 
enough and j an integer such that ^2^,5 is the set of points u for which (j — 1/2)6 < fiu) < 
(j + 1/2)5. Since / is continuous with compact support, it is bounded so that only a finite 
number of j's are needed here. Let then Y,j^s{ijj) be then set of n's in 7? such that T~'^uj G ^j,s. 
Using Birkhoff's ergodic theorem, for P-almost every a;, Sj_5(a;) has a finite density given by 
the probability 'P{VLjs)- 

We then slice Sj ,5 into a finite subpartition {Sj ,5,,} where S^ ,5^^ corresponds to those points 
n G Ej_5 for which (r — 1/2)5 < cun < ij + 1/2)5. Thus modulo an error of order 0{5) we get 
f(T~^uj)(f){n) = |apj(5 sin^ {r5){l + 0{5))/\n\^ on ^j,s,r- Moreover since this slicing concerns 
only a finite partition it permits to write F^ as a finite sum namely 

F°- = VF"|y^ 

hi" 
It is thus sufficient to compute the Dixmier trace of the restriction to Sj ,5,. of F^. But up to 
an error of order 0[5) this restriction is nothing but |apj5sin^ (r(5)/?s i^- Using the Lemma ^ 
we then get 

'Yx^.XKKJ = ^ \a\'j5sm' {r5){l + 0{5)) Dens(S,-5,.) . 

Due to the slicing of the angles one gets Dens(Sj,5,,) = 5/27rDens(Sj 5) = P{Qjs)5/2'k. Plugging 
everything together, summing up over j, r and letting 5 converge to zero, the sum over r gives 
the averaged value of sin^(a), namely 1/2, whereas the sum over j reconstructs the integral of 
/• n 

Proof of Proposition |^ (end). Thanks to Lemma 0, it is enough to compute the diagonal 
elements of I^A^^p because A E Ao so that the number of non-zero diagonals is finite. The 
diagonal elements are 

<n\\5AJ^\n>= Y. \A{T~^uj,a)\^\^ - ,^^^\\ 



^r,o \n\ \n — a\ 

aGZi-' III I 



The number of terms in this sum is finite. Using Lemma H we find 
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Tro..(|Mj2) = ^ ^ |a|2 / dP{u)\A{u,a)\' , 

for P-almost all cj's. On the other hand, the definition of the differential on A gives (see 
eq. (|33D) VA{u!,a) = zdA{uj,a). In particular 

\VAWuj,0)= Y. \a\^\A{uj,a)\^ . 

aGZ2 

To get the trace per unit volume, we just have to integrate both sides of this equation giving 
the Connes formula. 

Since the left hand side of the Connes formula is the dominated by the Sobolev norm of 
A, one can extends this formula to A E S. In particular the finiteness of the right hand side 
implies that SA^^ G C^^ P-almost surely. D 

4.5 Non-commutative Chern character 

We denote by V{A) the set of orthogonal projections in the C*-algebra A, namely the set of 
elements P in ^ such that P = P^ = P*. If in addition P is differentiable, we define its Chern 
character as 

Ch(P) = 2z7r T{P[diP, diP]) . (49) 

If we work on a lattice, T is normalized such that T(l) = 1. In the continuum case, we 
will normalize it in reference to the projections onto Landau levels, namely the eigenvalues of 
the Landau Hamiltonian. (^. An elementary calculation gives for the lowest Landau level 
projection the following integral kernel 

no(f , y) = ^-^e-^^^^-y^'-^^^^^y . (50) 

In particular, IIq defines an element of ^ = C*{Q, R^, B) for any choice of the hull Q. In much 
the same way, we denote by n„ the projection on the n*^ Landau level {n = 0, 1, 2, . . .). We 
then deduce the following results. 

Lemma 5 The trace and the Chern character of the Landau levels are given by 

r(n„) = ^ , ch(n„) = -i . (51) 

Proof. Let us prove this result for IIq first. Its trace per unit volume is given by the space 
average of Ilo{x,x) = qB/h. This gives the first formula in (|5TD. To compute its Chern 
character, we remark that Vno(x, ^) = — z(x — y)Ilo{x,y). Introducing the complex variables 

X = xi + 1x2 and y = yi + iy2, we get 

Ch(no) = nA' f d'x d^y e-ii(l^l'+l2'l'-^^)(xy - yx) . 

It J L-" X \-j 

To compute this integral, we develop the exponential in powers of xy and notice that all 
contributions vanish except the one involving the term |xp|?/p. The corresponding integral has 
separated variables and can be computed explicitly. This gives —1. 
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For the other Landau levels, we remark that n„ ~ Ho in the sense of von Neumann equiv- 
alence (see below). More concretely, it is enough to exhibit an element Un E A such that 
Un = U*Un and Ho = UnU*. This implies that the traces are identical and we will show their 
Chern characters to be identical as well (see Lemma |lT| below). To construct Un, we introduce 
the annihilation operator 

a = (Pi - qAi + iP2 - qA2)/^J2hqB. 

Then [a^a*] = 1. Thus aa* is bounded below by 1 and is invertible. We set u = {aa*)~^^'^a 
and Un = Hqu"'. It is easy to check that uu* = 1 implying that Hq = UnU*. On the other 
hand, U*Un = a*"n„a"/n!. But it is a standard result that we obtain the n*'^ Landau level by 
applying the creation operator n times to the groundstate, namely \n >= {1/Vn\)a*'^\0 >. In 
particular a*"n„a"/ra! = 11^. 

It remains to show that Un & A. A straightforward calculation gives Un = {n\)~^^'^'\loa"' . Now 
using the explicit form of the matrix elements (|50|) and of a, we get the matrix elements of 
Un in the form of a polynomial in x and y times Ilo{x, y), showing that f/„(0, x) is absolutely 
summable in x. This is enough to show that it belongs to A (see Section |3.6| ). □ 

Our next step will be the von Neumann equivalence. Namely ii P,Q E V{A), then P r^ Q 
if there isU & A such that P = U*U,Q = UU*. In particular if P is trace class, it follows that 
T{P) = T{Q) then. The following results can be found in | 78[ 



Lemma 6 Let P,Q e V{A) be such that || P - Q ||< 1. Then P ^ Q. 



Lemma 7 If A is separable, then the set of equivalence classes of projections in A is at most 
countable. 



Lemma 8 Let P,Q & '^(-^) ^^ two mutually orthogonal projections (namely PQ = QP = 0). 
Then the equivalence class of their sum P ® Q depends only upon the equivalence classes of P 
and Q. This defines a commutative and associative composition law on the set of equivalence 
classes, which we denote by [P Q) Q] = [P] + [Q] . 



See [|79[ and [^T], Lemma 4.2.3] for a proof. 

Remark that we can add only mutually orthogonal projections in this way, because the sum 
of two projections is not a projection in general. So giving any pai P, Q of projections it is 
not always possible to find equivalent projections P' ^ P and Q' ^ Q such that P' and Q' be 
orthogonal. In other words the sum is not everywhere defined in V{A). In order to deal with 
this problem, we replace the algebra A by its stabilization, namely the tensor product A® K 
with the algebra of compact operators on a separable Hilbert space, which is nothing but the 
smallest C*-algebra containing all finite dimensional matrices. Then its possible to show that 
one can always choose pairs of projection as orthogonal up to equivalence. By the Grothendieck 
method one builts a group out of the equivalence classes of projections oi A®1C. This group 
is denoted Kq{A) (see []7U[ and [^ Theorem 10]). 



Lemma 9 For any P G 'P{A) there is a differentiable projection Pq G 'P{A) such that P ~ Pq. 

This last result is a consequence of the fact that the set of differential elements in A is norm 
dense. 
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Lemma 10 For any pair P,Q & 'Pi-^) '^^^ ^'^y ^ small enough, there are differentiable pro- 
jections Pe, Qe G V{A) and a differentiable element U^ such that \\ P^ — P \\< e, || Qe — Q ||< £ 
and Pe = UsrU* whereas Qe = U*Us. 

The proof is straightforward: it follows from the density of C^{A) and from the proof of 
Lemma p. We will say that P and Q are smoothly equivalent whenever the element U which 
connects them can be chosen differentiable. 

Lemma 11 For any pair P,Q ^ "^i^) of smoothly equivalent projections Ch(P) = Ch.{Q). 

The proof of this result is purely combinatorial provided we use the cyclicity of the trace T. 
We will omit it here |Q. 

Lemma 12 For any pair P,Q E V{A) of mutually orthogonal smooth projections we have 

ch(p©g) = ai(P) + ch(g). 

The proof of this result is standard and can be found in [Q for instance. To summarize this set 
of results we have 

Theorem 8 The Chern character Ch defines a group homomorphism from Kq{A) into a count- 
able subgroup of the real line. 

It remains to show that the image of this map is the set of integers. This will be done in the 



following Sections [4.6| and 4.7 



4.6 Connes formulae 

In order to compute eventually the Chern character of a projection, we need an intermediate 
tool, namely a cyclic cocycle. This is actually the heart of Connes work on the non- commutative 
extension of cohomology. Here, we shall actually only need a 2-cocycle T2, which is a trilinear 
form on the algebra Ao defined by: 

%iAo, Ai, A2) = 2m T{A^diA^d2A2 - Ao^aAi^i^) . (52) 

By the Schwarz inequality for the trace r, we see that T2 can be extended to the non- 
commutative Sobolev space iS, which is linear subspace of the von Neumann algebra W. The 
proof of the following lemma is algebraic and standard by now; it can be found in [Q . 

Lemma 13 T2 is a 2-cocycle, i.e. it satisfies the following algebraic properties: 
i) T2 IS cyclic: T2{Ao, Ai, A2) = T2{A2,Ao,Ai) 
a) T2 is closed under Hochschild's boundary operator, that is: b: 

{bT2)iAo,A,,A2,A,) 
= T2{AoAi, A2, A3) - T2iAo, A1A2, As) + T2{Ao, A^, A2A3) - T2{A3Ao, Ai, A2) = . 
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We will now give a formula which permits to compute the cocycle T2 by means of the the 
physical representations. For this purpose, we present the formalism introduced by Connes 
(we already gave some indications in Section p. 21 ). A graded Fredholm module is defined as 
follows. Let H+ and ?i_ be two separable Hilbert spaces; their direct sum H = 7Y+ © H^ 
becomes a graded Hilbert space through the graduation operator G equal to ±1 on 7i±. A 
representation n : V —^ ^0^) of an algebra V is said to be trivially graded, if [:n-(74),G] = 
for all A ^ v. One says that ti{A) is of degree 0; operators on Ti which anticommute with 
the graduation operator G are said to be of degree 1. Any operator on Id can be uniquely 
decomposed into the sum of an operator of degree with an operator of degree 1. 

Definition 6 A Fredholm module is a fam,ily ("D, 7i, tt, F) where T> is an algebra with a trivially 
graded representation tc on a graded, separable Hilbert space Ti and where F G C{l-C) is a 
selfadjoint operator such that the following three conditions are satisfied: 

i) FG = -GF li) F^ = 1 izi) [Tf{A), F] e K. VA e I? . 

Here K. is the ideal of compact operators on 7i. An element A eV is called p-summable (resp. 
p+-summable) whenever [ti{A)^F\ G C^{7{) (resp. G C^~^{7i)). 

The graded commutator of two graded operators T, T' G C(H) is defined by 

\T,T']s = TT' - [-iY°^'^)'^°^'^')t'T . 

This commutator extends to the whole algebra C{li.) by bilinearity. The non- commutative 
differential is given by 

dT=[F,T]s TeC{H). 

One can check that it obeys the graded Leibniz rule d{TT') = dTT' + {—Y°^'^'TdT' and that 
(i^ = . Finally, the graded trace or the supertrace is defined by 

Tr5(T) = ^Tr(G'F(ir) , 

whenever the right-hand side is well defined. Here, Tr is the usual trace in C{7i). Remark that 
Tr^ is linear and satisfies: 

Tr5(TT') = (-l)'^°(^')^°(^)Tr5(T'T) . 

Moreover, if c/°(T) = 1, then Tr5(r) = 0. 

Now we shall consider the concrete family of Fredholm modules which will be of interest 
to us. The algebra will be ^o- The Hilbert space is H = H^ ® H' = L'^{G) © L'^{G) (where 
G = R^ or G = Z^) and for A G ^0 the representation is given by 



'^Lu[Aj — Ai^ ~ 1 n A ] ' '^ ~ ^^^ -^ ' 

where vr^ is the family of representations defined in (|3lD . Next 

Xi + 1X2 





u 



\Xi + lX2\ 
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here Xi,X2 are the two components of the position operator on L^{G). 

The following result holds for G = Z^ or R^. However, we will prove it only in the discrete 
case and leave the continuous case for a forthcoming work ||7^ . 

Theorem 9 (First Connes formula) The Fredholm module {AojTijTfi^jF) defined above is 
2+-summable for Y* -almost all uj 's. Moreover for every A G ^o; the following formula holds: 

r(|VAH = -Tr^i^drfA^n , for P-almost all u . (53) 

vr 

IfS denotes the Sobolev space associated to T , this formula can be continued to elements A E S. 
In particular, if A E S, then dA^ G £^+ for P-almost all co. 

This result is an elementary extension of Proposition |^. Its proof is left to the reader. 

The following formula links the cocycle T2 defined in (|5^) to the previous Fredholm modules 
{AoyHyTTi^jF); it can already be found in 0]. This is a result specific to the algebra ^0 and 
depends upon the dimension D = 2 oi the physical space (or of the Brillouin zone). 

Theorem 10 (Second Connes formula) For Aq,Ai,A2 G Aq, we have the following for- 
mula: 

f dP{uj)TTs{Ao,^dAi,^dA2,^) = T2(Ao, Ai, A2) . (54) 

Jn 

Proof. We remark that the left hand side is well defined thanks to Theorem ^. For indeed, if 
A E Aq then dA^ G C^^ C C^ P-almost surely. The trace of an integral operator on L^{G) with 
continuous and compactly supported integral kernel is given by the integral of its diagonal. We 
can therefore evaluate the left hand side: 



dP{uj)TTs{Ao^^dAi^^dA2,u 



dPiuj) / d'^xo d^xi d'^X2 



_n _ ^0 ^1 \/-| _ ^i ^2 w. _ X2 xq , 

|xo| |a;i| |a;i| \x2\ \x2\ \xq\ 



(55) 

where A = qB/2h. The main ingredient of the proof is now the following lemma for which there 
are two different proofs in 0] and Jl^. We shall follow ||l^, but present a discrete version of 



the proof; for the continuous case G = R^ we refer to [[l7 |. 



Lemma 14 Let a,b E G which we write as a = ai + ia2, b = bi + tb2. Then we have: 



, , s s — a . . s — a s — b . . s — b s 

-2m aAb= (1 - —- r)(l - -)(1 - 



/seG |s| |s — a| |s — a||s — 6| |s — 6| |s| 

Proof. In the discrete case G = Z^, the integral appearing in the lemma is in fact a sum, let us 
denote it by G{a,b). Then G{a,b) = —G{b,a) and G{a,b) = —G{a,b). Hence G{a,b) is purely 
imaginary. Now we define: 



48 



e(s,t) = i-riTn - TTTrr) = -e(^, -5) = -e(s,t) . 

\s\ \t\ \t\ \s\ 
A direct calculation leads to 

C{a, &) = — X! (^('^ — a, s — b) + e{s — b,s) + e(s, s — a)) . 
sez2 

We introduce CAr(a,6) as the same sum in which s is restricted to be smaller or equal to 
A/^ e N. The finite difference operators are defined with help of the translation operators Ti, T2 
on functions on Z^ as: 

A,=T,-1 A, = l-i;ri ^- = 1,2. 
We consider the finite differences of of CN{a, b): 

{Aa,Ab,-Ab,Aa,)CN{a,b) 



v^x s — a ~ s — b ~ s — b r s — a . . 

= - E A^il -As,, Ti-^s-r ^A,,- + (1^2). 

I ~r.. s — a \s — b\ \s — b\ \s — a\ 

\s\<N I I I I I I I I 

A discrete analogue of Stokes' theorem allows us to transform the sum over the square into a 
sum over the border of the square. As A^ — ^ 00, this sum converges to the Riemann integral 
/q'^ e*'^(ie~"^. The term (1 ^^ 2) gives the same contribution and we obtain: 

{AaAb, - ^bAa,) C{a,b) = -Am . 
As this is true for every a,b & Z^, C{a, b) is of the form 

C{a, b) = a + /3(a, b) — 27ii a Ab , 

where a is a constant and (3 is linear in the vector (a, 6). As C(0, 0) = 0, we have a = 0. 
Because C{a, b) and a /\b are both odd under permutation of a and 6, (3 must be odd as well 
so that it vanishes. □ 

Proof of Theorem |T^ (end). We set s = Xq, Sq = Xq — Xi, Si = Xi — X2 in equation (^) 
and use the invariance of the measure P in order to replace T~''^uj by u. Applying Lemma |1^ 

we get 



dP{uj)TTs{AoJA^JA2j 
n 

The right hand side is precisely the formula for '/^(Aq, Ai, A2). □ 
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4.7 Chern character and Fredholm index 

The main interest of Connes' theory of non-commutative Fredholm modules is the fol- 
lowing in our context: for a given 3-summable projection P in an algebra A, the expression 
a{P, P, P) = ^TisiPdPdP) can be related to the index of a Fredholm operator. In order to 
make this article self-contained, we will reproduce here the main steps relevant for us. 



First we need the following formula due to Fedosov ||80|. A proof can be found in the 
appendix of ^. 

Proposition 8 (Fedosov's formula) Let F be a bounded operators on a Hilbert spaceTi. We 
suppose that (1 - F*F) e Cp{H) and (1 - FF*) e Cp{H) for some p E [1, oo). Then F is a 
Fredholm operator and for every integer n > p its index satisfies 

Ind(F) = Tr((l - F*F)") - Tr((l - FF*)") . 

Proposition 9 Let {V, Ti., tt, F) be a Fredholm module and P E V be a 3-summable projection. 
Then Fp~ = PF+^|p-^- is a Fredholm operator and 

Ind(F+-) = TrsiPdPdP) , P = (^ ""^p ^ °^^ ^ . 

Proof. Suppose Ti'^ = Ti^ for simplicity. Then 

-P[F,PfP = P-PfPfP 

( (l-F+-Fp+)|p^+ \ 

^ (i-Fp+F+-)|p^- ; • 

By hypothesis [F,P] G C^{H). H51der's inequalities imply (1 - F/-Fp+)|p^+ G £^(7^+) and 
(1 - Fp+F^-)|p^- G C\n ). By Fedosov's formula we get 

Ind(F+-) = Trp^-((l-Fp+F+-)2)-Trp„+((l-F+-Fp+)2) 

= -Tt^{G{P - Pppppf) . 

We can check that this is equal to TTs{PdPdP) by using the following algebraic identities 

[F,P] = P[F,P] + [F,P]P , P[F,Pf = [F,PfP, 

pyp^Pfn+i ^ -[F,Pf''+^F , FG = -GF . 

The proof can easily be completed. □ 

We shall now extend this result to stochastic operators. 



Theorem 11 ( ||16|] ) Let P E W be a projection belonging to the non- commutative Sobolev 
space S. Then for P-almost every uj eVL, P is 2+-summable and 

Ch(P)=Ind(P^w|p„(w-)) , 
where P^ = n^{P) and -u = t^. In particular, Ch(P) is an integer. 
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Proof. Using Theorem ^, Theorem ^ and Proposition ^, it just remains to show that the 
index is P-almost surely independent of uj. Using the ergodicity of P, it is enough to show 
that the index is translation invariant. For indeed, translating P^ by a G G just changes uj into 
T~°-uj. On the other hand, translating m by a changes it into u + 0(1/|X|). Thus P^^u\p^(J^-^| 
is changed into Pt-'^u^uIp _„ (u-) modulo a compact operator. Since a compact perturbation 
does not change the index of a Fredholm operator, the result is achieved. □ 

Remark. The above theorem is true independently of the choice of the probability measure 
P. However, changing P is equivalent to changing the disorder. Therefore we cannot expect 
the Sobolev condition to hold independently of P. In particular, if if is a given bounded 
selfadjoint operator, the spectrum of its representative n^^{H) is P-almost surely constant [pl|] , 
but changing P may change it. Therefore, P has a physical content. We will see in the 
Section |6.3| some of the consequences of changing the disorder. D 



4.8 Quantization, Fredholm and relative index 

In this section, we will discuss the links between the Laughlin argument as presented in Sec- 



tion 2.5 and our approach. The essential ingredient for that will be the relative index of two 



projections as defined by Avron, Seller and Simon |jT^. It turns out that the singular gauge 
transformation of Laughlin corresponds to the unitary operator u = X/\X\. Thus the charge 
transported after changing the fiux by one quantum is then exactly given by the index we 



computed in Theorem |Tl|. The main improvement upon Laughlin's argument is that we control 
completely the effect of the disorder now, since this index is a topological invariant. Let us 
remark that the topology we are talking about is the one of the Brillouin zone and not of the 
sample as can be erroneously derived from a superficial understanding of Laughlin's argument. 



More precisely, following Section |2.5| , let the varying fiux be 0(t) = ht/er with r so large as 



to produce an adiabatic change. Then at time t = t the new Hamiltonian in (|]) is given by 

HBir) = uHb{0)u* , w = ^ , (56) 

I I 

namely the phase of the wave function changes by e*^ where 6 is the polar angle of the position 
X. Formally, one gets an equivalent formula at each intermediate times, but the domains change 
with time. The eq. (|56[) implies for the corresponding Fermi projections: 



Pf{t) = uPf{0)u* , Ppit) = x<EAHB{t)) . (57) 

The charge transported to infinity after this adiabatic change is the number of states in Pf{t) 
that are not in Pp(0). Since both projections are infinite dimensional, we must be careful in 



computing this number. It is the purpose of the following definition |]T^ to take care of this 



difiiculty. Namely, given two projections P and Q on a Hilbert space we set 

Index(P, Q) = dim(Ker(P - Q - 1)) - dim(Ker(Q - P - 1)) , 
whenever the right side is well defined. Then (see [|l7l for a proof) 

Proposition 10 Let P and Q be two projections on Hilbert space. If {P — Q) G O* , then for 
every integer m such that 2m + 1 > p, the relative index can be calculated as 

Index(P, Q) = Tr((P - g)2™+i) . 
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We then apply this forinahsm to the 2+-sunimable Fermi projections Pf,w and uPf,loU* . It is 
an immediate consequence of the proof of Proposition ^ that: 

Index(Pp,^, uPf,uU*) = -TrsipF,udpF,udpF,co) , (58) 

which is P-almost surely equal to — Ch(Pp). The suitable p = 2m + 1 is p = 3 here. We remark 
that it is the smallest possible one giving a non-zero relative index because, if Pf — uPfu* were 
trace class, the relative index would vanish. Hence the Chern character of the Fermi projection 
can be identified up to sign, with the charged transported at infinity during the adiabatic flux 
change. Let us then call this index the charge deficiency index. 

Owing to the stability of the Fredholm index, if one pierces the flux tube at some other 
place than the origin and then uses another unitary operator than u = e*^, the Fredholm index 
will not change. Moreover, adding some disorder potential to the Landau Hamiltonian will not 
change the index as well as long as the Fermi projection belongs to the Sobolev space (we will 
discuss this condition more precisely in the next section). 



5 Localization and non- commutative Sobolev space 

In this chapter, we relate the Sobolev condition on Fermi projections to the Anderson localiza- 
tion. We will give several mathematical tools to describe rigorously localization in terms of our 



formalism. The main results are Theorems |T3| and |T^ below. As a consequence, the remaining 
results in the main Theorem |I] follow. 



In Section 5.1 we review the well-known Anderson- Past ur criterion for localization. Then 



we give a new definition of the localization length and discuss its relation to other notions of 



localization. Notice, however, that a good part of Section 5.2 has already been published in 



81 1 . Our non-commutative localization length allows us to formulate a mathematically precise 



sufficient condition for the existence of the plateaux. All the tools we introduce fit well within 
the non-commutative framework developed so far. For technical simplicity, however, we will 
restrict ourselves to the lattice case namely for G = TP (tight binding representation) in this 
chapter. Even though we believe that most of these results hold for the continuum case as well, 
the proofs are more difficult and will be postponed to a future work W\ . 



Our first definition of the localization length starts by demanding that the mean square 
distance that a particle moves from a given point in an infinite amount of time be finite. We 
will give a precise definition of the word "mean" we use here. In Theorem |T3|, we show that 
provided this condition holds, the localization length can be defined as a L^(R, (iA/") -function of 
the energy, where M is the DOS of the Hamiltonian under study. Furthermore we give another 
definition of the localization length, based upon the Sobolev norm of the eigenprojections, and 
show in Proposition |TB| that it is equivalent to the first one. Then Theorem |T3| shows that 
if this localization length is finite in some energy interval, the spectrum of the Hamiltonian 
in this interval is pure point for almost all disorder configurations. Theorem |1^ shows that 
under the same condition, the spectral projection Pf on energies lower than or equal to E' is a 
continuous function of E with respect to the the Sobolev norm. This technical result, together 
with the results of the previous chapter, implies that the Chern number of Pg is constant on 
that interval. This is the reason why we get the plateaux of the Hall conductance. 

Let us remark that our localization length can be computed as a disorder average of a 
product between two Green functions, or equivalently by means of a current-current correlation. 
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In particular, let m be the positive measure on R^ defined by [Q 

/ dm{E, E') f{E)g{E') = T {V H f {H)V H g{H)) , f,g e Co(R) . 

Using spectral theory in Li^^AiT), one can indeed show that such a measure exists and that 
it can be calculated using Green's functions. Then the localization condition that we define 
below is given by 

^.(^) o / dm{E,E') 



/axr {E - E'Y ' 

whenever the integral exists. We will not develop this point of view here, but we emphasize 
that our approach is equivalent to the one used in solid state physics. 

5.1 The Anderson-Pastur localization criterion 

Most results in this section are due to Pastur [^, |3^. The underlying physical idea can be 



traced back to Anderson p8[. Since this theory holds in any dimension, we will assume that 
the lattice in space is Z^. Let then fi be a compact space endowed with an action of Z"^ by 
homeomorphisms. The magnetic field B in Z)-dimension is an antisymmetric bilinear form on 



Z^ written as Ba A h. We define the C*-algebra A= C*{Vtyi Z^, B) as in Section |]6[ In this 



discrete case, this algebra has a unit. Given a Z^-invariant ergodic probability measure P on 
n, we get a trace T on A, which is actually normalized. We will denote by W the von Neumann 
algebra L°° {A, T) of the corresponding GNS representation. Recall that the representations tt^ 
extend to W and give random operators. Then, to avoid inessential difficulties, the Hamiltonian 
H will be a selfadjoint element of A in this section. We will denote by cr(co') the spectrum of 
the operator Ti^{H) = H^^. As a preliminary, let us recall the 

Lemma 15 (Wiener criterion) Let ^ be a finite complex measure on the real line, i.e. fi is a 
linear combination of four finite positive measures fj, = fj,i—fj,2+if^3—if^4. Let E^(t) = J e'^^^dfiit) 
be its Fourier transform, then: 



hm r ^\FM' = T. HE)\' 



BeR 

A proof for a positive measure can be found in [|7B[; it can be completed to a proof of Lemma |TB 
without any difficulties. 

The next result is given by 

Proposition 11 Let P be a projection in W, the von Neumann algebra L°°{A,T) generated 
by A in the GNS representation of the trace per unit volume T associated to the probability 
measure P. Then P-almost surely, 7r^(P) is a projection on £^(Z^). Moreover its dimension 
dim.{TTi^{P)) = Tr(7r^(P)) is either zero or infinity P-almost surely. 

Proof. Because of the covariance of P, diin{TT^{P)) is translation invariant. Let us define 
n = J dP {u;)'7i^{P) as a weak integral acting on £^(Z^). By construction, 11 commutes to 
the translation group. Since the measure P is invariant and ergodic, using the monotone 
convergence theorem, we get for P-almost all uj E Q: 
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dim(7r^(P))= /"rfP(cj)Tr(7r<^(P)) = Tr(n) = J2 < a;|n|x >= J2 < 0|n|0 > , 

x&ZO x&ZD 

and this last expression is either zero or infinity. □ 

Let now A C R be a Borel subset of the spectral axis. Let then -P(A) be the corresponding 
spectral projection of H. Let A„(A,a;) be the time-averaged probability for a particle initially 
at site |n >,n G Z^, to stay at this same site \n > the time evolution being governed only by 
the restriction of the Hamiltonian H^^ to the interval A. More precisely 

rT rif 



An{A,u;) = lim / - | < n\7r^{e''''PA)\n > 



The covariance of the Hamiltonian leads to A„(A, u) = Ao{A, T~'^uj). After averaging over the 
disorder configurations, it will therefore be sufficient to consider only A{A,u) = Ao{A,u!). We 
introduce the (disorder average) return probability: 

aA) = JdPiuj)AiA,uj). 
^(A) can also be expressed as follows 

as can be checked by a direct calculation. Here V = (9i, . . . idn) is the derivation on A. If 
^(A) > 0, we expect that there are some localized states corresponding to energies within the 



interval A, this is contained in Theorem |T2| below. Using Wiener's criterion we obtain: 
i{A) = /rfPH hm r^l /<0|7r^(t/PEPA)|0>e'^f 

J T-^oo Jo I J 

= jdP{u) y: i^-,i^(o)r, (59) 

Eeapp{uj)nA 

here Pe G W is the spectral projection of H on the interval {—oo,E], ip^^E is the eigenstate 
of H^^ corresponding to the eigenvalue E and crpp{uj) is the set of eigenvalues of H^. Hence 
we see that ^(A) is related to the so-called "inverse participation ratio". The main result is 
summarized in the following 

Theorem 12 Let A G H be an open interval. Then the following results hold: 

i) Let M he the DOS of H. If X is a growth point of N' then A G cr^uj) P-almost surely. If X 
is not a growth point of JV then X is outside (j{uj) P-almost surely. 

ii) The number of eigenvalues of H^ counted with their multiplicity and contained in A is 
either zero or infinity P-almost surely. 

ill) A given real number X is P-almost never an eigenvalue of H^ with finite multiplicity. 

iv) The number of eigenvalues in A is infinite if and only if ^(A) > 0. 



54 



Proof, i) Suppose A is a growth point of A/". This imphes that < 0\7!'^^{Pa)\0 » on a set of 
positive measure for any neighborhood A of A. Thus, for any such neighborhood, the set of u^s 
such that (J^uj) fl A 7^ has probabihty one, because it is measurable and translation invariant. 
Since 

{u;en\Xe(r{u;)}= f]{u; e n\a{uj) n {X - -, X + -) 7^0} , 
j>o 3 J 

then A G c^u) P-almost surely. 

Conversely let us suppose that A is not a growth point of Af. Then there is a small open 
interval A containing it such that < 0|7rt^(PA)|0 >= on a set Qq of probability one. We set 
^00 = Clnezo T^^Qq to get a translation invariant subset of probability one. On this subset 
< n\TTi^{PA)\n >= for all n G Z^ by covariance. Therefore on this subset, cr{uj) fl A = 0. This 
implies that A is outside cr(a;) with probability one. 

ii) Let Ppp{uj) be the projection onto the subspace spanned by the eigenvectors of H^. It 
defines a covariant, measurable family of bounded operators |^|. So by Theorem ^ it defines an 



element Ppp of W. Thus the number of eigenvalues counted with their multiplicity and contained 
in the Borel set A is the dimension of vr^(-Ppp-PA); namely it is zero or infinity P-almost surely 



by Proposition 11 



iii) A is an eigenvalue of H^^ of finite multiplicity if and only if < Tt{'k^{Px)) < 00. By 



Proposition [TT] this happens with probability zero. 

iv) Let JIa be the set {u E Q \ appiu) fl A 7^ 0}. Thanks to eq. (^) this set has positive 
probability if and only if ^(A) > 0. Since this set is translation invariant and measurable it has 



probability one if and only if ^(A) > 0. By Proposition |TI| again, the result is achieved. □ 



Remark 1: ii) says that, with probability one, no eigenvalue of H^ with finite multiplicity is 
isolated, whereas iii) shows that such eigenvalues are fiuctuating with the disorder. There 
are however examples of models having non-fiuctuating eigenvalues of infinite multiplicity. 

Remark 2: The criterion iv) does not eliminate the occurrence of some continuous spectrum. 



5.2 Non-commutative localization criterion and localization length 

We introduce a second physical idea of localization: the average mean square displacement 
8X{T). Let X be the position operator in G and X^{t) = e^-'^'-'^Xe"*^"* its time evolution. 
Then we consider: 

rT dt 



nl (if 

5X^,„(T)2 = / - < n\{X^{t) - Xf\n > . 
Jo I 



The covariance relation leads to 6X^^n{T) = 6Xt~"^q[T), so that it is sufficient to examine the 
behavior of the state |0 >. Averaging over the disorder and using again our algebra we get 

6X{Tr = £ ^TiMe-'n^X) . (60) 

Boundedness of SX(T) in time T will be an indicator for localization and behavior proportional 
to T'^, a G (0, 1] will be interpreted as diffusive or ballistic quantum motion. In order to localize 
in energy we restrict the motion to energies in a Borel subset A of the real line and get 6X^{T) 
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in much the same way. On the other hand taking T -^ oo leads us to define the A-locahzation 
length as 

f{A) = hmsup r ^r(|V(e-^^*PA)r) • (61) 

Boundedness of /^(A) will be our localization condition for the spectral subset A. 

Theorem 13 Suppose that /^(A) < oo. Then H^ has pure-point spectrum in A for almost 
every uj E Vt. Moreover, if N'{E) = T{Pe) is the density of states, there is a M -measurable 
function I on A such that for every Borel subset A' of A: 

/2(A') = / dU{E) l{Ef . (62) 

J A' 

Finally, if f'{A) < oo, /^(A') is given by 

l\A') = jdV{uj) Y.\n\^ Y. I < ^K{P{E})\n > r . (63) 

neZ2 Eeapp{uj)nA' 

Theorem 14 If the localization condition /^(A) < oo is satisfied, then at every regularity point 
of the density of states, the application E E A ^ Pe E S is continuous. 

Remark 1: Letting A shrink to one point E., the function /^(-E) represents a kind of average 
of the quantity J2n&zD |?^PlV'E,a;(^)p5 where ipE,uj is an eigenvector of H^ corresponding 
to energy E; here the average is taken over the disorder and a small spectral set around 
E. This quantity measures the extension of this eigenstate. For this reason the function 
/ will be called localization length for H. Note that no exponential decay of the wave 
functions is needed for our localization length to be finite. However, such behavior may 
be studied within the present framework; we postpone the details to future work. 

Remark 2: The index theorem which we proved in Section [4.7| only requires that Pe be in S in 
order to insure the integrality of the corresponding Chern character (compare to results 
of Section [4.4| ). Thus boundedness of /^(A) is sufficient but not necessary to prove the 



index theorem. 

Remark 3: The localization condition /^(A) < oo implies pure-point spectrum in A almost 
surely. In mathematical physics, pure-point spectrum of a Hamiltonian if in a certain 
region of the density of states has been considered as a the criterion for localization. The 
eigenstates being square integrable are thus localized. The classical RAGE-theorem [[76| 
permits to make such a statement more accurate. 



Proof of Theorem |T^. i) The basic argument we will use here is due to Guarneri [p^] . 
The spectral projection on the continuous part of the spectrum is then given by vr(^(Pc) = 
1 — ni_j{Ppp). Using the definition of the trace, then applying the theorems of Fubini and 
monotone convergence, we find (with the notation |n|oo = maxi<j<£) \nj\): 



5X^{Tf = fdP{uj) lim V \n\'^PT{uJ,n) 

J N—*oo — ' 



\n\<N 

where 
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rT dt 
/o T 
Then pT{uj,n) satisfies 



PTiuJ,n)= / |<0|7r^(e^^*PA)|n>P 
JO 1 



(64) 



0<PT{uJ,n)<l , Y. PT{oo,n) = < 0|7r^(PA)|0 > < 1. 

nez^ 

Now, we use the spectral theorem for iJ^ and the Wiener criterion 

\impT{uj,n) = hm / —\ < 0\7r^{dPEPA)\n > e'^'\^ 

T^oo T^oo Jo I J 

E \<0\n^{P{E})\n>\\ 

E£a-pp{uj)nA 

In particular, for fixed positive integer L: 

^lim^ ^ PT(w,n) = E Y. < ^\^-^iP{E})\n >< nKiP{E})\0 > 

\n\<L E<^app(u})nA \n\<L 

< < 0\n^{PA)\0 > - < 0\n^{PcPA)\0 > . 
Let us now introduce the following notations: 

Prin) = J dP{uj)pT{uj,n) , 

r = JdP{uj) < 0|7r^(Pei'A)|0 > s = J dP{uj) < 0|7r^(PA)|0 > . 

By the dominated convergence theorem it gives 

/ dP{uj) hm Y, PTi^,^) = lim ^ Prin) < s — r. 

~*°^|n|<L ^°^|n|<L 

Since r > 0, one can find Tl > such that, if T > T^, one has J2\n\<LPT{n) < s — |. Thus, for 
T>Tl: 

5X^{Tf> fdP{u) Y \n\'pT{oo,n) > ^. (65) 

\n\>L ^ 

Taking the lim sup over T we get Lr < P{A) for all L, implying that r = 0. Thus for almost all 
a; e fi: < 0|7r(^(PcPA)|0 >= 0. Using the covariance relation and since Z^ is countable, there 
is r2' C fi of probability one such that < n\7Ti^{PcPA)\n >= for all n G Z^ . In particular 
7r(^(Pc)|n >= for all n G Z^ and all u G Q', namely the continuous spectrum in A is empty. 

ii) Given two Borel subsets Ai, A2 C A, we define the following expression: 

#2(Ai, A2) = / 7^ E 1^1' < 0|7r^(e^^*PAj|n > <0|7r^(e^^*PAj|n> . 
^0 T \n\<L 

This expression gives a Borel function in uj. We use the Wiener criterion: 
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lim42(Ai,A2)= E Nl' E |<0|7r^(P|^})|n>p = 4^)(AinA2) 

' — >oo — — 

|n|<L EecTpp{ui)nAinA2 



From this definition of Sl^\^') for a Borel set A' C A it follows that: 

(a) < E^J^K^') < L2<0|7r^(PppPA')l^> < ^' 

(b) If Ai n A2 = 0, then 4^)(Ai U A2) = ^'^^(Ai) + ^'^^(Aa) 

(c) eHA') < e+^HA') 

(d) Sj^\^') is a Borel function of uj as pointwise limit of Borel functions. 

(e) If (Aj)jgN is a decreasing sequence of Borel subsets such that flj-gN Aj = 0, then Sj^\^j) 

decreases to zero. 

Averaging over the disorder, we obtain £^*^'^)(A') = j dP{uj)S^\^') which fulfills (a),(b),(c), 
and also (e) thanks to the dominated convergence theorem. Moreover, since S'^^\/S!) < L^, we 
can use the dominated convergence theorem, Fubini's theorem and the definition of /^(A) to 
get 



£(^)(A') = lim r^ I dP{uo) V InH <0|7r,(e*^*PA')|n> P < /'(A') . 



rT dt 

I df[uj) 

\n\<L 

In much the same way, we get, thanks to (b), £^^\A') < P{A) < 00 , V L G N. As £(^)(A') is 
bounded and increasing in L, it follows that ^(A') = limi^oo^'-^HA') exists. £ defines a non- 
negative set function on the set of Borel subsets of A. Because of property (e), it is continuous 
from above and since it is moreover finite, we can conclude its a-additivity. Therefore £^ is a 
Radon measure. Moreover, using the monotone convergence theorem: 



^(A') = fdPico) J2 E \n\''\<OKiP{E})\n> 



(iii) Let us show now /^(A') < <^(A'). Recall that: 

t rT df 

S{A') = /rfPMsup Y: H'lim / 7^1 <0| 



n,.W'P^,)\n> P < 00. 



We can replace the limit by the lim sup. Hence for a fixed Tq large enough, the following 
expression is finite and we may apply Fubini's theorem: 

r /"^ dt 

/ dP{uj) sup V |np sup / — I < 0|7r<^(e*-^*PA')|n > ^ 



/f dt 
dP{u) / — sup V InH < 0|7r,;(e*^*PA')|n > ^ 

> lim sup / — / dP{iu) sup V \n\'^\ < 0|7r^(e'^*PA')|« > 

T>To Jo I J ^^'^\n\<N 

= l\A') . 
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As this is true for all Tq, we obtain P{A') < S{A'), and with (ii) their equality. 

(iv) To finish the proof we use the Radon-Nikodym theorem. It is thus sufficient to show 
that the measure /^(A') = S{A') is absolutely continuous with respect to the density of states. 
Let A' C A be such that J\f{A') = T{Pa') = 0. From the definition of the trace it follows 
that < 0\7ii^{PA')\0 >= almost surely. By covariance and because Z^ is countable, this gives 
n^(P^i)\n >= for all n E TP almost surely. Then the definition of £^j.^(A', A') implies that 
it is zero for any /,T and almost all uj] consequently = ^(A') = /^(A'). D 

Let us give another useful expression for /^(A). We consider finite partitions P of A into 
disjoint Borelian subsets: 

p = { Aj C A , Aj Borehan | j = 1 . . . g, |J Aj = A, A^ n A^ = 0} . 

i 

The set 2(A) of such finite partitions is ordered by refinement 

V <V' ^ y A! eV ^A!' eV such that A' c A" . 
This gives an ordered net. We define: 

^''(A)= lim E^IVPa'H, (66) 

where the limit is understood to be the one under refinements of partitions. 
Proposition 12 /^(A) is well defined m R = RU {oo}. Moreover: 

P{A)= sup Y. mVPAf) . 

PeZ{A) A'GP 

Proof. It is sufficient to show that refining the partition results in the increase of the quantity 

Ea'gp ^(IVPa'H- Take A = Ai U A2, Ai n A2 = 0, and suppose TdVPAp) < 00. We have: 

TdVPAp) = T(|VPaJ' + IVPa.T + VPa, ■ VPa, + VPa, ■ VPaJ . 

Now, either one of the Sobolev norms of Pai and Paj is infinite and in this case the inequality 
^(iVPAiuAaP) < ^(IVPaiH + TdVPAaP) is trivially satisfied, or the components of 
VPai and VPa2 are in L'^{A,T). The Holder inequality implies VPa^ ■ VPa^ e L^{A,T). We 
may therefore treat each term separately. 

T(VPa, ■ VPaJ = T(Pa,VPa,(1 - PaJ ■ VPaJ + T((l - PaJVPa.Pa, ■ VPaJ . 
We apply the cyclicity of T and the formulae 

Pai < 1 — Pa2 ) Pa2 < 1 — Pai , VPajPai = — PA2VPA1 , 

to get: 

r(vPA, -vpaJ = r(PA,vPA, •vPAj + r(vPA, -(-vpa.PaJ) 

= -T(|VPa2PaJ')-T(|VPa,Pa.P) < 0. 
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The same calculation for the other term then implies the result. □ 

Because of equation (|63|), one expects /^ to be equal to /^ and this is what we shall prove in 
the sequel. We first need the following technical lemma: 

Lemma 16 Let fi, v he two finite positive measures on a Borel subset A of the real line and 
fipp, Upp their pure-point parts. For a finite partition V of A we set: 

AfiA,V)= 5:MA')KA')- 
A'ev 

Then: 

\im AfiA,V)= inf AfiA,P) = Y. ^pp({i5;})z/,,({E}) . 

Proof. First, we show that A/'(A,P) decreases as the partition is refined. We introduce 
^CP) = Ua'sp A' X A'. Suppose that V < V\ then T{V') C J^(P). As /xOz/ is a finite positive 
measure on a Borel subset of R^, we have: 

M{A, V)=ii® iy{r{V)) >ii® iy{T{V')) = Ar{A, V') . 

Let us now consider the Lebesgue decomposition of the measures /i and v in continuous and 
pure-point parts: 

fi ® U = fipp ® Upp + fie ® ^pp + ^J^pp (S) l^c + fJ-c (^ l^c ■ 

As the measures are finite, the continuous parts satisfy: 



V e > 3S>0 such that /ic(A') < e, u^A') <e V A' with diam(A') < S . 

Here, the diameter is defined by: diam(A') = sup^,^g^/ \x — y\. Now we choose and fix a 
sequence of Vn of finite partitions which satisfies: 

lim ( max diam(A') ) = . 
For such a sequence Vn the following holds: 

fl I U (A'x A')] ={(x,a;)|xG A} = Diag(Ax A) . 

n=l \A'eP„ / 

In the limit, the contribution containing a continuous part vanishes for indeed for instance 



limuc®z/p„ 11 (A' X A') < lim e(n) V z/„„(A') < Mm e(n)iy(A) = 0. 
\A'eP„ / A'ePn 

Therefore, using the a-additivity of the finite measure fi®v (which is equivalent to continuity 
from above), we find: 
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Jim y: -"(aoka') = Jim y: f^pp{A>ppiA') = /ipp®^pj n u (a'x a') 

A'ePn A'ePn \n=i A'ePn 

= f^pp ® z^pp(Diag(A X A) = ^ yUpp({^}) ® '^ppd^}) • 

EeA 

This is true for every such refining sequences of partitions, leading to the result. □ 

Proposition 13 If P{A) < oo, then P{A') = r(A') for every Borelian subset A' C A. 
Proof. As /^(A) < CX3, Theorem O shows that 



£(^)(A) = jdV{u) Y.\n? E I < ^K{P{E})\n > 



\n\<L _EGo-pp(ij)nA 

We decompose the complex measure yu(A') =< 0|7rj^(PA')l^ > by polarisation into four positive 



measures: 



/i(A') = /Xi(A') - /X2(A') + Z/X3(A') - 2^4(A') . 

With these notations, we apply Lemma |16| to each terms to get 

^^^^A) = fdP{uj)Y.H^ E ( E /^fc ® A^fc- 2/^1 ®/^2- 2/13 ®/i4)({^xE}) 

\n\<L E(^cTpp{uj)nA \k=l / 

= fdPiu;) E \^\\}^^„ E I <0|7r.(PA")l^> P ■ 

•' \n\<L '^^(^>A"eV 

The dominated convergence theorem now gives: 



f(^)(A)= hm E /rfPH E \n\' \ < OKiPA.)\n > \' = i'^^^A) , 
I e^(A) ^„^^j 1^1^^ 

by definition of £^^'{A). Obviously we have: 

^(^)(A) < ^lim E /rfPHEl^ri<0|vr.(PA")l^>r = I'i^') ■ 
' '^^(^ > A"ev nez 

Moreover, £^^\A') is increasing in L. Since ^^^•'(A) is bounded by /^(A), so is £^^\A). Its limit 
therefore exists. Actually, it converges to r{A'). For indeed, with the dominated convergence 
theorem, we find: 

P(A') = sup E [ dP{uj) Mm E I^H < 0|^^(^A")I^ > T 

= sup sup E dPi^) E l^n < Ok<^(-PA")l'^ > P 
= lim £^^\A') . 

L— >oo 

This finishes the proof. □ 
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Proof of Theorem [Til. Fix E,E' eA, E < E', then: 



Pe'-Pe\\1 = \\P[E,E']\\l< i\[E,E'])+U{E')-U{E) 



,E' 



= f {l{E"f + l)dM{E") . 

JE 

Now, in the hmit E' -^ E, this is zero if i? is a regularity point of A/". □ 

5.3 Localization in physical models 

In this section we give an example of a physical model for which the localization condition 
Z^(A) < cxD is satisfied both for weak disorder at the band edges and for high disorder all over 
the spectrum. For the mathematical treatment we will, once again, restrict ourselves to the 
discrete case. Our line of arguments will use results of Aizenman and Molchanov |^, ^ . They 



give a simple proof of mathematical results proved earlier [^ |SB|, ^ |^, ^ . We will conclude 
with some remarks about the continuous case. 

As a preliminary, let us remark that if the spectrum has a finite gap A, then the condition 
/^(A) < oo is satisfied for the simple reason that Pa = 0. In dimension two, the Chern character 
corresponding to every energy band is therefore an integer and the Hall conductivity an integer 
multiple of e^/h. However, this integer may be zero. 

The analysis by Aizenman and Molchanov consists of two steps: first, one shows exponential 
decay of low moments of Green's function for concrete classes of models and in specified regions 



A of the spectrum; then, this decay is used pSi to show exponential decay in \n — m\ of 



the quantity / t/P(c<j) supj>o < n\7i{e^^^PA)\m >. This, in turn, will allow us to show that 
/^(A) < oo is satisfied. 

Here, we fix our attention to the Z)-dimensional Anderson model with constant magnetic 
field. The Hamiltonian acting on £^(Z^) is given by: 

Huj,\ = Hq + XV^ a G R , 

where Hq is the Z)-dimensional analog of Hamiltonian given in equation (10). K) is the disorder 
potential: at every site of the lattice it takes a random value in the interval [—1, 1]; the prob- 
ability distribution is supposed to be uniform and the sites are independent. In ^, ^ much 
more general situations are considered: the kinetic part Hq may have non-zero elements farther 
off the diagonal as long as they decay exponentially in distance from the diagonal; the random 
potential may also have gaussian or Cauchy distribution; moreover, correlations between the 
sites are allowed in a sense to be made precise. For us, however, the Anderson model captures 
the essential of the localization phenomenon and it has the advantage that the hull (fi, T, P) 
is easily constructed as topological product of intervals [—1,1]; by Tychonow's theorem, f2 is 
compact; the action T on fi is given by the translations in physical space and the product 
measure of the uniform probability on [—1, 1] is a T-invariant and ergodic measure P on Vt. 

We will now describe precisely in which situations Aizenman and Molchanov prove expo- 
nential decay of low moments of Green's functions for the Anderson model: 

j dV{uj)\<n\'n^{j^^)\m>Y < De"^!""'"! s G (0, 1) . (67) 

The first situation is that of strong disorder: for every s G (0, 1) there exists an Ac = \c{D-, s) 
such that for A > Ac the bound holds on all of the spectral axis. 
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The second situation considers weak disorder, that is small A. Recall that the free Haniilto- 
nian Ho has energy bands of continuous spectrum. In the interesting case for us of dimension 
two, the energy bands and their dependence on the magnetic field are given by Hofstadter's 
butterfly; Hq then gives Harper's equation. For the Anderson model the bound (|67D is obtained 
for energies E situated at band edges: 

E e R{\) = \J {Een- a{Ho) I lxgs{E) < 1} , 

s6(0,l) 

where 




gs{,E)= sup ( Y. I < ^l^^( jj-_^ )l"^> 



(For the Anderson model, the constant k{s) appearing in p5[ is equal to \ independent of 



s.) By a Combes-Thomas argument, one gets [^ the asymptotic behavior of gs{E) in ^ 
dist(E,o-(i/o)): 



S-'^i-^^+T) e ^ 



^^(^)-1 i-\i + oa)) e 



£// V. ■ oo 



Therefore, because (x^Hx) C {E \ E G {a {Ho) + v}]u < A}, the intersection of -R(A) with 
(t{Hx) is non-empty for A sufficiently big. In other words, the bound (|67D, is satisfied for some 
energies belonging to the spectrum and situated at the band edges. This is probably the case 



as soon as A 7^ 0; dimension two is critical [28 



Let us now come to the second step of Aizenman's analysis. If the bound (pTl) is satisfied 
for an interval A of the spectral axis, he proves that the unitary evolution operator filtered 
with the spectral projection Pa satisfies: 

fdP{u}) sup I < n\'K^{PAe'^')\m >\ <C e^^l""™! C,DeR+ . (68) 

Let us now use this result to show that the non-commutative localization condition is satisfied. 
Because 11 ^uiiPA^^^^) II ■, n < 1, we have: 

r(|a(PAe*^*)n = Y. 1^1' /^P(^)l < 0|7r^(PAe^^*)|n > 1^ 

nez^ 

< ^ |n|2 /'c/P(cj)| <0|7r^(PAe'^*)|n> I 
nezo 

< Y. l^rC7e-^l"l < C" < 00 . 

neZ-D 

Now, taking the time-average we see that /^(A) < 00. This is, of course, only true if either 
A > Ac or A c P(A). 

Let us now comment on the case of continuous physical space. First of all, if the Fermi 



energy Ep lies in a finite gap, then it can be shown that Pp G S [^. For strong disorder, 
one expects a finite localization length, but near the set {Tiojcin + \)\uJc = ^, n E Z} the 
localization length probably diverges. This would be the analog of the discrete case. 
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6 Applications and complements 

6.1 Low- lying states do not contribute to the IQHE 

In this section, we address the following question: consider the Landau Hamiltonian Hi and 
add a periodic potential Vp of varying strength 

H = Hl + <tVp, < Vp < 1 . 

How do the Chern numbers evolve as the coupling parameter a is increased? In particular, 
what happens ii a ^ Tiojc 



Halperin et al [£8| made corresponding numerical calculations. They took a finite size 
sample and computed the Chern numbers for various values of a by the method described in 
Section |6.3| . For small a, the weak periodic potential approximation is valid and the Chern 
numbers are those given by the Diophantine equation ([70|) . For intermediate values of a, the 
energy bands cross each other in a complicated way and it is difficult to follow the Chern 
numbers. However, for a ^ hujc-, Halperin et al observed the following striking result: all low 
energy bands of the spectrum carry zero Hall current. This result was put on a rigorous basis 



by Nakamura and one of the authors [37 . 



Theorem 15 Consider the Hamiltonian H = Hl + aV where V is a potential satisfying the 
conditions below. Then for a sufficiently large the spectrum has gaps between his low energy 
bands (energy smaller than ehuJc, where e appears in the conditions below). If the Fermi energy 
lies in one of these gaps, the Hall conductivity vanishes. 

Hypothesis on the potential V: 

i) inf^eR2 V{x) = , sup^.gR2 V{x) <C <oo . 

ii) There is a countable set of R^ , {x„, n = 1, 2, . . .} such that 
\xn — Xm\ > d ii n ^ m with d > 0. 

iii) There are e > and V G C^{B^ where B^ = {x E R^||x| < e} such that d > 2e and 
V{xn + x) = V(x) for X E B^ and for all n 

iv) is the unique non-degenerate minimum of V G -B^ 

v) If \xn — x\ > e for all n, then V{x) > S > 0. 

Note that this potential is not strictly periodic, although the bottoms of the wells need to 
be identical. The framework of the proof is that of Non- Commutative Geometry completed by 
semi-classical analysis. Let us just describe the strategy here. The existence of the gaps follows 
from semiclassical analysis: the energy levels of a single, isolated quantum well are enlarged by 
the tunneling effect; the band width can be estimated by e~"^ , a > 0. In order to show that 
the Chern character is zero, one proceeds as follows. First consider the situation of quantum 
wells separated by infinitely high barriers. The projection P^ on the energy level of quantum 
well number j is one-dimensional. Its Chern character therefore vanishes. Now we consider 
the projection Poo = ®jPe) where the sum runs over all wells. Because of the additivity of 
the Chern character (compare Lemma |12|), the Chern character of P^o is zero as well. Finally, 



semiclassical analysis allows us to show that Pqo is Murray-von Neumann equivalent to the 
projection Pg onto the energy band which arises as the barriers are reduced to finite height. In 
view of Lemma 0, this finishes the proof. It is likely that the result also holds if the projection 
is not in the gaps. 
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6.2 Where and how does the locahzation length diverge? 

In this section, we review the most important results on the behavior of the locahzation length 
at the center of the broadened Landau level. Discussion of theoretical ideas is followed by 
a brief presentation of numerical techniques and results. As before, we consider a system of 
independent fermions in a disorder potential described by a one-particle Hamiltonian. 

In numerical calculations and scaling theory, a rather strong notion of localization is often 
used, that of a Lyapunov exponent X{E). If one supposes exponential decay of the correlation 
function of the wave function ipE corresponding to the energy E, X{E) is defined by: 

<Mr)Mr') > oc e-^(^)l^-'''l . 

Here the average is taken over disorder configurations. Physicists often call the inverse of X{E) 
the localization length. However, except for in the one-dimensional case, there is yet no clear 
mathematical formulation. 

In 1979, the gang of four [^] used scaling theory (ideas due to Thouless |^) and renor- 



malization group calculations to show that, at absence of magnetic field, two is the critical 
dimension for localization in the following sense: for dimension smaller than two, all states are 
exponentially localized at any strength of disorder, that is, for all wavefunction there are posi- 
tive Lyapunov exponents. For dimension greater than two, there exist extended states for low 
disorder. At dimension two, all states are localized except for states corresponding to isolated 
critical energies at which the localization length diverges. In the case of perpendicular magnetic 



field in two dimension, the same result holds as shows a development of higher order |^ . 

In the QHE, the Hall conductivity does not vanish. This led Halperin to postulate the 
existence of extended states |5l| . A corollary of the Theorem |I] is the following: 



Corollary 3 /i5/ Between the plateaux, there has to be a spectral interval A such that the 
localization condition P{A) < oo is violated. 

Remark that the interval may be very small. It is even likely that there is, in fact, a single 
energy E such that P<e is not in the Sobolev space S. This implies, of course, /^(A) = cx) for 
every interval containing E. Corresponding numerical results are presented at the end of this 
section. 

We now come to the question how the localization length diverges. If one considers a disorder 
potential varying on a length scale much bigger than the magnetic length, the motion can semi- 
classically (in the magnetic field) be well approximated by a motion along equipotentials of the 
disorder potential. Percolation theory now indicates that there is one critical energy E^ to 
which correspond extended states. Moreover, at this energy the localization length diverges as 

(xlE-E^l-". (69) 



X{E) 



Trugman evaluated i^ = ^ [iy], but more thorough examination of the effect of tunneling 



between the wells and mountains led Mil'nikov and Sokolov to predict ^ = | [^^l- ^^ spite of 
the simplicity of the theoretical approach, the agreement with numerical results obtained by 
three different methods is very good. These latter are the Thouless number method, the Chern 
number method and the third one is a scaling theory approach; the original ideas behind are 
due to [0, IQ and |^3|, |^ respectively. 



The Thouless number method is the most elementary because only energies and no eigen- 
vectors need to be calculated. In the us interesting case of strong magnetic field, calculations 
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were made by Ando |Q. The Chern number method IQ is described in some more detail in 



the next section. Let us concentrate on the scahng theory approach. One considers a cyhnder 
of (small) circumference M and extended to infinity in one direction (denoted by the variable 
z). The magnetic field is perpendicular to the cyhnder surface. Chosing the Landau gauge 
A = (0, 0, Bz), the eigenstates of the free Hamiltonian are centered on circles around the cylin- 
der corresponding to fixed values of z. This discrete set of eigenstates will serve as base for the 
Hilbert space; in this base one also expresses Hamiltonian perturbed by disorder. Now for fixed 
circumference and fixed energy, one calculates numerically the disorder averaged Lyapunov ex- 
ponent Lm of a wave function by the transfer matrix method. Finite size scaling theory ^3 



then allows to estimate the behavior of the localization length near the critical energy; in par- 
ticular, the exponent z/ in equation (|69D can be calculated. In presence of magnetic field, this 



method was first used in 45, 431; more recent and more extended numerical studies 46, 95 



confirm the theoretical prediction u = ^. It seems that the exponent is independent of the 



3 

disorder distribution and the Landau level index 



What is the behavior of the wave function at the critical energy E^ itself? Pichard and 



Sarma |43| suggest the following in order to calculate the exponent x for the powerlaw decay 



of the wave function ipE^'- 



\'^Ec{r)'ipEc{r')\ oc \r-r 



One calculates as before the Lyapunov exponent Lm for different M at fixed energy E^. 
With the hypothesis of conformal invariance, x can then be shown to be equal to the slope of 



the linear relation between M and 1/Lm- Chalker and Coddington [Q estimate x ~ 0.27. We 
remark that this implies that the weak localization condition P<Ec € iS is violated. Therefore, 
at critical energies (that is at the centers of the broadened Landau levels) the Chern number 
may change. 



6.3 Chern numbers and localization in Harper's equation 

The Chern characters or Chern numbers of the energy bands of Harper's equation corresponding 
to rational flux have already been computed by Thouless, Kohmoto, Nightingale and den Nijs 
Q; we gave some insight in Section p.6| . Ando ^^ ^ and more recently Tan |]98| have made 



numerical studies on the influence of an added disorder potential on these Chern numbers. We 
describe some of the results because they constitute a good illustration of theoretical ideas. 

Let us first recall that Harper's equation (|Tl|) is obtained in two different limiting cases 
p9| , p8|] . One consists in adding a weak periodic potential Vp to the Landau Hamiltonian Hi, 
for example, a two-dimensional sinusoide; projecting Hi + V^ on a subspace corresponding to 
one Landau level then leads to Harper's equation. On the other hand, if one adds a strong 
periodic potential (which means that its well is much deeper than huc), the low energy region of 
the spectrum is well described by a tight binding model (equation |l^). As shown in Section p.6|, 
Bloch's representation also leads to Harper's equation. 

The spectrum of Harper's equation is well known and its dependence on the magnetic flux 
per unit area is represented in Hofstadter's butterfly [ |10CI| ]. The parameter a appearing in 
equation (0) is proportional to (f) in the tight-binding case and proportional to the inverse of 
(f) in the case of weak periodic potential. Let us now suppose the rationality of a = -, where p 
and q are relatively prime integers. The spectrum then consists of p energy bands of extended 

th 

states. The Chern number of the r band can be calculated with the Diophantine equation 
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in the following way P, |101|| . There exist unique integers s^, tr (except for even p and r = p/2) 
such that \sr\ < | and satisfying: 

Srq + trP = r r = 0, . . . ,p . (70) 

The Chern number of the r band and therefore its contribution in units of ^ to the Hall 
conductance is given by t^, — t^-i in the weak periodic potential case and by Sj. — Sr-i for the 
tight-binding model. Remark that as to = and tp = 1, the sum of all p Chern numbers in the 
weak periodic potential limit is equal to one. This reflects the fact that the bands result from 
the splitting of one Landau level with Chern number 1. On the other hand, sq = Sp = 0, such 
that the sum of the Chern numbers in the tight binding approximation is zero. 

Now a disorder potential will be added. This could be, for example, an Anderson-type 
perturbation (compare Section |5]^), but in the numerical simulations a densely distributed 
potential with two values V^, —Vd was used. The fractal structure of Hofstadter's butterfly 
disappears as soon as disorder is turned on. Because many gaps close, it becomes possible to 
label the remaining energy bands by a finite number of rationals for any (maybe non-rational) 
parameter a. What are the Chern numbers of the bands and how do they evolve as a function 
of the disorder strength? 

Before starting the discussion, let us comment on the physical relevancy of the two approx- 
imations. In the usual QHE samples presented in Section p.2| , neither seems to describe the 
reality: the disorder potential is much stronger than the weak periodic potential due to the 
ions within the surface. On the other hand, because the mobility is very high, the tight-binding 
approximation is not appropriate. However, lately so-called anti-dot superlatices seem to open 
the possibility of the experimental realization of parts of Hofstadter's spectrum (of course, only 
in a very approximative way). In these systems the ideas we present here could be tested. 

For all details on how the numerical results were obtained, we refer to the original works 
[]97| , |98[| , but let us describe the principal steps. Of course, the calculations are made in finite 
size samples and therefore disorder is just treated in an increased unit cell, causing thus a 
splitting of the Harper bands corresponding to a fixed rational flux into subbands. Chern 
numbers for these subbands can be calculated by the contour integrals given in equation (|T^). 
Remark that because we are in dimension two, these subbands have corresponding Hilbert 
sub-spaces generically well separated except for in isolated points. If this happens, the Chern 
number is simply calculated for a two dimensional fiber bundle (this is also the case for the 
closed central gaps of the Harper equation with even denominator). In this way, one gets a 
distribution of Chern numbers over the subbands of every Harper band. Now, changing the 
disorder configuration will give a different distribution of the Chern numbers over the subbands. 
This reflects the fact that changing the disorder configuration is equivalent to chosing another 
probability measure on the space of disorder configurations; we already commented that this 
will change the Chern numbers (see Section [4.7| ). However, the sum of all Chern numbers of 



the subbands is constant and equal to the Chern number of the Harper band, as long as the 
disorder is weak. As it becomes stronger, the Harper bands overlap energetically; the Hofstadter 
spectrum loses more of its structure as indicated above. But the distribution of Chern numbers 
may still be calculated. 

Now, for every strength of disorder, averages of the Chern numbers over the disorder con- 
figurations can be calculated for every subband. The numerical results then indicate that, on 
average, the Hall current is essentially carried by states corresponding to one central band sit- 
uated somewhere at the center of the Harper band. This belief is supported by the calculation 
of the localization length by the Thouless number method |^ : it diverges at about the energy 



67 



Disorder 








-3 2 



2 -3 2 



Energy 
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corresponding to the Hall current carrying subband |^ . In the limit of an infinite sample, one 



therefore expects the Hall current to be carried by one critical energy, just as in the case with- 
out periodic potential. This is, in fact, not very surprising because the semiclassical argument 
presented in Section |6.2| applies to the present situation as well. 



In summary, the numerical calculations suggest the phase diagram |98| shown in an idealized 



form in Figures |6.3| and |6.3| . Starting at rational fiux without disorder, the evolution of the 
current carrying states and their Chern numbers is given as a function of increasing strength 
of disorder. At some critical values, they amalgamate; the Chern number of the resulting 
state is then given by the sum of the merging ones (this sum may also be zero. Hall currents 
then annihilate each other). To conclude, let us comment on the limiting behavior for strong 
disorder. In the weak periodic potential approximation, one eventually obtains a Chern number 
one as for the Landau level. It is an interesting question what happens to the Chern number 
as the quantum Hall regime breaks down, certainly the corrections to the Kubo formula will 
be important at very high disorder. In the tight binding approximation, the system evolves to 



zero Hall current, this refiects the results of Section lO 



7 Introduction to the FQHE 

7.1 Overview 

While in 1982 Thouless, Kohmoto, Nightingale and den Nijs made an important step towards 
the understanding of the IQHE |^, D.C. Tsui, H.L. Stormer and A.C. Gossard from Bell 
Laboratories made the surprising discovery that plateaux of the Hall conductivity could be 

2 

observed not only at integer, but also at fractional multiples of ^ ||1U2|| . Using a high-quality 
GaAs — AlxGai_xAs heterojunction, they observed a plateaux at cfh/^ = \- The corresponding 
minimum of the direct resistance was also observed. The filling factor was varied by changing 
the magnetic field while keeping the charge-carrier density fixed. In order to obtain the filling 
factor 1/ = |, they needed a magnetic field strength of about 15 T. The experiment was 
undertaken at different values of temperature (0.48, 1.00, 1.65 and 4.15 K) and even at the 
lowest of them, the accuracy of an on the plateau was much more modest than in the integer 
effect. 

Since this first experiment, the FQHE has been observed for many values of a^j- The most 
striking results are the following: 

1. Let (Th/\ = - with p, q relatively prime integers. Then, q is odd ('odd-denominator rule'). 
Recently, some deviations of this rule have been observed, but this seems to correspond 
to non spin-polarized electrons. In the sequence, we will not consider this case. 

2. The observation of the fractional effect requires very clean samples. It is easily destroyed 
by impurities causing inelastic scattering. Moreover, the stability of the plateaux with 
respect to temperature and impurity effects is very dissimilar for different fractions. In 
general, fractions with small denominator seem to be more stable. Thus, the plateaux 
appear in a well-defined order when lowering the temperature. This gives rise to series of 
fractions, also called 'hierarchy of states'. 



69 



3. In all devices, the accuracy of the Hall conductivity on the plateaux is much less for 
fractional than for integer values. Whereas in the IQHE we have -^^ ~ 10"^, this is 
about 10"^ for the FQHE. 

Although we have given a rigorous explanation of the IQHE, the situation is much less sat- 
isfactory for the FQHE. The QHE contains two distinct aspects: the existence of plateaux 
(which, for the IQHE, follows from localization of the states at the Fermi edge) and the inte- 
grality (respectively fractionality) of the Hall conductivity on them (which has a topological 
origin). Usually, these two points are not well separated in the literature on the FQHE. In 
fact, there are conceptual difficulties in understanding the role of localization in the case of 
interacting electrons. The mechanism which gives rise to the existence of plateaux is quite un- 
clear. The main problem is that no Kubo formula is known to calculate the Hall conductivity 
for interacting Fermions. Instead, one usually follows another strategy: since we know that 
non-interacting electrons only exhibit an IQHE, it is clear that the FQHE might appear when 
electron-electron interaction can no longer be neglected. It is generally believed that at certain 
values of the filling factor, the ground state becomes incompressible owing to the repulsive 
two-particle interaction. Thus one tries to find explicit ground states which exhibit a fractional 
Hall conductivity. Then arguments are given to explain that this value remains constant while 
the filling factor is slightly changed. The most promising states that have been found show 
a remarkable property: they exhibit particle-like excitations with some charge q ^ —e. The 
states are only defined for specific values of the filling factor, but it seems that the existence of 
quasi-particles is stable against perturbations if - is a fractional number. Apart from numerical 
evidence, there is no serious explanation of this property. The stability of quasi-particles is due 
to an energy gap separating the ground state from the excitations. It is conceptually not clear 
why the existence of such a gap is related to the fractionality of the quasi-particle charge. No 
rigorous results are known here. 

7.2 Lauehlin's ansatz for the v = — groundstate 

Up to now, most attempts to explain the FQHE - in particular the classical hierarchy of 
Haldane and Halperin [|103|, |104|| as well as Jain's composite-fermions approach ||105|| - are based 



on Laughlin's wave functions which describe a finite number of electrons in a finite volume, 
yielding a filling factor z/ = — with m an odd integer. In the following, we will review the 
ideas which led Laughlin to his proposal. Following Haldane, we will then construct a class of 
Hamiltonians for which Laughlin's states are exact ground states. 

Before doing this, we want to clarify the role of these states and try to justify our interest in 
them. As we have already mentioned, there does not exist a satisfying theory of the FQHE at 
the moment. In particular, the universality of the effect and its insensitivity to the exact form 
of the Hamiltonian is not quite clear. The main advantage of our approach to the IQHE is that 
we do not need to know the exact form of the one-particle potential, but only have to impose 
some mild conditions on it. To explain the existence of plateaux and to prove the integrality 

2 

of (Th/^, we did not need to calculate the ground state explicitly. In our opinion, despite the 
wide use of Laughlin's wave functions in the literature, the explicit form of the ground state 
does not lie at the heart of the theory, neither for the integer nor for the fractional effect. 

Nevertheless, in our theory of the IQHE we were not able to give the explicit values of 
the Hall conductivity just by means of general considerations. In fact, the value given by the 
Kubo-Chern formula might depend on the choice of the probability measure P on the hull of 
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the one-particle Hamiltonian. To derive the 'right' values, we computed explicitly the Chern 
character of the eigenprojections of the Landau Hamiltonian. A homotopy argument then 
shows that these values remain constant when the disorder is switched on, at least for small 
values of the potential. 

The states which arise from Laughlin's wave functions in the thermodynamical limit might 
play the same role for interacting particles: They yield the explicit values of the conductivity 
for a Hamiltonian with a restricted two-particle interaction and without one-particle potential. 
The main task which remains is to show that these values stay constant as the Hamiltonian is 
deformed continuously. This is going beyond the scope of our article. 

We will now come to the to the description of Laughlin's states. The general Hamiltonian 
for N electrons is given by 

1 ^ e - ^ 

^m^ j=l C l<j<k<N j=l 

where in our case V x A{x) = Bcz, B = const. Thus we assume that the two-particle interaction 
depends only on the distance between the particles. Under certain conditions on the functions 
U and V, H^ is essentially self-adjoint on C^(R^^). Furthermore, the A^-Fermion space 
A^-^^(R^) is an invariant subspace. Choosing suitable boundary conditions (we will choose 
Dirichlet boundary conditions on the disk Aji = {x G R^ : |a^| < R}) yields a self-adjoint 
operator H^^^ in the A^-fold tensorproduct 0^ -^^(A). We will denote its restriction to A^ -^^(^) 
by the same symbol. Since the Laughlin states should play the same role for interacting electrons 
as the ground-state of the Landau-Hamiltonian does for non-interacting ones, we will set V = 0. 
In fact, up to now it is unclear what happens when the one-particle interaction is switched on, 
even for small values of the potential. The basic idea which leads to Laughlin's states is the 
following: for a strong magnetic field, the two-particle interaction can be treated perturbatively; 
if we start from a system of non-interacting electrons in the lowest Landau level, the excitations 
to higher Landau levels in the perturbation series can be neglected. Thus, for filling factors 
z/ < 1, we assume the groundstate to be a linear combination of Slater determinants of one- 
particle wavefunctions belonging to the lowest Landau level. 

For V = 0, the remaining operator Hn is rotationally invariant. Thus, we can look for 
joint eigenvectors of energy and total angular momentum. It is clear that, for a repulsive 
two-particle interaction, the particles will escape to infinity when not restricted to a finite 
volume. Classically, the trajectory of an electron in a constant magnetic field with angular 
momentum / and energy E is a, circle around the origin with radius R = J-^^l^- Thus, for 
E = ^, restricting the particles to a disk with radius R is essentially equivalent to restricting 
the angular momentum to values less or equal than "'"^^'^'^ = ^-4 = ^ = ^^^.h (with £g the 
magnetic length). Thus, instead of calculating the ground state of H^/^, we are looking for an 
eigenstate of H]\f built of lowest Landau level wavefunctions with maximal angular momentum 
A^$. Laughlin, strongly influenced by the theory of liquid ^He, proposed the ansatz 



tpizi, . . . , zn) = n dizk 



z 



]l 



i<j<k<N 

for the A^-particle groundstate. Here and in the following, we set the magnetic length to unit 
and we use complex coordinates z = x + ty in the plane. The lowest Landau level is spanned 

1 I I 2 

by the functions 0m oc ^"-e'zi^i , thus 

ij{z,,...,ZM)= n /(^fc-z,)e-^^"-l^'l' 

^<j<k<N 
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with an analytic function /. Since ip should be anti-symmetric under particle exchange, we have 
f{—z) = —f{z). Furthermore, as we are looking for eigenstates of total angular momentum 
M and the operator of angular momentum for one particle is given by zd^ — zdz, the product 
in the above formula has to be a polynomial of degree M in each variable Zj. From this, 
it follows that f{z) = z^ with an odd integer m. We will denote the resulting normalized 
wave function by 4'N,m- The total angular momentum is M = ^ ^~ ' m and the particles are 
essentially restricted to a disk with radius R ~ J2m(N — 1) (in units of the magnetic length) 
or, equivalently, z/ ^ — . 

Before we investigate the question for which Hamiltonians Laughlin's ansatz yields an exact 
groundstate, we want to make some remarks on the thermodynamical limit. The QHE (the 
integer as well as the fractional) is assumed to be exact in the case of infinite volume and zero 
temperature, thus necessarily we have to involve statistical mechanics. Laughlin's states are 
states for a finite number of particles in a finite volume, i.e. we are working in the canonical 
ensemble. The Gibbs states are 

tr(e"'^^'^'Aa) 



tr(e"'3^iv.A) ' 

where a is a bounded linear operator in /\'^ L'^{A). At zero temperature, they reduce to the 
vector states 

(a)oo,iV,A = {'ipN,A,ailJN,A), 

where 4'n,a is the (non-degenerate) groundstate of H]\r/^. To investigate the limit N ^ oo, 
A -^ H? with fixed particle density i^, we have to include vector states with an arbitrary 
number of particles, i.e. we have to work in Fock space. The algebra of observables should be 
a suitable subalgebra of bounded linear operators in this space. We will choose the canonical 
anti-commutation relations (CAR) algebra. We recall the definition: Let 

oo 
N=0 

denote the Fermi-Fock space. For (f> G L^(R^) and ^ E 0^L^(R^), we shall now define 
a*{(j))tlj := VN+l(j) ® ^p. Let P^ denote the projection onto A^ L^{R^) in ^^ L^{R^), i.e. 

where the sum runs over all permutations of (1,...,A^). Let al(0) = P^'^^a*{(j))P:^ on 
A^'L^R^), hence, for ^ e A'^^'(R') 

1 ^+^ 

a*_{(f))-^{Zi,...,ZN+l) = r^^—r I] (-)-'^V(%-)V'(^l, • • • , 4'' • • • ' ^^+l)' 

where Zj denotes that the jth variable has to be omitted. This yields a bounded linear operator 
a*_{(f)) in JF_(L^(R^)) with ||al(0)|| = ||0||. The adjoint operator is given by 



a-{(j))^lj{zi,...,ZN-i) = VN / d z(j){z)tlj{z,zi,...,ZN-i) 

JR2 
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for ip G A^-^^(R.^)- These operators satisfy the CAR 

{a_(0),a„(V^)} = {a*_(0),a*_(^)} = 0, {a_(0), a*_(V^)} = (0,^). 

The CAR algebra over L^(R^), which we will denote by A, is the C*-algebra of bounded linear 
operators in JF_(L^(R^)) generated by the identity and the set of operators {a_(0) : cf) G 
L2(R2)}. 

We regard a state as a linear mapping (■) from this algebra into the complex numbers such 
that {A* A) > VA G .4 and (1) = 1. Given a sequence {tpN)Ne-N of vectors ipN e A^-^^(R-^), 
we can ask for the limit of the corresponding vector states. Here, the limit has to be understood 
in the weak-*-topology, i.e. (■) = limiv-+oo(-)Ar if and only if {A) = limAr_^oo(^)Af VA G A. It 
is clear that the limit is not a vector state. 

As an example, we will treat the case m = 1. Then, Laughlin's wave function is a single 
Slater determinant (namely Vandermonde's determinant). The normalized wave functions are 

-L T — r 1 V^^ I |2 

\/7l^V2 ^ A\...N\i<J<k<N 

We will denote the limit state by (■)oo,i- It is completely determined once the correlation 
functions {IYjLi^*-i'Pj)Ilk=i^-ii'k))oo,i are known. Obviously, these vanish for m ^ n. We 
will calculate the two-point function exphcitly: We have (al(0)a_ ('?/'))oo,i = ijp, Tcj)), where the 
integral kernel of T is given by 

{z', Tz) = lim A^ / d^Z2 . . . d'^ZN ipN,i{z, Z2, . . . , zn)iPn,i{z', Z2, . . . , zn) = 

1 N-1 (zz'\j 1 



Af^oo 27r j^Q j\ 2n 



e 41 



Thus, the operator T is just the projection onto the lowest Landau level. As a consequence, we 
get the correct filling factor u = 1. Furthermore, it can be shown that (■)oo,i is a quasi- 
free state, i.e. the 2?7,-point functions are certain sums over products of two-point func- 
tions, for example (a!l(0i)al(02)a-(^i)Q -(V'2 ))oo,i = (a!l(0i)a-(^i))cxD,i(a!l(02)a-(^2))oo,i - 
(al(0i)a_('?/'2))oo,i(ol(02)fl-(V'i))oo,i; see ||106|| for an exact definition. But this is just the limit 



Gibbs state (for temperature T = and filling factor z/ = 1) for non-interacting electrons in a 
constant magnetic field. Thus in the thermodynamical limit we get the exact groundstate. 
This can also be seen from the work of Haldane ||107|| which we will now review. Let 11 denote 



the projection onto the lowest Landau level in the one-particle space. We will investigate the 
Hamiltonian H^^eff '■= H^^HnU'^^ . The kinetic part just equals ^^j^. In order to calculate 
the integral kernel of 11 (8> IlU{\zj — Zk\)Il ® 11, we perform the unitary transformation 

U : L\R^) (g) L2(R2) ^ L2(r2) (g) ^^(R^), 

(U4'){Zjk, Zjk) = ip{zj = Zjk + Zjk, Zk = Zjk — Zjk). 
Then, 

u{ii ® ii)u{\z, - Zk\){n ® Ti)u-^ = n ® {Yiu{\zjk\)^). 

Since U{\zjk\) is rotationally invariant, the second factor in the tensorproduct is diagonalized 
by the simultaneous eigenstates of the Landau Hamiltonian and angular momentum. In this 
basis, 

Wm^U {\Zjk\)(j)mi) = UmSm,m' 
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with 

" Jo 22™+im! ^4^ ^ ^' 

where we have reintroduced the magnetic length £g. Thus, in some sense, the projection onto 
the lowest Landau level yields a 'quantization of inter-particle separation'. For example, for 
the Coulomb interaction f/(r) = ^ we have 

e^ r(m + I) 

"^ ~ 47r 2m! ' 

For m ^ oo, this falls off as -^. It is remarkable that the coefficients Um are independent of 
ij3 exactly for the Coulomb interaction. 

As a result, we state: The set of functions 

^{Z, z) oc Z^^^™e-^(l^l'+l^l') oc {z, + zu)''{z, - ^,)™e-3(l^.l^+l^^n 

with M,m G {0, 1,2, . . .} forms a basis of eigenfunctions of 11 IlU{\zj — 2:^1)11 IT. Each 
eigenvalue Um is infinitely degenerate (according to an arbitrary angular momentum of the 
center-of mass motion). 

The general two-particle wave function in the lowest Landau level with angular momentum 
maximal equal to A^$ is given by 

^(zi,Z2) = P^Jzi,Z2)e-i(l^^l'-*-l^^l') 

where P/v^ is a polynomial of degree N<^ in each of its variables. These states form a ( 2* 
dimensional subspace, but this subspace is not invariant under the two-particle interaction. Let 
us assume that the sequence Um decreases monotonically. Then the desired ground state would 
be a state with relative angular momentum A^$, or, since we are dealing with Fermions and 
therefore the relative angular momentum has to be odd, N^ — 1. In the latter case, we have 
a degenerate ground state since we can choose M = or M = 1. Let us remark that these 
functions are not equal to Laughlin's ansatz for N = 2. 

In the case of more than two particles, it is in general not possible to give the ground- 
state. The general A^-particle wave function in the lowest Landau level with maximal angular 
momentum equal to A^$ is given by 

1 v^A'^ I |2 

ij{zi,...,ZN) = Pn^{zi, . . . , ZN)e~'i ^i-i™ , 

where again PjVi, is a polynomial of degree Nq> in each of its variables. These states form a 
f *j^ j -dimensional subspace. Let us require that only pairs with a relative angular momentum 
of at least mo appear in the wave function. This is possible if and only if mo < j^ -^ i, and 
the general form of the wave function is given by 

_ — _ 1 v-^N" I |2 

ij{zi, ...,zn)= PN^-mo(N~i) n (^fc - ^,)™"e-^ ^'=^ '"' • 

^<j<k<N 

Let us now consider the truncated Hamiltonian which results from setting Um = for m > mo. 
Then, the above ip is an eigenstate with energy 0. Furthermore, if we assume Um = for 
m, < mo, it is a ground state. In general, for u < — , the ground state is degenerate due 
to the polynomial pre-factor PAr,j,_mo(Af-i)- Let us assume that m,o is odd. Then, to get a 
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Ferinionic state, this pre-factor has to be symmetric under particle exchange. Thus, if we 

assume N^ — mo{N — 1) = which is equivalent to :^ = u — -^ > z/, we get a unique ground 

state which is equal to Laughlin's state 4'N,m with m = rriQ. For ttiq even, PN^-mo(N-i) has to 
be antisymmetric. There is (up to a constant factor) an unique possibility if the polynomial 
is of degree one, namely Pi{zi, . . . , zn) = I\i<j<k<N{^k — Zj). In this case, we get ipN,m as a 
unique ground state for u = — and m = itiq + 1. 

Let us recapitulate the result: if we project the Hamiltonian onto the lowest Landau level 
and assume that all coefficients Um vanish for m greater or equal to some rriQ and Um > 
otherwise, then the unique (infinite volume) groundstate with filling factor u = — for mo odd 
respectively u = ^'', ^^ for rriQ even is the limit state which arises from Laughlin's wave functions 
ipN,mo or ^Ar,mo+i respectively. 

In fact, we have not proved this rigorously. First, it remains to show that Laughlin's wave 
functions indeed yield the correct filling factor in the thermodynamic limit. This has been 
shown by Laughlin by means of a somewhat heuristic argument [[108|| . Then, it has to be shown 
that the restriction of angular momentum indeed yields the same thermodynamic limit as a 
classical boundary condition. This should not be too hard to prove. 

7.3 The elementary theory of the u = ^-FQHE 

The elementary theory of the u = — -FQHE is solely based on Laughlin's states. Thus we 
assume in the following that the correct ground state is indeed given by their thermodynamical 
limit. We then have to explain that 



e2 



1. the Hall conductivity for the state given by ipN,m equals ^^, 

2. plateaux of the Hall conductivity arise when the one-particle potential is switched on. 

We remark that the first point is just the classical result cth/^ = v- Indeed, for l^ = 0, the 
classical result should hold for any value of the filling factor. 

To calculate a^, Laughlin |p.08|| proposed a Gedanken-experiment similar to that we have 



already described in Section |2~5|. We consider the iV-particle wavefunction ipN,m- We introduce 
a infinitely thin solenoid at the origin, perpendicular to the disk in which the particles are 
confined. Then we force adiabatically one fiux quantum $o = — through the fiux tube. We 
assume that il'N,m is the non-degenerate groundstate for some Hamiltonian and that there is 
a gap between the groundstate energy and the rest of the spectrum. Let us mention that 
this 'gap-condition' is a sufficient condition for our considerations, but in general need not be 
necessary for cfhI\ to be equal to ^. After the fiux quantum has been added, the Hamiltonian 
has not changed up to a gauge transformation. Since we can apply the adiabatic theorem, we 
know that the system is still in an eigenstate of the Hamiltonian. The single-particle wave 
functions ^"^e"*'^' evolve, up to a phase factor and normalization, to 2;™'+^e~4l^l during this 
operation. The particles which were essentially confined in a disk with radius R = a/2A^$ at 



the beginning are now restricted to the radius R' = J2{N^ + !)• These two disks differ by an 
area 2it, and since the charge density equals uniformly 2^, a charge ^ has left the original 
disk while adding the flux quantum. Similar to the IQHE, this 'gauge argument' does not prove 
directly the fractionality of the Hall conductivity, but shows that the classical result cth/\ = v 

2 

is valid. Fractionality of (Th/j[ follows from the fractionality of the fllling factor. 
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Let us discuss the excited state which appears in Laughhn's Gedanken-experiment in the 
case of Haldane's truncated Hamiltonian. After the flux quantum has been added, the single- 
particle angular momentum is restricted to values less than or equal to A^$ + 1. Since we are 
looking for an excited state, the relative angular momentum for each pair is still at least m. 
Thus, 

1p{zi, ...,Zn)= Pl{zi, ..., ZN)lpN,m{zi, ..., Zn)- 

In general, Pi(2;i, . . . , zjq) oc Y{^=i{zj — Uj) with arbitrary Uj. But as Pi has to be symmetric 
under particle exchange, all Uj have to be equal. This yields the excited states already proposed 
by Laughlin p8| : 

N 
Szo'^N,m{Zi, ...,Zn) = Y[{Zj - Zo)i:N,m{Zi, ...,Zn)- 

i=i 
The adjoint operator is given by 

j=l l<k<l<N 

The operators Szo and S*^ are called 'quasiparticle' respectively 'quasihole' creation operators 

1 I 1 2 

for the following reasons: Denote by ^N,m,z the A^-particle state e"*;^'^' SzipN,m- Then, one can 
show that 

Thus ^z,m looks like the wavefunction of a 'quasi-particle' with charge — in the lowest Landau 
level. More formally: Let ip G L2(R^) (a quasi-particle wave function). Define the A^-electron 
state 



,|2 



QN,mi' ■= /- j <frj%l){ri)e 4™l''l^5'^^;v,m- 

Then, for ^'^(r/) = J^ — e~i^^^^^ "'"'''I -^z?^)^ ^g have QN,m'ipz = ^N,m,z and moreover {■tpz^tpz') = 

{QN,m4'z, QN,m'ipz')- 

Of course, the state $2 ^ is in general not an exact eigenstate of the truncated A^-particle 
Hamiltonian. But as A^ — >■ 00, its overlap with an eigenstate becomes larger and larger. Let us 
assume, for simplicity, that it is an eigenstate. Since the quasi-particle state ipz is an eigenstate 
of the Hamiltonian Hq = (-V,, + —A{ri)Y, we have 

{^N,m,z,HN,eff^N,m,z) OC {^Jz, Hg^pz) , 

where, for simphcity, we have denoted also the truncated Hamiltonian by H^^eff- One is led 
to define a quasi-particle potential Vg by the equation 

N 1 /• 

{^N,m,z,J2yi^:i)^N,m,z) OC (t/-,, l^,^,) = ^^ J A e"^'^"'^! V,(r/). 

Thus, we have defined an effective Hamiltonian which governs the quasi-particle motion. This 
yields a qualitative explanation for the existence of plateaux: If the quasi-particles which are 
created by an electric field are localized due to randomness of the potential Vg, they do not 
contribute to the Hall conductivity. The FQHE looks like the IQHE for quasi-particles. 
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But this only concerns the mechanism which causes the existence of plateaux. The frac- 
tionality of (Th/j^ is due only to the fact that the charge density of Laughlin's states is exactly 
■2^- As we have shown, these states are indeed the non-degenerate ground states for some 
Hamiltonian, but it is not clear at all why the quasi-particle charge should be invariant under 
slight deformations of the ground state. Thus, the topological origin of the FQHE remains 
undiscovered. 



7.4 The role of gauge invariance and incompressibility 

We have tried to emphasize that the main problem of the common approach to the FQHE is 
the lack of a Kubo formula to calculate the Hall conductivity. To calculate it for Laughlin's 
states, we used gauge-invariance and assumed the existence of a gap separating the ground 
state energy from the rest of the spectrum. These two assumptions - gauge-invariance and 
incompressibility - also form requirements for a completely different approach which is based 
on a work of Halperin |^ and has been developed further e.g. by Wen ||109|| and Frohlich et 



al. ||110| , |111| , |112| , |113| , |114| , |115|| . We will particularly refer to the work of the latter group. 



Let us first look at incompressibility from a classical point of view. How does the system 
respond to external "small" fields £ and B (where the latter does not include the strong magnetic 
field which causes the quantum Hall state)? In (2 + l)-dimensional space-time (we will only 
consider cases where space is a subset A C R^ - in particular a disk around the origin - with 
the usual 'fiat' metric and space-time will then be R x A), we only have two components 
of the electric and one component of the magnetic field. In covariant notation, we have the 
electromagnetic- field tensor F = ^F^^dx^ A dx'^ with x° = ct and 

(U Ox Oy 

—Ex —Bz 
-£y B, 

The homogeneous Maxwell equations dF = reduce to Faraday's induction law 

c 

Let j denote the electric current density. Then we have the Ohm-Hall law (respectively the 
definition of the conductivity tensor a) j = aS. Furthermore, a = p~^, where p is the symmetric 
resistivity tensor. Incompressibility means that the diagonal components of p vanish and hence 

a = 

Furthermore, we have the continuity equation (i.e. charge conservation) 

-%° + Vj=0, 
c 

where cj° denotes charge density minus the uniform charge density of the unperturbed quantum 
Hall state. From these three basic equations, we get dtj'^ = audtB. Since, with our above 
definition of j*^, we have j° = for i3 = 0, the integration of this equation yields j° = (ThB. 
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We introduce the 2-form J dual to {j^,j): 

J = -^Jfiudx A ax , [Jfiu)fj,,ue{o,i,2} = 

Then, the continuity equation reads dJ = and we can combine the Ohm-Hall law and the 
connection between j^ and B^ into one equation: J = —afjF. Since we work in simply connected 
space-time, we have 

J = db, F = da 

with 1-forms a and b. Hence, db = —anda. This is the Euler-Lagrange equation for some action. 
Whether we regard it as an equation for a or for b is just a matter of taste, corresponding to 
the choice that we can either create a current by applying an external field or, equivalently, 
can create a Hall voltage by forcing some current. For our purpose, we will choose the latter 
possibility. Thus, our 'dynamical variable' is the gauge potential whereas we keep the potential 
for the current density fixed. A possible choice for the action is: 



^A,6(a) = -Y 



/ a A da + b A da + b.t., 

JAxR JAxR 



where, since we work in a finite volume, we have to include some (unknown) boundary term 

'b.t:. 

Can these phenomenological equations be derived from a microscopical description? And is 
there a connection between the quantization of the classical theory (by means of path integrals) 
and our original quantum mechanical problem? Formally, the Hamiltonian H\ = 0*^=0 Hn,a 
(with H]\fj>^ as in the previous sections) can be described by means of 'fields' \E' and \E'* which 
satisfy at each time t the CAR 

{^{x,t),^{y,t)} = {^*{x,t),^*{y,t)} = 0, {^{x,t),^*{y,t)} = 6{x-y). 

The fields act on the Fermi- Fock space JF_(L2(R^)). The particle-number operator is given by 
-^A = /a ^*(a^, t)\E'(x, t). We will work in the grand-canonical ensemble. We have 

HA{t)-fiNA = -Y- J d^x[V ~—dtot{x,t)]^*{x,t)[\/ + —dtot{x,t)]^{x,t) + 
+ f d^x'^*{x,t)[V{x)- fi + ea°]'${x,t) + 

J A 

+ 11 d^xd'y^*{x,t)^*{y,t)U{\x-y\)^{y,t)^{x,t). 

Here, atot represents the potential for the total magnetic field (i.e. the potential A for the 
constant field Bcz plus some - in general time-dependent - source term a) and a° also is a 
source term. 

We then have the Heisenberg equation of motion 

dt^{x,t) = ^[H^{t) - fiNj,,^{x,t)]. 

But this is the Euler-Lagrange equation derived from the action 

5^(^*,^) = f dtU d'^X^*{x,t)thdt^{x,t) - HA{t) + fMNA 
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The partition function (at zero temperature and chemical potential /i) is given by the path 
integral 

We define the effective action as S'^ := Mn ^a- Regarded as a functional of the gauge potential 
{O'tot) ij.dx^ = a^dx^ — Stotdx, it is the generating functional for connected time-ordered Green 
functions of the current- density operators 



\1/*\1/, 



5|^{|(V-£a«)*1*-r(V + £a„)*}. 



We have 



{Tf^{x^)...r'{Xr,)r = ^^^^'^- 



6 



6 



-SI" {a', a) 



We now perform a scale transformation, which means the following: We enlarge the domain 
A by multiplying the space-time R x A with a scalar A > 1. We introduce new coordinates 
C, = (^°,0 = A~^(x°,x) (thus ^ G A) and define the rescaled sources a^(A^) := \~^a^{C,), such 
that a^(A^)c/(A^) = a^{^)d^. The 'scaling-limit' A — ;> oo corresponds to a thermodynamical, 
but also adiabatic, limit. We assume some 'strong-clustering property', namely we require that 
for n < 3 the distributions 



„n+l 



l^tJ.l...tJ.nrf 



• Q 



\imy-{Tf^{x^,)...j^-{x^^)y 



are local distributions, i.e. 



supp 



LfJ.l...fJ.n 



{(6 



n\ A^oo 



,en)e(RxAr:ei 



^n}. 



Then, assuming some differentiability condition on 5*^ , we have up to second order (with 
A = -Adx): 



Sli^{A + a') = Sl{f{A)+ I d'x 

+ / d^xd^y 



A(RxA) 5a^(x) 



a=A 



[A(RxA)]2 



6a^{x)6a^{y) 



Sl{^ia)a'^ix) + 
Sl^I{a)al{x)al{y). 



a=A 



Now, 



5a^ (x) 



:eff 



a=A 



Sli'ia) = - {nx)ri^^. We define m) ■= ^^x-^oo X' {nXOr\a=A and J, 



as the 2- form dual to jc- 
Furthermore, we have 



6' 



2 5af,{x)6a,y{y) 



S 



AA I 



a=A 



\ {j,{^)ny))\=A "^ 



r^{X-'x,x-'y). 



By the above clustering property, we can split 0'^'' in the manner 

where a is a constant and R^'^ consists of second or higher derivatives of the (5-function. 



It can be shown [113] that neither the part with R^^'^ in the second-order term nor any of the 
terms of higher order contribute to the scaling limit up to some boundary term. In the language 
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of field theory, they are irrelevant, whereas the first-order term is relevant and the remaining 
part of second order (the 'Chern-Simons term') is marginal. Hence, a simple substitution in 
the remaining integrals leads to 

Sl^^{A + a) - Sl^^{A) ^-=^/ J,Aa + a f a A da + b.t. 

JAxR jAxR 

This is exactly our phenomenological action; thus the clustering conditions indeed correspond 
to incompressibility. By comparison, we expect that a = ^. Let us check this explicitly. We 
reintroduce the non-rescaled variables: 

S'/j/iA + a) ^-=^ S'/JiA) - f d^x {3''{x)Y\,al{x) +a [ d^xe^'''Pal^u{x)d,a]{x). 

JA(AxR) JA(AxR) 

By calculating the functional derivative, we obtain 

(/W)lA+a^ = {3'{x)YU-2e^^''d,a>;{x). 

For example, let us take as source the potential for a constant electric field in ?/-direction: 
Qq^x^^x) = —Sx2, d^ = 0. Then {j^{x)y\^_^_^x = {j^{x)y\^ — 2a£. Thus indeed we have 
a = ^. 

The boundary term in the effective action is in general unknown, but we have the require- 
ment that the whole action is gauge- invariant, i.e. S*^ (a + dx) = S^/ {a) for arbitrary smooth 
functions x- From this, it follows that 



b.t. {a + dx) — b.t. (a) = -— / dx A a. 

2 JdAxn 



We introduce light-cone coordinates u± = -j={vt ± ^L) on dA x R, where v is some (a priori 
arbitrary) velocity, 6 is the angle of polar coordinates in the plane and L is the circumference of 
the disc A. Then the above functional equation for the boundary term has the general solution 

b.t.(a) = f d^u ia+a^ - 2a+^a_) + Wia) =: A(a) + W(a) 
JdAxn. □ 

where a|aA = a+du+ + a_du_, d± = -^, □ = 29+(9_ and W is some arbitrary gauge-invariant 
function. 

A is known as the generating functional for time-ordered connected Green functions of chiral 



?7(l)-currents. We just state the result and refer the reader to ||116| , |117|| : Assume that there are 
current-density operators in two-dimensional space-time whose time-ordered connected Green 
functions are given by the functional derivatives of A. These operators, then are sums of 
derivatives of N massless scalar fields 0i, . . . , (p^ for some A^ G {1,2,...}. The dynamics of 
these fields is determined by the action 



5,(01,..., 0„) = 11 d'ud+<f)'Kd^<l)- 

Z JdA 

- / d^u Wd+d) - (d-d) - K-^a)^a,] + - f d 

JdA 2 JdA 



ua^a^, 



where i^ is a positive-definite A^ x A^-matrix, 0* = (0i, . . . , 0„), a* = (a, . . . , a) and a^. 
(a±, . . . ,a±), K = J2fk=ii^~^)jk ^i-iid we have to impose the chirality constraint 

9_0 - K~^a- = 0. 
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Indeed it can be checked that the effective action of this theory, which is defined analogously to 




5*2 , is given by |A. Thus after performing the 'Abelian Bosonization', we see that the total 

2 

effective action is gauge-invariant if and only if (^h/\ = i^- 

Up to now, we have (apart from positive-definiteness) no requirements on the matrix K and 
thus cfh/^ could be any real number. What are the requirements for a rational value of the 
Hall conductivity? It turns out that we have to impose some further condition. We will sketch 
the idea and omit the details. For simplicity, we will restrict ourselves to the case n = 1. Thus, 
we have a single massless scalar field. We introduce the variables a± := ■^a± (which have the 
dimension of an inverse length) and K := ^K (which is dimensionless) and we normalize the 
field such that the action is given by 

SM = 7^ I d\ |:|(9+0)(9_0) - [a_9+0 - (9„0 

2tC J-RxdA [ I 

We can break down into its chiral components: 

0(m+, M_) = 0l(m+) + 0ij(M-)- 

'L' and 'i?' correspond to the sign of the charge of the (quasi-)particles described by 0. Let 
us treat the part (pL- We have the current-density operator Jl{u^) oc (9+0l(m+), where at the 
moment it is not quite clear how to define the constant of proportionality. We will choose it 
such that the corresponding charge operator is given by 

e / 1 

Ql = -^ f dO ((9+0L - -^d+). 

The form of Ql is clear for a = and the general case follows from gauge invariance. Of course, 
this determines Ql only up to some constant factor. We will come back to this point which 
seems to be crucial to us. We define 'vertex operators' 

Vn{u+) =: e*"'^^("+) :, 

where ': :' denotes normal ordering. These operators obey the important relations 

TIC 

[QL,Vniu+)]=-^Vn{u+). 

K 
Thus, Vn{u+) creates a quasiparticle of 'charge' Ql = — -^^ If we create two such quasiparticles 

and interchange them, we get a phase factor e k , i.e. these quasiparticles obey fractional 
statistics in general, li K = n, then we get a phase factor e^*'^^. We now come to the 
announced additional assumption: Assume that the quasiparticles with 'charge' — e created by 
vertex operators are electrons and therefore obey Fermi statistics. Then we get the celebrated 
result 

i^ = 2/ + l^(7H = -i-^, leZ. 

21 + 1 h 

Let us just mention that in the general case, i.e. with more than one field, cth/^ is still a 
rational number, namely 

N 2 

aH= J2 (^"')ifcT' ^Jk e Z, Kjj G 2Z + L 
j,k=i '^ 
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Furthermore, if we also take quasiparticles with positive charge into account, cfh/\ is the 
difference between two such rational numbers. 

We end this section with two remarks: First, it seems to us that the last assumption is more 
or less ad hoc. Why should we exclude quasiparticles with fractional statistics and charge — e? 
Even if we assume that the electron is among the possible excitations, it is not clear that we 
can create this excitation by just applying one vertex operator. 

Second, there is still the question of normalization of the charge operator. Assume that we 
change its normalization by multiplying it with a factor 7. Then, the vortex operator Vn{u+) 
creates a quasiparticle with 'charge' — ^. Interchanging two such quasiparticles with 'charge' 



itlTT 



— e would then yield a phase factor e '^ . Thus, by varying 7, K could be any real number! 
So what is the reason for the above normalization of QlI From its definition and Stokes' 
theorem, we see immediately that adding one flux quantum 0o = ^ creates a 'charge' — ■§ at 
the boundary. On the other hand, we know from Laughlin's gauge argument that a charge 
—ev leaves the bulk when the flux quantum is added (and, in this framework, then circles 
along the boundary). But this means that K = v, 01 an = i^^, which is the classical result! 
Our feeling is that we find the same problem as in our discussion of Laughlin's argument in 
Section |2]^. Without localization, the classical result an = t-'\ should hold for any u. When 
will the value of an be robust against turning on some disorder? Presumably, this is exactly 
the case if v belongs to the set of fractions which yield an incompressible state. But no serious 
arguments are known for this hypothesis, apart from the fact that it fits well with numerical 
and experimental results. Thus to explain the existence of plateaux in the FQHE is still an 
open problem, but we believe that the reviewed results form an important step to its solution. 
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